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2.  PURPOSE:  To  compute  the  longitudinal  and  transverse  radii  of 
curvature  of  the  elliptic  cylinder  at  a  given  point. 

3.  METHOD:  The  longitudinal  radius  of  curvature  of  the  elliptic 
cylinder  (In  the  plane  of  Incidence)  at  the  point  defined  by 
elliptical  angle  VR  (as  shown  In  figure  1)  Is  given  by: 


=  (A2sin2VR  +  B2cos2VR)3/2 

pq  2 

*  AB  sin  a 

S 


The  transverse  radius  of  curvature  at  the  point  defined  by 
elliptical  angle  VR  Is  given  by: 

p  =  (A2sin2VR  ♦  B2cos2VR)3/2 

1  AB  sin2(a  -  ir/2) 

s 

where 

as  *  SAS 
Pg  "  I* 

Pt  *  RT. 


Y 


Figure  1.  Illustration  of  Elliptical  Cylinder  Geometry 
Used  In  Computing  the  Radii  of  Curvature  at 
Point  X 
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RADCV  (6TD) 

INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

A 

Cylinder  radius  along  the  x  axis 

B 

Cylinder  radius  along  the  y  axis 

R6 

Radius  of  curvature  in  the  plane  of 
incidence 

RGT 

Radius  of  curvature  of  the  elliptic 
cylinder  in  the  principal  x-y  plane 

RT 

Radius  of  curvature  transverse  to  the  plane 
of  Incidence 

SAS 

The  sine  of  AS,  where  AS  is  n  minus  THSR  (THSR 
is  the  theta  angle  of  the  observation 
direction  In  the  reference  coordinate 
system  (RCS)  relative-  to  the  cylinder  axis 
in  radians) 

SASP 

The  absolute  value  of  the  sine  of 
AS  -  ir/2,  where  AS  is  n  minus  THSR  (THSR  is 
the  theta  angle  of  the  observation 
direction  In  the  RCS  relative  to  the 
cylinder  axis  In  radians) 

VR 

Elliptic  angle  defining  the  desired  point 
on  cylinder 

I/O  VARIABLES: 

A. 

INPUT 

LOCATION 

A 

/GEOMEL/ 

B 

/GEOMEL/ 

SAS 

/GTD/ 

SASP 

/GTD/ 

VR 

F.P. 

B. 

OUTPUT 

LOCATION 

K*  «* 


NAME:  RCLDPL  (GJO) 

PURPOSE:  To  compute  the  unobstructed  electric  field  from  a  unit 
source  reflected  by  a  cylinder  and  diffracted  by  edge  ME  of  plate  MP 
Into  the  given  far-field  observation  direction  or  to  a  given  near¬ 
field  point. 

METHOD:  RCLDPL  is  the  driver  routine  which  directs  all  the  ray 
tracing,  physics  and  field  calculations  for  determining  the  electric 
field  from  a  unit  source  reflected  from  a  cylinder  and  then 
diffracted  by  a  plate  in  a  given  far-field  direction  or  to  a  given 
near-field  observation  point.  Pertinent  geometry  is  shown  in 
figure  1. 
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Figure  1.  Ray  Reflected  by  Cylinder  and  then  Diffracted 
by  Plate  Edge 


The  code  first  checks  the  wedge  angle  number  of  edge  ME  of  plate  MP. 
If  It  Is  greater  than  2,  indicating  It  is  part  of  a  wedge,  and  the 
edge  has  already  been  considered,  the  fields  are  set  to  zero.  If 


RCLDPL 


(GTO) 


debug  information  is  requested,  it  is  printed  on  file  LUPRNT.  Then 
control  is  returned  to  the  calling  routine.  If  the  wedge  angle 
number  Is  less  than  2,  the  code  then  checks,  for  far  field  only,  if 
diffraction  Is  possible.  If  it  is  not  possible,  a  flag  is  set  which 
Indicates  that  reflection  point  starting  data  are  not  available  for 
the  next  time  RCLDPL  is  called.  The  field  is  set  to  zero.  Debug 
information  (if  requested)  is  printed  on  file  LUPRNT,  and  control  is 
then  returned  to  the  calling  routine.  This  check  is  not  made  for 
near  field  at  this  point  in  the  code.  Now  for  both  near  field  and 
far  field  subroutine  RFDFPT  is  called.  RFDFPT  computes  the  ray  path 
and  checks  for  near-field  calculations  if  diffraction  was  possible. 
After  returning  from  RFOFPT,  the  code  makes  three  checks  to  deter¬ 
mine  if  reflection  and  diffraction  points  are  legal.  The  first 
check  is  to  make  sure  that  the  reflection  satisfies  Snell's  Law. 
Then  the  reflection  point  is  checked  to  make  sure  it  is  on  the 
curved  surface  of  the  cylinder.  The  diffraction  point  is  checked  to 
make  sure  it  is  on  edge  ME  of  the  plate.  If  any  of  these  checks 
fails,  the  fields  are  set  to  zero,  debug  information  is  printed  if 
it  was  requested,  and  control  is  returned  to  the  calling  routine. 
If  reflection  and  diffraction  have  occurred  properly,  then  the 
complete  ray  path  is  checked  for  obstructions.  If  it  is  obstructed 
at  any  location,  the  fields  are  set  to  zero.  Debug  information  is 
printed  if  it  was  requested,  and  then  control  returns  to  the  calling 
routine.  If  a  ray  path  is  unobstructed,  the  field  computations  can 
begin. 

The  polarization  unit  vectors  for  the  rays  incident  and  diffracted 
on  the  plate  and  incident  and  reflected  at  the  cylinder  are 
computed.  The  source  field  pattern  factor  is  found  by  calling 
subroutine  SOURCE.  The  first  field  computed  is  that  which  is 
incident  on  the  cylinder.  It  is  computed  in  components  perpen¬ 
dicular  and  parallel  to  the  plane  of  incidence.  Then  the  cylinder 
reflected  field  is  computed.  Following  this,  the  field  incident  on 
the  plate  can  be  computed  In  parallel  and  perpendicular  components. 
The  caustic  distances  and  ray  spreading  factors  are  computed  for  the 
reflected-dlffracted  ray.  The  phase  factor  is  computed.  Diffrac¬ 
tion  coefficients  are  found  by  calling  subroutine  DW.  Now  the  total 
diffracted  field  can  be  computed  and  converted  to  theta  and  phi 
components  In  the  reference  coordinate  system  (RCS).  Subroutine 
XYZFLD  Is  called  to  compute  the  x,  y,  and  z  components  of  the  field 
and  to  accumulate  them  with  the  fields  from  other  interactions. 

If  debug  information  was  requested,  the  total  field  magnitude  Is 
computed.  The  total  field  magnitude,  and  theta  and  phi  components 
are  printed  on  file  LUPRNT.  Control  is  then  returned  to  the  calling 
routine. 
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INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

BO 

Diffracted  field  polarization  unit  vector 
parallel  to  edge 

BOP 

Incident  field  polarization  unit  vector 
parallel  to  edge 

DO 

Normal ization  constant  for  cylinder  tangent 
vector 

DH 

Edge  diffraction  coefficient  for  hard  field 
components 

DHIT 

Distance  from  source  to  hit  point  (from 
PLAINT) 

DOTP 

Test  parameter  used  to  determine  if  reflec¬ 
tion  is  legal 

DPH 

Slope  diffraction  coefficient  for  hard 
boundary  condition 

DPS 

Slope  diffraction  coefficient  for  soft 
boundary  condition 

OS 

Diffraction  coefficient  for  soft  field 
components 

DV 

Dot  product  of  edge  unit  vector  and 
diffracted  ray  propagation  direction 

EDPH 

Phi  component  of  diffracted  field  in  RCS 

EDPL 

Diffracted  field  component  parallel  to  edge 

EOPR 

Diffracted  field  component  perpendicular  to 
edge 

EDTH 

Theta  component  of  diffracted  field  in  RCS 

EF 

Theta  component  of  source  field  pattern 
factor 

EG 

Phi  component  of  source  field  pattern 
factor 
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EIPL 

EIPR 

EIX,EIY,EIZ 

ERX.ERY.ERZ 

EXPH 

FN 

FNP 

GAM 

LHIT 

LRDC 

ME 

MP 

PH 

PH  ICR 

PHO 

PS 

PSOR 


Component  of  field  incident  on  cylinder  (or 
plate)  parallel  to  plane  of  incidence  (or 
edge) 

Component  of  field  incident  on  cylinder  (or 
plate)  perpendicular  to  plane  of  incidence 
(or  edge) 

Source  pattern  factors  for  x,y,z  components 
of  incident  E-field 

X,Y,Z  components  of  cylinder  reflected 
field  in  RCS 


Complex  phase  and  spreading  factor 

Wedge  angle  number 

2ir  minus  the  wedge  angle 

Dot  product  of  vector  to  the  diffraction 
point  with  the  observation  unit  vector 

Set  true  if  ray  hits  plate  (from  PLAINT) 

Set  true  if  reflection  data  are  available 
from  previous  pattern  angle  (or  for  next 
pattern  angle  (when  leaving  routine)) 

Edge  on  plate  MP  where  diffraction  occurs 

Plate  where  diffraction  occurs 

Diffracted  field  polarization  unit  vector 
normal  to  edge 

Phi  component  of  field  incident  on  cylinder 
In  RCS 

Incident  field  polarization  unit  vector 
normal  to  edge 

Diffracted  ray  phi  angle  in  diffraction 
point  coordinate  system  in  degrees 

Incident  ray  phi  angle  in  diffraction  point 
coordinate  system 
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PSR 

RHI1 

RHI2 

RHIE 

RH01 

RH02 

SBO 

SMAG 

SNF 

SP 

THICR 

TPP 


Diffracted  ray  phi  angle  in  diffraction 
point  coordinate  system 

Caustic  distance  of  cylinder  reflected 
field  incident  on  edge  in  the  direction 
perpendicular  to  the  edge 

Caustic  distance  of  cylinder  reflected 
field  incident  on  edge  in  the  direction 
parallel  to  the  edge 

Edge  caustic  distance 

Ray  spreading  radius  at  cylinder  in  plane 
normal  to  plane  of  incidence 

Ray  spreading  radius  at  cylinder  in  plane 
of  incidence 

Sine  of  the  diffraction  angle 

Length  of  ray  from  reflection  point  on 
cylinder  to  source  and  distance  between 
reflection  and  diffraction  points 

Distance  between  diffraction  point  on  plate 
and  near-field  observation  point 

Distance  between  reflection  and  diffraction 
point 

Theta  component  of  incident  ray  direction 
on  cylinder  in  RCS 

Distance  parameter  for  edge-diffracted 
field 


UB  Unit  binormal  of  elliptic  cylinder  at  phi 

angle  at  which  reflection  occurs  (2-D) 

UIPPX,UIPPY,UIPPZ  X,Y,Z  components  of  incident  polarization 

unit  vector  parallel  to  plane  of  incidence 


UIPP.X,UIPRY,UIPRZ  X,Y,Z  components  of  incident/reflected 

polarization  unit  vector  perpendicular  to 
plane  of  incidence 


UN 


Unit  normal  of  elliptic  cylinder  at  phi 
angle  at  which  reflection  occurs  (2-D) 
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URPPX , URPPY , URPPZ 

X,Y,Z  components  of  reflected  polarization 
unit  vector  parallel  to  the  plane  of  inci¬ 
dence 

VI 

X,Y,Z  components  of  ray  propagation 
direction  of  ray  incident  on  diffraction 
point 

VIC 

X,Y,Z  components  of  ray  propagation  direc¬ 
tion  of  ray  incident  on  cylinder 

VR 

Elliptical  angle  defining  reflection  point 
on  cylinder  (2-D) 

XD 

X,Y,Z  components  of  diffraction  point  in 
RCS 

XDMAG 

Normalization  constant  for  vector  from  RCS 
origin  to  diffraction  point 

XDP 

Modified  diffraction  point  location  for 
shadowing  test 

XI  MAG 

Normalization  constant  for  vector  from  RCS 
origin  to  second  corner  on  edge  ME 

XMAG 

Normalization  constant  for  vector  from 
origin  to  first  corner  on  edge  ME 

XR 

X,Y,Z  components  of  reflection  point  on 
cylinder 

XRR 

Reflection  point  location  on  cylinder 

XSS 

I/O  VARIABLES: 

Source  location 

A.  INPUT 

LOCATION 

A 

/GEOMEL/ 

8 

/GEOMEL/ 

BCD 

/BNDRCL/ 

CTC 

/GEOMEL/ 

0 

/OIR/ 

820 
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OP 

OPR 

DT 

FLOPT 

FN 

LDEBUG 

LNRFLD 

LRDC 

LUPRNT 

ME 

MEP 

MP 

PHSR 

PI 

THSR 

TPI 

V 

VN 

VP 

VXS 

X 

XS 

ZC 

8.  OUTPUT 
EDPH 


/THPHUV/ 

/PIS/ 

/THPHUV/ 

/NEAR/ 

F.P. 

/TEST/ 

/NEAR/ 

/CLRDC/ 

/ADEBUG/ 

F.P. 

/GEOPLA/ 

F.P. 

/DIR/ 

/PIS/ 

/DIR/ 

/PIS/ 

/GEOPLA/ 

/GEOPLA/ 

/GEOPLA/ 

/SORINF/ 

/GEOPLA/ 

/SORINF/ 

/GEOMEL/ 

LOCATION 

F.P. 
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RCIDPL 


(GTD) 


EDTH 

F.P. 

CALLING  ROUTINE: 

GTDDRV 

CALLED  ROUTINES: 

ASSIGN 

RFDFPT 

BEXP 

SMAGNF 

BTAN2 

SOURCE 

CYLINT 

STATIN 

OH 

STATOT 

NANOB 

TPNFLO 

NFD 

WLKBCK 

PLAINT 

XYZFLO 

J 

> . 
.v.v. 
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COMPUTE  POLARIZATION 
UNIT  VECTORS  FOR  INCI¬ 
DENT  AND  DIFFRACTED 
FIELD  ON  PLATE 


COMPUTE  REFLECTED- 
OIFFRACTED  RAY 
CAUSTIC  DISTANCES 
AND  SPREADING 
FACTORS 


COMPUTE  SOURCE  FIELD 
PATTERN  FACTOR 
SOURCE 


COMPUTE  PHASE 
FACTOR 


COMPUTE  POLARIZATION 
UNIT  VECTORS  FOR 
INCIDENT  AND 
REFLECTED  FIELDS  AT 
CYLINOER 


COMPUTE  DIFFRACTION 
COEFFICIENTS 

OH 


CALCULATE  COMPO¬ 
NENTS  OF  FIELD 
INCIDENT  ON  CYLINOER 
PERPENDICULAR  AND 
PARALLEL  TO  PLANE 
OP  INCIDENCE 


COMPUTE  DIFFRACTED 
FIELDS 


COMPUTE  COMPONENTS 
OF  CYLINDER 
REFLECTED  FIELD 


COMPUTE  THETA  AND 
PHI  COMPONENTS 
OF  DIFFRACTED  FIELO 
INRCS 


COMPUTE  COMPONENTS 
OF  FIELD  INCIDENT  ON 
PLATE  EDGE  PARALLEL 
ANO  PERPENDICULAR  TO 
EDGE 


COMPUTE  x,  y.  *  COM¬ 
PONENTS  OF  FIELD 
ANO  ACCUMULATE 
XYZFLO 


PURPOSE:  To  compute  the  geometrical  optics  field  reflected  from  the 
elliptic  cylinder  and  then  reflected  by  plate  MP. 


METHOD:  RCLRPL  functions  as  a  service  routine  for  subroutine  SCLRPL 
where  the  total  cylinder-plate  scattered  fields  are  computed.  The 
field  components  computed  in  RCLRPL  which  are  used  in  SCLRPL  are  the 
hard  (EHTHJ,  EHPHJ)  and  soft  (ESTHJ,  ESPHJ)  theta  and  phi  components 
of  the  source  field  incident  on  the  cylinder  at  the  reflection 
point.  These  components,  along  with  several  other  useful  parameters, 
are  passed  to  subroutine  SCLRPL  through  common  block  /FUDGJ/. 

The  geometrical  optics  reflected  field  components  ETH  and  EPH  are 
computed  In  RCLRPL.  These  are  calculated  for  the  cylinder- 
reflected,  plate-ref lected  fields  from  a  unit  source  in  the  given  far- 
fleld  observation  direction  or  to  a  given  near-field  observation 
point.  These  components  are  not  used  presently.  The  pertinent 
geometry  for  this  routine  is  shown  in  figure  1. 
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CYLINDER  REFLECTION 
POINT  IMAGE 


OBSERVATION 

FIELD  POINT 

IMAGE  FLDPTl 

(FOR  NEAR  FIELD  ONLY! 


Figure  1.  Illustration  of  Ray  Reflected  by  Cylinder 
and  Then  Reflected  by  Plate. 

The  code  first  determines  the  ray  path  cylinder  and  plate  reflection 
points.  The  procedure  followed  Is  different  for  near-field  and  far- 
fleld  calculations.  The  flowchart  shows  near-field  and  far-field 
paths. 
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The  reflection  point  on  the  cylinder  is  found  by  imaging  the 
observation  direction  for  far  field  or  the  observation  point  for 
near  field  through  plate  HP.  This  point  Is  checked  to  make  sure  it 
lies  on  the  cylinder  within  the  end  caps.  If  it  does  not,  the  field 
is  set  to  zero  and  control  Is  returned  to  the  calling  routine.  If 
It  Is  a  legal  point,  the  reflection  point  on  the  plate  is  checked  to 
make  sure  it  is  a  legal  point  also.  If  plate  reflection  did  not 
occur,  the  field  Is  set  to  zero  and  control  is  returned  to  the 
calling  routine.  If  the  point  Is  legitimate,  the  ray  path  from  the 
source  to  the  cylinder,  to  the  plate,  to  the  far-field  observation 
direction  or  near-field  observation  point  is  checked  for  shadowing. 
If  the  path  Is  shadowed  anywnere,  again  the  field  is  set  to  zero  and 
control  returns  to  the  calling  routine.  If  the  path  is  clear,  then 
at  this  point  It  Is  known  that  the  cylinder-reflected,  plate- 
reflected  field  does  exist. 


The  physics  and  field  computations  begin  by  computing  the  source 
field  pattern  factor  from  the  source  by  calling  subroutine  SOURCE. 
Then  the  spreading  radii  needed  FOR  including  the  effect  that  the 
curved  cylinder  wall  has  on  spreading  the  cylinder-reflected  wave  is 
computed.  Other  parameters,  polarization  vectors,  and  Image  loca¬ 
tions  are  calculated  for  the  field  computations.  The  code  computes 
the  fields  Incident  on  the  cylinder,  the  cylinder-reflected  field 
and  the  plate-reflected  field.  The  plate-reflected  field's  phase 
factor  Is  based  on  the  cylinder  reflection  point  imaged  through 
plate  HP.  For  far  field,  this  refers  the  field  to  the  origin  of 
the  reference  coordinate  system.  For  near  field,  the  phase  factor 
includes  the  spherical  wave  spreading  factor.  The  code  ends  by 
computing  the  hard  (theta  and  phi)  and  soft  (theta  and  phi)  compo¬ 
nents  of  the  field  Incident  on  the  cylinder. 


4.  INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 


Al,  A2 


Field  components  of  ray  Incident  on  plate, 
normal  and  tangent  to  the  plate 


A3 


Determinant  of  polarization  transformation 


CI1,C12,C21,C22  Coefficients  used  to  convert  polarization 

from  theta  and  phi  components  in  RCS  to 
components  normal  and  tangent  to  plate  (and 
vice-versa) 

CTHW  2-D  dot  product  of  unit  normal  at  cylinder 

reflection  point  and  ray  propagation  direc¬ 
tion  between  reflection  points 
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Propagation  direction  after  plate  reflec¬ 
tion  in  x,y,z  RCS  components 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  1  (2-D) 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  2  (2-D) 

Distance  from  source  to  hit  point  (from 
PLAINT) 

Distance  between  cylinder  and  plate  reflec¬ 
tion  point  (from  subroutine  PLAINT) 

Distance  to  hit  point  (from  PLAINT  and 
CYLINT) 

X,Y,  and  Z  components  of  incident  ray 
direction  on  cylinder  in  RCS 

X,Y,Z  components  of  propagation  direction 
of  ray  incident  on  plate 

Oistance  between  plate  reflection  point  and 
near-field  observation  point 

Test  variable 

Phi  unit  vector  for  observation  direction  D 

Theta  unit  vector  for  observation  direc¬ 
tion  D 

Pattern  factor  of  theta  component  of  inci¬ 
dent  field  in  RCS  (also  theta  component  of 
cylinder-reflected  field  In  RCS) 

Pattern  factor  of  phi  component  of  Incident 
field  in  RCS  (also  phi  component  of  cylin¬ 
der-reflected  field  in  RCS) 

Phi  component  of  hard  component  of  field 
Incident  on  cylinder 

Theta  component  of  hard  component  of  field 
incident  on  cylinder  (parallel  to  plane  of 
incedence) 
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EIPP 

Incident  field  component  parallel  to  plane 
of  incidence  on  cylinder 

EIPR 

Incident  field  component  perpendicular  to 
plane  of  incidence  on  cylinder 

EPH 

Phi  component  of  total  cylinder-reflected, 
plate-ref lected  field 

ERPP 

Component  of  cylinder-reflected  field 

parallel  to  plane  of  incidence 

ERPR 

Component  of  cylinder-reflected  field  per¬ 
pendicular  to  plane  of  incidence 

ERX,  ERY,  ERZ 

X,Y,Z  components  of  cylinder-reflected 
field  in  RCS 

ESPHJ 

Phi  component  of  soft  component  of  field 
incident  on  cylinder 

ESTHJ 

Theta  component  of  soft  component  of  field 
incident  on  cylinder 

ETH 

Theta  component  of  total  cylinder-reflec¬ 
ted,  plate-reflected  field 

EX,  EY,  EZ 

X.Y.Z  components  of  source  field  pattern 
factor  in  RCS 

FLDPTI 

X.Y.Z  components  of  the  location  of  the 
near-field  observation  point  image  through 
plate  MP 

GAM 

Phase  constant 

LHIT 

Set  true  if  ray  hits  plate  (from  PLAINT) 

LRFS 

Set  true  if  reflection  data  is  available 
from  previous  pattern  angle  (or  for  next 
pattern  angle  (when  leaving  routine)) 

LTRFJ 

Set  true  if  geometrical  optics  reflected- 
reflected  fields  do  not  exist 

MP 

Plate  on  which  reflection  occurs  after 
cylinder  reflection 

PH 

Complex  phase  constant 

RH01J 


SMAGJ 


SXN,  SYN,  SZN 


Phi  component  of  incident  ray  direction  on 
cylinder  in  RCS 

Phi  component  of  ray  propagation  direction 
between  cylinder  and  plate  in  RCS 

Radius  of  curvature  of  cylinder  at  reflec¬ 
tion  point 

Ray  spreading  radius  in  plane  of  cylinder 
curvature  at  reflection  point 

Ray  spreading  radius  in  plane  normal  to 
plane  of  incidence  at  reflection  point 

Distance  between  reflection  points  on 
cylinder  and  plate 

Distance  between  cylinder  reflection  point 
and  near-field  observation  point  image 
through  plate  MP  (therefore  distance  of 
complete  ray  path  between  reflection  point 
on  cylinder  and  the  near-field  observation 
point) 

Length  of  ray  from  reflection  point  on 
cylinder  to  source 

Distance  from  plate  reflection  point  to 
near-field  observation  point 

X,Y,Z  components  of  unit  vector  of  ray 
from  reflection  point  on  cylinder  to  source 
location  in  RCS 

Theta  component  of  incident  ray  direction 
on  cylinder 

Theta  component  of  ray  propagation  direc¬ 
tion  between  cylinder  and  plate 

Unit  binormal  at  the  cylinder  reflection 
point 


UIPPX.UIPPY.UIPPZ 


X,Y,Z  components  of  incident  polarization 
unit  vector  parallel  to  plane  of  incidence 
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UIPRX.UIPRY.UIPRZ 

UN 

UR 

URPPX.URPPY.URPPZ 

VR 

VT 

VXS 

XRJ 

XRR 

XRS 

XSS 

5.  I/O  VARIABLES: 

A.  INPUT 


X,Y,Z  components  of  incident/reflected 
polarization  unit  vector  perpendicular  to 
plane  of  incidence 

Unit  normal  at  the  cylinder  reflection 
point 

The  z  component  of  the  location  of  the 
reflection  point  on  the  cylinder 

X,Y.Z  components  of  reflected  polarization 
unit  vector  parallel  to  plane  of  incidence 

Phi  angle  used  to  define  the  x  and  y  com¬ 
ponents  of  the  reflection  point  on  cylinder 

X,  Y,  Z  components  of  polarization  unit 
vector  tangent  to  plate  and  normal  to  ray 
incident  on  plate 

Matrix  defining  source  coordinate  system 
axes  in  RCS  components 

X,  Y,  Z  components  of  reflection  point 
location  on  cylinder 

Cylinder  reflection  point  location 

Reflection  point  on  plate  (also  cylinder 
reflection  point  image  location  in  plate) 
Also  cylinder  reflection  point 

Source  location 


LOCATION 


A 

B 

BTS 

CTC 

D 


/GEOMEL/ 

/GEOMEL/ 

/BNDSCL/ 

/GEOMEL/ 

/DIR/ 

/THPHUV/ 


DP 


RCLRPL 


OT 

/THPHUV/ 

DTS 

/BNDSCL/ 

FLDPT 

/NEAR/ 

LNRFLD 

/NEAR/ 

LRFS 

/CLRFS/ 

MP 

F.P. 

PHSR 

/DIR/ 

PI 

/PIS/ 

THSR 

/DIR/ 

TP  I 

/PIS/ 

VN 

/GEOPLA/ 

VXS 

/SORINF/ 

XS 

/SORINF/ 

ZC 

/GEOMEL/ 

OUTPUT 

LOCATION 

EHPHJ 

/FUDGJ/ 

EHTHJ 

/FUDGJ/ 

EPH 

F.P. 

ESPHO 

/FUDGJ/ 

ESTHJ 

/FUDGJ/ 

ETH 

F.P. 

LRFS 

/CLRFS/ 

LTRFJ 

/FUDGJ/ 

RGJ 

/FUGDJ/ 

RHOIJ 

/FUDGJ/ 

SMAGJ 

/FUDGJ/ 

SNFF 

/DIST/ 

TRANJ 

/FUDGJ/ 

XRJ 

/FUDGJ/ 

CALLING  ROUTINE: 

SCLRPL 

CALLED  ROUTINES: 

ASSIGN 

NFD 

SOURCE 

BEXP 

PLAINT 

STATIN 

BTAN2 

REFBP 

STATOT 

CYLINT 


RFDFIN 


TPNFLD 
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— RAY  PATH  T0\ 
REFLECTION  POINT 
ON  PLATE 

^V^SHORTEST^^ 


-^RAY  PATH\. 
•''FROM  CYLINDER^ 
TO  PLATE  TO 
OBSERVATION 
DIRECTION  OR  POINT^ 
^^HADOWED/^ 


-KD 


COMPUTE  POLARIZA 
TION  UNIT  VECTORS 
FOR  INCIDENT  AND 
REFLECTED  FIELD  AT 
CYLINDER 


COMPUTE  FIELD  COM 
PONENTS  OF  RAY  INCI 
DENT  ON  CYLINDER 
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COMPUTE  DIRECTION 
OF  RAY  INCIDENT 
ONTO  CYLINDER 


RAY  FROM*' 
SOURCE  TO 
CYLINDER 
.SHADOWS  Dx 


COMPUTE  SOURCE 
FIELD  PATTERN  FAC¬ 
TOR 

SOURCE 


COMPUTE  SPREADING 
RADII  FOR  CYLINDER 
REFLECTION  POINT 


COMPUTE  LOCATION 
OF  CYLINDER  REFLEC¬ 
TION  POINT  IMAGE  IN 
PLATE  MP 

IMAGE 


COMPUTE  PHASE  FAC 
TOR  FOR  CYLINDER- 
REFLECTED.  PLATE  - 
REFLECTED  WAVE 


COMPUTE  COMPONENTS 
OF  CYLINDER-RE 
FLECTED  FIELD  PARA 
LLEL  AND  PERPEN¬ 
DICULAR  TO  PLANE 
OF  INCIDENCE _ 


COMPUTE  THETA  AND 
PHI  COMPONENTS  OF 
CYLINDER-REFLECTED 
FIELD 


CALCULATE  POLAR¬ 
IZATION  VECTORS  AND 
DOT  PRODUCTS  NECES¬ 
SARY  TO  COMPUTE 
PLATE-REFLECTED 
FIELD 


COMPUTE  FIELD 
REFLECTED  FROM 
PLATE 


MOTS' VVW1VA\.V  V  V.VA',%  A  /-  .•<  > 


LTRFJ-.TRUE. 


SET  FIELD  TO  ZERO 


1.  NAME:  RDEFIL 


(GTO,  INPUT,  MOM,  OUTPUT) 

2.  PURPOSE:  Read  data  from  the  logical  unit  specified  and  increment 
the  internal  file  pointers  to  Indicate  the  current  file  position. 

3.  METHOO:  The  number  of  words  between  the  current  file  position  and 
the  end  of  file  Is  determined  and,  if  less  than  the  number  of  words 
requested,  a  fatal  error  Is  generated.  Otherwise,  the  file  is  read 
in  the  binary  mode  and  the  current  file  pointer  is  incremented  to 
point  at  the  last  word  read  from  the  file. 

4.  INTERNAL  VARIABLES: 


VARIABLE 

LUNIT 

NUMLFT 

NWORDS 

X WORDS 

5.  I/O  VARIABLES: 

A.  INPUT 
DBGPRT 
IOCKPT 
IOFILE 
ISON 
LUNIT 
LUPRNT 
MODCHK 
NDFILE 
NWORDS 


DEFINITION 

Input  argument  designating  logical  unit  to 
read 

Number  of  words  left  before  the  end  of  the 
file  from  current  pointer  position 

Input  argument  designating  number  of  words 
to  be  read 

Input  array  into  which  the  data  will  be 
read 


LOCATION 

/ADEBUG/ 

/SYSFIL/ 

/IOFLES/ 

/ADEBUG/ 

F.P. 

/AOEBUG/ 

/SYSFIL/ 

/IOFLES/ 

F.P. 


P  RDEFIL 

•*.4 

(GTD,  INPUT,  MOM,  OUTPUT) 

•  « 

B.  OUTPUT 

LOCATION 

1  IERRF 

/ADEBUG/ 

v:  IOFILE 

/ IOFLES/ 

V-  X WORDS 

F.P. 

1  6.  CALLING  ROUTINES*: 

y  BUBBLE  (1) 

RESTRT  (1) 

DECOMP  (3) 

RWCOMS  (1,2, 3, 4) 

P  GEODRV  (1) 

RWFILS  (1,2, 3,4) 

|  GETSYM  (1,2, 3, 4) 

SOLDRV  (3) 

jgj  MOVFIL  (1,2, 3, 4) 

STRTUP  (2,3,4) 

I  PUTSYM  (1, 2,3,4) 

SUBPAT  (1) 

js]  7.  CALLED  ROUTINES: 

I  ASSIGN 

1  ERROR 

g  STATIN 

STATOT 

|  WLKBCK 

1  *1- INPUT 

S  2-GTD 

ffl  3-MOM 

JJ  4-OUTPUT 

840 

RDEFIL 


(GTD,  INPUT,  MOM,  OUTPUT) 
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PURPOSE:  To  reblock  the  interaction  matrix  Into  several  square 
submatrlces  when  structure  symmetry  Is  present. 

METHOD:  When  structure  symmetry  Is  present  the  full  square  inter¬ 
action  matrix  Is  not  generated.  Instead,  only  an  NR  x  NC  matrix  is 
needed,  where  NC  is  the  number  of  elements  per  symmetry  cell,  and  NR 
is  the  total  number  of  elements.  NR  is  always  an  integer  multiple 
of  NC,  this  integer  being  the  number  of  symmetry  cells.  Since  the 
matrix  problem  will  be  solved  In  (NC  x  NC)  blocks,  the  data  must  be 
reblocked  Into  that  format. 

Each  column  of  NC  elements  is  read  into  core  from  the  input  symbol 
(the  matrix  is  stored  in  transposed  form),  and  the  proper  elements 
stored  In  the  columns  of  each  submatrix  of  the  output  data  set,  as 
shown  in  figure  1. 


NC  NC  NC 


NR/NC  =  3  SYMMETRY  CELLS 


Figure  1.  Illustrating  the  Reblocklng  of  an  NRxNC  Matrix 
Into  an  NCxNR  Matrix 


INTERNAL  VARIABLES: 


VARIABLE 


DEFINITION 

Internal  variable  equal  to  NC 

Pointer  to  column  number  being  read  from 
input  data  set 

Pointer  to  column  number  being  written  to 
output  data  set 


PREVIOUS  PAGE 
IS  BLANK 


LOCNAM 


Pointer  to  input  data  set  name  in  symbol 
table 


MORE 

N 

NAMEZ 

NAMEZ1 

NBIAS 

NBITS 

NC 

NPRELM 

NR 

NUMMAT 

Z 

I/O  VARIABLES: 

A.  INPUT 
ISON 
KBCPLX 
KOLB IT 
KOLNAM 
LUPRNT 
NAMEZ 
NAMEZ1 


Flag  indicating  a  complex  data  set 

Loop  index  over  symbol  table  entries  and 
submatrices 

User-assigned  name  of  input  data  set  (NR  x 
NC) 

User-assigned  name  of  output  data  set  (NC  x 
NR) 

Pointer  to  beginning  of  output  for  a  column 
of  NAMEZ1 

Attribute  word  of  input  data  set 

Number  of  columns  of  input  data  set 

Number  of  data  words  per  matrix  element 

Number  of  rows  of  input  data  set 

Number  of  symmetry  cells  in  input  data  set; 
number  of  submatrices  in  output  data  set 

Temporary  storage  for  matrix  reblocking 

LOCATION 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

F.P. 

F.P. 


NDATBL 


6. 


7. 


NPDATA 

NR 

Z 

B.  OUTPUT 
IERRF 
Z 

CALLING  ROUTINE: 
ZIJDRV 

CALLED  ROUTINES: 


/PARTA8/ 

/PARTAB/ 

F.P. 

F.P. 

LOCATION 

/ADEBUG/ 

F.P. 


ASSIGN 

GETSYM 

STATIN 

CONVRT 

IBITCK 

STATOT 

ERROR 

PUTSYM 

WLKBCK 
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PURPOSE:  To  calculate  the  incident  ray  direction  needed  in  order  to 
obtain  the  reflected  ray  in  a  given  direction  from  a  specified 
plate. 


METHOD:  The  incident  ray  unit  vector  (VI)  is  found  by  imaging  the 
reflected  ray  unit  vector  (DR)  through  the  plate  (MP).  Figure  1 
shows  the  important  geometry.  The  equation  for  VI  is: 

A  A  AAA 

VI  =  DR  -  2 ( VN  •  DR)  VN 


A 

The  theta  and  phi  angles  which  define  VI  are  sent  back  to  the 
calling  routine. 


PLATE 

NORMAL 


Figure  1.  Illustration  of  Incident  and  Reflected 
Rays  on  Plate 


INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

CPS 

Cosine  of  PHSR 

CTS 

Cosine  of  THSR 

ON 

Cross  product  of  DR  and  VN 

DR 

Reflected  ray  propagation  direction 
x,y,z  RCS  components 

in 

ERD 

Error  detection  variable 

LUPRNT 

File  on  which  warning  message  will 
printed 

be 

MP 

Plate  upon  which  reflection  occurs 

PHIR 

Phi  component  of  incident  ray 
direction  in  RCS 

propagation 

PHSR 

Phi  component  of  reflected  ray 
direction  in  RCS 

propagation 

SPS 

Sine  of  PHSR 

STS 

Sine  of  THSR 

THIR 

Theta  component  of  incident  ray 
direction  in  RCS 

propagation 

THSR 

Theta  component  of  reflected 
gation  direction  in  RCS 

ray  propa- 

VI 

XtY,Z  components  of  incident 
gation  direction  in  RCS 

ray  propa- 

VIN 

Dot  product  of  plate  normal  and 

VI 

VN 

Array  which  includes  unit  vector  normal  to 
plate  MP 

I/O  VARIABLES: 

A. 

INPUT 

LOCATION 

LUPRNT 

/ADEBUG/ 

MP 

F.P. 

PHSR 

F.P. 

THSR 

F.P. 

VN 

/GEOPLA/ 

9.  OUTPUT 


LOCATION 


1.  NAME:  REFCAP 


(GTD) 


m 


PURPOSE:  To  calculate  the  unobstructed  electric  field  resulting 
from  the  reflection  of  a  unit  source  off  a  given  cylinder  end  cap. 

METHOD:  REFCAP  is  the  driver  routine  which  directs  all  the  ray 
tracing,  physics  and  field  calculations  for  determining  the  electric 
field  reflected  by  an  elliptical  cylinder  end  cap  in  a  given  far- 
field  direction  or  to  a  given  near-field  observation  point  from  a 
unit  source.  The  important  geometry  is  shown  in  figure  1. 


OBSERVATION 
DIRECTION  FOR 
FAR  FIELD  OR  TO 
NEAR-FIELD 
OBSERVATION  POINT 


SOURCE  _ 

LOCATION  XS 


SOURCE  IMAGE 
LOCATION  XIC 


Figure  1.  Illustration  of  Source  Ray  Reflection  from  End  Cap 


First,  the  ray  from  the  source  Image  location  is  checked  to  make 
sure  It  passes  through  the  given  cylinder  end  cap.  If  it  does  not, 
the  theta  and  phi  components  of  the  field  are  set  to  zero,  and  no 
other  computations  except  debug  functions  (If  requested)  are  per¬ 
formed  In  this  routine.  If  reflection  can  occur,  the  ray  path  from 
the  reflection  point  in  the  far-field  observation  direction  or  to 
the  near-field  observation  point  is  checked  for  obstructions.  If  It 
Is  blocked  by  a  plate,  the  field  components  are  set  to  zero,  and  no 
other  computations  except  debug  functions  (if  requested)  are 
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REFCAP 
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performed  in  this  routine.  If  the  path  is  unobstructed,  the  ray 
path  from  the  source  location  to  the  reflection  point  is  checked  for 
obstructions.  If  a  plate  blocks  this  path,  the  theta  and  phi  field 
components  are  set  to  zero,  and  no  other  computations  except  debug 
functions  (if  requested)  are  performed  in  this  routine.  If  this 
pa*h  is  clear,  it  is  then  known  that  reflection  off  the  given  end 
cap  did  occur  and  that  the  complete  ray  path  is  unobstructed. 

The  source  field  pattern  factor  from  the  source  image  location  is 
computed  by  calling  subroutine  SOURCE  and  multiplying  the  returned 
field  values  by  the  reflection  coefficient.  Then  the  phase  factor 
is  computed.  For  far  field  this  factor  refers  the  field  to  the 
reference  coordinate  system  (RCS)  origin.  For  near  field,  the  phase 
factor  includes  the  spherical  wave  spread  factor.  Now  the  theta  and 
phi  components  of  the  field  can  be  computed.  The  electric  field  is 
given  by: 


The  x,  y,  z  components  of  the  field  are  then  computed  and  added  to 
the  previous  components  due  to  other  reflection-diffraction  interac¬ 
tions.  The  values  are  stored  in  common  block  /FLDXYZ/. 


If  the  debug  capabilities  have  been  requested,  the  end  cap  reflected 
field  magnitude  Is  computed.  The  magnitude,  theta  and  phi  complex 
components  are  printed  on  file  LUPRNT. 


REFCAP 


(STD) 


INTERNAL  VARIABLES: 

VARIABLE 

D 

OHIT 

DHT 

DI 

DN 

DNI 

EF 

EG 

EIX 

EIY 

EIZ 

EPH 

ETH 

EX 

FLDMAG 

FLDPT 

GAM 


DEFINITION 

Observation  direction  unit  vector 

Distance  from  source  to  hit  point  on  plate 
(from  PLAINT) 

Distance  from  source  to  hit  point  on  end 
cap  (from  CAP I NT) 

Unit  vector  of  incident  ray  propagation 
direction 

Dot  product  of  reflected  ray  propagation 
direction  and  end  cap  unit  normal 

Dot  product  of  incident  ray  and  end  cap 
unit  normal 

Pattern  factor  for  theta  component  of 

incident  E-field 

Pattern  factor  for  phi  component  of 

incident  E-field 

X  component  of  source  pattern  factor  of 
incident  E-field 

Y  component  of  source  pattern  factor  of 
incident  E-field 

Z  component  of  source  pattern  factor  of 
incident  E-field 

Phi  component  of  reflected  E-field  in  RCS 

Theta  component  of  reflected  E-field  in  RCS 

Phase  term  and  for  near  field  It  also 
contains  the  spherical  wave  spread  factor 

The  electric  field  magnitude 

The  near-field  observation  point  in  x,y,z 
components 

Phase  term  parameter 
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LHIT 

MC 

N 

NC 

HI 

NJ 

PHSR 

SNF 

THSR 

VAX 

VN 


XIS 

XS 

XSS 

I/O  VARIABLES: 
A.  INPUT 
CNC 
D 

FLDPT 

LDEBUG 

LNRFLD 


Set  true  if  ray  hits  plate  (from  PLAINT) 

End  cap  where  reflection  occurs 

DO  loop  variable 

Sign  change  variable 

DO  loop  variable 

00  loop  variable 

Observation  direction  phi  angle 

Distance  from  field  observation  point  to 
source  image  location 

Observation  direction  theta  angle 

X,Y,Z  components  defining  the  image  source 
coordinate  system  in  x,y,z  RCS  components 

Unit  normal  to  end  cap  in  RCS  x,y,z 
components 

Source  image  location 

Source  location  in  x,y,z  components 

Source  location  in  x,y,z  components  in  RCS 


LOCATION 

/GEOMEL/ 

/DIR/ 

/NEAR/ 

/TEST/ 

/NEAR/ 

/ADEBUG/ 


LUPRNT 


REFCAP  (GTD) 


PHSR 

/DIR/ 

SNC 

/GEOMEL/ 

THSR 

/DIR/ 

TP  I 

/PIS/ 

VXIC 

/IMCINF/ 

XIC 

/IMCINF/ 

XS 

/SORINF/ 

B.  OUTPUT 

LOCATION 

EPH 

F.P. 

ETH 

F.P. 

CALLING  ROUTINE: 

GTDORV 

CALLED  ROUTINES: 

ASSIGN 

BEXP 

CAP I NT 

NFD 

PLAINT 

SMAGNF 

SOURCE 

STATIN 

STATOT 

Wl KBCK 

XYZFLO 


855 


'£Jt 
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SET  FIELDS 
TO  ZERO 


COMPUTE  FIELD 
MAGNITUDE 


PRINT  FIELD 
PARAMETERS 
ON  FILE  LUPRNT 


w. 


NAME:  REFCYL 


(GTD) 


PURPOSE:  To  compute  the  geometrical  optics  field  due  to  reflection 
of  a  unit  source  off  the  curved  surface  of  the  elliptical  cylinder. 

METHOD:  REFCYL  functions  as  a  service  routine  for  subroutine 
SCTCYL,  where  the  total  cylinder-scattered  field  is  computed.  The 
field  components  computed  In  REFCYL  which  are  used  in  SCTCYL  are  the 
hard  (EHTH,  EHPH)  and  soft  (ESTH,  ESPH)  theta  and  phi  components  of 
the  field  incident  on  the  cylinder  at  the  cylinder  reflection  point. 
These  components  along  with  several  other  useful  parameters  are 
passed  to  subroutine  SCTCYL  through  common  block  /FUDG/.  The  geome¬ 
trical  optics  reflected  field  components  ETH  and  EPH  are  computed  in 
REFCYL.  These  are  calculated  for  the  cylinder-reflected  field  from 
a  unit  source  in  a  given  far-fleld  observation  direction  or  to  a 
given  near-field  observation  point.  These  components  are  not  used 
presently.  Pertinent  geometry  is  shown  in  figure  1. 


«  SOURCE 
LOCATION 


Figure  1.  Geometry  of  Ray  Reflected  From  Cylinder 


The  code  first  checks  to  see  If  reflection  is  possible.  If  it  Is 
not  possible  the  code  sets  flags  to  Indicate  that  starting  point 
data  are  not  available  for  the  next  time  REFCYL  is  called  and  that 
reflection  did  not  occur.  The  fields  are  set  to  zero  and  control  Is 
returned  to  the  calling  routine.  If  reflection  is  possible,  the 
reflection  point  is  computed  by  calling  subroutine  RFDFIN  for  near¬ 
field  calculations  and  by  calling  subroutine  RFPTCL  for  far-field 
calculations.  Then  the  far-field  reflection  point  is  checked  to 
make  sure  It  satisfies  the  law  of  reflection.  This  is  based  on  a 
returned  value  from  RFPTCL.  Then  for  both  far-field  and  near-field 
cases,  the  reflection  point  is  checked  to  make  sure  it  lies  on  the 
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curved  sides  of  the  cylinder  within  the  end  cap  boundaries.  Next 
the  ray  path  is  checked  to  make  sure  it  is  not  shadowed.  If  it  is 
shadowed,  the  code  sets  the  flag  which  indicates  a  reflected  field 
did  not  occur  and  sets  the  field  to  zero.  Control  is  then  returned 
to  the  calling  routine.  If  the  ray  path  is  unobstructed,  then  at 
this  point  it  is  known  that  reflection  did  occur  and  the  fields  can 
be  computed. 

The  physics  and  field  calculations  begin  by  computing  the  source 
field  pattern  factor  from  the  source  and  computing  the  cylinder- 
reflected  wave  spreading  radii.  Other  parameters  and  polarization 
vectors  are  then  calculated  for  the  field  computations.  (See  the 
flowchart).  The  code  computes  the  field  incident  on  the  cylinder 
and  the  field  reflected  from  the  cylinder.  The  phase  factor  refers 
the  far-field  cylinder-reflected  field  to  the  origin  of  the 
reference  coordinate  system  (RCS),  and  includes  for  near-field 

calculations  the  spherical  wave  spreading  factor.  The  code  ends  by 
computing  the  hard  (theta  and  phi)  and  soft  (theta  and  phi)  compo¬ 
nents  of  the  field  incident  on  the  cylinder  at  the  reflection  point. 

Additional,  in-depth  details  to  this  solution  are  given  on  pages 

105-107  of  reference  A. 

INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

CSV  Cosine  of  VR 

CTHI  Dot  product  of  cylinder  normal  and  reflec¬ 

tion  propagation  direction  unit  vector 

CW  Cosine  of  WR 

D  Propagation  direction  after  reflection  in 

x,y,z  RCS  components 

D1  Variable  used  in  various  near-field  calcu¬ 
lations  as  a  unit  vector  to  indicate  the 

direction  between  two  points 

D2  Variable  used  in  various  near-field  calcu¬ 
lations  as  a  unit  vector  to  indicate  the 

direction  between  two  points 


D12 


Dot  product  of  source  vectors  tangent  to 
cylinder  (2-0) 
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DICOEF 


Normalization  constant  for  reflection  point 
unit  normal  (from  RFPTCL) 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  1  (2-D) 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  2  (2-D) 

Distance  from  source  to  hit  point  (from 
PLAINT) 

X,Y,  and  Z  components  of  incident  ray 
direction  in  RCS 

Difference  of  dot  products  returned  from 
subroutine  RFPTCL  (2-D) 

The  phi  unit  vector  for  observation  direc¬ 
tion  D 

The  theta  unit  vector  for  observation 
direction  D 

Dot  product  of  vector  from  origin  to  source 
and  propagation  direction  (2-D) 

Pattern  factor  of  theta  component  of 
incident  field  in  RCS 

Pattern  factor  of  phi  component  of  incident 
field  in  RCS 

Phi  component  of  the  hard  component  of 
field  incident  on  cylinder 

Theta  component  of  the  hard  component  of 
field  incident  on  cylinder 

Incident  field  component  parallel  to  plane 
of  Incidence 

Incident  field  component  perpendicular  to 
plane  of  incidence 

Phi  component  of  reflected  E-field 


ERPP 

Reflected  field  component  parallel  to  plane 
of  incidence 

ERPR 

Reflected  field  component  perpendicular  to 
plane  of  incidence 

ERX.ERY.ERZ 

X,Y,Z  components  of  reflected  field  in  RCS 
(Also  used  to  define  components  incident  on 
cylinder) 

ESPH 

Phi  component  of  the  soft  component  of 
field  incident  on  cylinder 

ESTH 

Theta  component  of  the  soft  component  of 
field  incident  on  cylinder 

ETH 

Theta  component  of  reflected  E-field 

EX,EY,EZ 

Pattern  factor  of  x,y,z  components  of  inci¬ 
dent  field  in  RCS 

FPTXY 

The  equivalent  near-field  observation  point 
in  the  x-y  plane 

LHIT 

Set  true  if  ray  hits  plate  (from  PLAINT) 

LRFC 

Set  true  if  reflection  data  are  available 
from  previous  pattern  angle  (or  for  next 
pattern  angle  when  leaving  routine) 

LTRF 

Set  true  if  geometrical  optics  reflected 
field  does  not  exist 

ORIGIN 

The  x,y,z  components  of  the  origin  of  the 
reference  coordinate  system  (0.,  0.,  0.) 

PH 

Phase  and  magnitude  constant  for  incident 
or  reflected  field 

PHIR 

Phi  component  of  incident  ray  direction 

PHSR1 

Variable  used  in  various  near-field  calcu¬ 
lations  to  indicate  the  phi  angle  between 
two  points 

PHSR2 

Variable  used  in  various  near-field  calcu¬ 
lations  to  indicate  the  phi  angle  between 
two  points 
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RG 

RH01 

RH02 

S 

SMAG 

SNF 

SNV 

SNX 

SNY 

SQRH 

SW 

SXN.SYN 

THIR 

THSR1 

THSR2 

TRAN 

TX1 

TX2 


Parameter  used  in  transition  function 

Ray  spreading  radius  in  plane  of  cylinder 
curvature  at  reflection  point 

Ray  spreading  radius  in  plane  normal  to 
plane  of  incidence  at  reflection  point 

Distance  from  source  to  reflection  point  in 
x-y  plane 

Distance  from  source  to  reflection  point 

Distance  between  near-field  observation 
point  and  the  reflection  point 

Sine  of  VR 

X  component  of  normal  at  cylinder  reflec¬ 
tion  point 

Y  component  of  normal  at  cylinder  reflec¬ 
tion  point 

Spreading  factor 

Sine  of  WR 

SZN  X,  Y,  and  Z  components  of  unit  vector  of 

ray  from  reflection  point  to  source  in  RCS 

Theta  component  of  incident  ray  direction 

Variable  used  in  various  near-field  calcu¬ 
lations  to  indicate  the  theta  angle  between 
two  points 

Variable  used  in  various  near-field  calcu¬ 
lations  to  indicate  the  theta  angle  between 
two  points 

Parameter  used  in  transition  function 

X  component  of  source  vector  tangent  to 

tangent  point  1  (2-D) 

X  component  of  source  vector  tangent  to 

tangent  point  2  (2-D) 
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Y  component  of  source  vector  tangent  to 

tangent  point  1  (2-D) 

Y  component  of  source  vector  tangent  to 

tangent  point  2  (2-D) 

X,Y  components  of  unit  vector  tangent  to 
cylinder  reflection  point  in  RCS  (2-D) 

X,Y,Z  components  of  incident  field  polar¬ 
ization  unit  vector  parallel  to  plane  of 
incidence 


X,Y,Z  components  of  incident/reflected 
field  polarization  unit  vector  perpen¬ 
dicular  to  plane  of  incidence 

X,Y  components  of  unit  normal  to  cylinder 
reflection  point  in  RCS  (2-D) 

X,Y,Z  components  of  reflected  field  polari¬ 
zation  unit  vector  parallel  to  plane  of 
incidence 


Elliptical  angle  defining  reflection  point 
in  RCS  x-y  plane 

X,Y,Z  components  of  unit  vectors  defining 
source  coordinate  system  axes  in  RCS 

Phi  angle  defining  propagation  direction  in 
cylinder  reflection  point  coordinate  system 

Point  which  lies  on  the  infinite  cylinder 
whose  z  component  is  equal  to  the  inter¬ 
section  point  on  the  z  axis  with  the  more 
positive  end  cap  and  whose  x  and  y  values 
are  based  on  the  elliptical  cylinder's 
radii  at  the  phi  angle  of  a  vector  from  the 
RCS  origin  to  the  field  point 

Point  which  lies  on  the  infinite  cylinder 
whose  z  component  is  equal  to  the  inter¬ 
section  point  on  the  z  axis  with  the  more 
negative  end  cap  and  whose  x  and  y  values 
are  based  on  the  elliptical  cylinder's 
radii  at  the  phi  angle  of  a  vector  from  the 
RCS  origin  to  the  field  point 


I/O  VARIABLES: 
A.  INPUT 


FLDPT 


LNRFLD 


Location  of  reflection  point  in  x,  y,  z 

reference  coordinate  system  (RCS) 

coordinates 

Source  location  in  x,y,z  components 

The  x,y,z  components  of  tangent  point  1 

(2-D) 

The  x,y,z  components  of  tangent  point  2 

(2-0) 

LOCATION 

/GEOMEL/ 

/GEOMEL/ 

/BNDSCL/ 

/DIR/ 

/GEOMEL/ 

/DIR/ 

/DIR/ 

/THPHUV/ 

/THPHUV/ 

/BNDSCL/ 

/NEAR/ 

/NEAR/ 

/CLRFC/ 

/DIR/ 

/PIS/ 

/DIR/ 

/DIR/ 


VTS 

/BNDSCL/ 

VXS 

/SORINF/ 

XS 

/SORINF/ 

zc 

/GFOMEL/ 

B.  OUTPUT 

LOCATION 

EHPH 

/FUDG/ 

EHTH 

/FUDG/ 

EPH 

F.P. 

ESPH 

/FUDG / 

ESTH 

/FUDG/ 

ETH 

F.P. 

LTRF 

/FUDG/ 

RG 

/FUDG/ 

RHOl 

/FUDG/ 

SMAG 

/FUDG/ 

TRAN 

/FUDG/ 

XR 

/FUDG/ 

CALLING  ROUTINE: 

SCTCYL 

CALLEO  ROUTINES: 

I** 


,  v,'«v*  • 


PLAINT 


TPNFLD 


RFDFIN  WLKBCK 

REFERENCE: 

A.  R.  J.  Marhefka,  "Analysis  of  Aircraft  Wing-Mounted  Antenna 
Patterns,"  Report  2902-25,  June  1976,  The  Ohio  State  University 
ElectroScience  Laboratory,  Department  of  Electrical  Engi¬ 
neering;  prepared  under  Grant  No.  NGL  36-008-138  for  National 
Aeronautics  and  Space  Administration. 


PURPOSE:  To  execute  a  coordinate  reflection  for  the  geometry  pro¬ 
cessor. 

METHOD:  The  axis  along  which  the  reflection  Is  to  take  place  is 

determined,  that  coordinate  sign  Is  changed,  and  control  is  returned 
to  the  calling  subroutine. 

INTERNAL  VARIABLES: 


VARIABLE 


DEFINITION 


JAXIS 


I/O  VARIABLES: 
A.  INPUT 


Input  argument  determining  axis  of  reflec¬ 
tion:  IRF=1  for  X  axis,  -2  for  Y  axis,  =3 
for  Z  axis 

Variable  IRF  minus  2  for  arithmetic  IF 
statement 


LOCATION 


*  •  **-  ***  •  «  «  _  W  _  j  m  «  «,  V«  "  a  ***  •  a  •  »  *  •  *  k  *  *  *  .  *  *  •  ,  •  .  *  .  .*  •  •  ,  •  #  •  .  •  *  '  .  " 


B.  OUTPUT 


LOCATION 


CALLING  ROUTINES: 
HYRDRV 


SPWDRV 

CALLED  ROUTINES: 
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RETURN 


PURPOSE:  To  calculate  the  unobstructed  electric  field  due  to  single 
reflection  from  a  given  plate  from  a  unit  source. 


METHOO:  REFPLA  is  the  driver  routine  which  directs  all  the  ray 
tracing,  physics  and  field  calculations  for  determining  the  electric 
field  resulting  from  single  reflection  off  a  given  plate  in  a  given 
far-field  direction  or  to  a  given  near-field  point  from  a  unit 
source.  The  significant  geometry  is  shown  in  figure  1. 


SOURCE  . 

LOCATION  XS 


INCIDENT 

RAY 


D  -  OBSERVATION 
DIRECTION  FOR 
FAR  FIELD  OR  TO 
NEAR-FIELD  OBSERVATION 
POINT 

-  REFLECTED  RAY 


PLATE  MP 


SOURCE  IMAGE 
LOCATION  XIS 


Figure  1.  Geometry  for  Source  Ray  Reflection  from  Plate 


First  the  ray  path  from  the  source  image  location  in  the  desired 
direction  Is  checked  to  make  sure  it  does  pass  through  the  plate 
being  considered.  If  It  does  not,  the  theta  and  phi  components  of 
the  electric  field  are  set  to  zero,  and  no  other  calculations  except 
for  debug  functions  (If  requested)  are  performed  in  this  routine. 
If  reflection  can  occur,  the  ray  path  from  the  reflection  point  on 
the  plate  in  the  given  far-field  observation  direction  or  to  the 
given  near-field  observation  field  point  is  checked  for  obstruc¬ 
tions.  If  the  path  is  obstructed,  the  fields  are  set  to  zero,  and 
no  other  calculations  except  for  debug  functions  (if  requested)  are 
performed  In  this  routine.  If  the  path  Is  clear,  then  the  ray  path 
from  the  source  to  the  reflection  point  Is  checked.  If  this  path  is 
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blocked  by  another  plate  or  a  cylinder,  the  fields  are  set  to  zero, 
and  no  other  calculations  except  for  debug  functions  (if  requested) 
are  performed.  If  this  path  is  clear,  reflection  from  the  given 
plate  did  occur  and  the  complete  ray  path  is  unobstructed. 

The  source  field  pattern  factor  from  the  source  image  location  is 
computed  by  calling  subroutine  SOURCE  and  multiplying  the  returned 
values  by  the  reflection  coefficient.  Next  the  phase  factor  is 
computed.  For  far  field,  this  will  refer  the  field  back  to  the 
origin  of  the  reference  coordinate  system.  For  near  field,  the 
phase  factor  includes  the  spherical  wave  spread  factor.  Now  the 
theta  and  phi  components  of  the  electric  field  are  computed.  The 
electric  field  in  theta  and  phi  components  is  given  by: 


E  =  (EF  9 


EG  *) 


theta  component  phi  component 
of  source  of  source 

factor  factor 


j.?ir(Xrs  •  D) 


EX  -  phase  factor 


,  for  far  field 


E  =  (EF  0 


EG  ♦) 


-  j  2ttSNF 


•,  for  near  field. 


theta  component  phi  component 

of  source  of  source 

factor  factor 


EX  -  phase  factor 
where  SNF  =  I FLDPT  -  XlS I 


The  x,y,z  components  of  the  electric  field  are  computed  and  added  to 
the  previous  components  due  to  other  reflection-diffraction  inter¬ 
actions.  The  values  are  stored  in  the  /FLDXYZ/  common  block. 


If  the  debug  capabilities  have  been  requested,  the  field  magnitude 
is  computed.  The  magnitude,  theta  and  phi  components  are  printed  on 
file  LUPRNT. 

INTERNAL  VARIABLES: 


VARIABLE 


DEFINITION 


CPHI 


Cosine  of  PH IR 
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CTHI 

Cosine  of  THIR 

D 

Unit  vector  x,y,z  components  of  ray  propa¬ 
gation  direction  after  reflection  in  RCS 

DHIT 

Distance  from  source  to  reflection  point 
(from  PLAINT) 

DHT 

Distance  from  source  to  hit  point  (from 
PLAINT  and  CYLINT) 

DICOEF 

Unit  vector  <,y,z  components  of  incident 
ray  propagation  direction  in  RCS 

EF 

Pattern  factor  for  theta  component  of 
source  field  in  RCS 

EG 

Pattern  factor  for  phi  component  of  source 
field  in  RCS 

EIX 

X  component  of  source  factor 

EIY 

Y  component  of  source  factor 

EIZ 

Z  component  of  source  factor 

ERP 

Phi  component  of  reflected  field  in  RCS 

ERT 

Theta  component  of  reflected  field  in  RCS 

EX 

Complex  phase  factor 

FLDMAG 

The  electric  field  magnitude 

FLDPT 

The  x,y,z  components  of  the  field  point 
location 

FX.FY.FZ 

The  x,y,z  components  of  the  accumulated 
electric  field  from  all  geometry 
interactions 

GAM 

Phase  distance  to  origin  (dot  product  of 
image  location  and  reflected  ray  propaga¬ 
tion  direction) 

LOEBUG 

Logical  variable  set  true  if  debug 
requested 
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LHIT 

LNRFLD 

LUPRNT 

MP 

N 

NI 

NJ 

PHIR 

PHSR 

SNF 

SPHI 

STHI 

THIR 

THSR 

TP  I 
VAX 

VXI 

XI 


Set  true  if  ray  intersects  a  plate  or 
cylinder  (from  PLAINT  or  CYLINT) 

Flag  to  indicate  if  near-field  ( LNRFLD= 1 ) 
or  far-field  (LNRFLD=0)  calculations  were 
requested 

Output  file  number 

Plate  from  which  reflection  occurs 

DC  loop  variable 

OO  loop  variable 

DO  loop  variable 

Phi  component  of  incident  ray  propagation 
direction  in  RCS 

Phi  component  of  ray  propagation  direction 
after  reflection  in  RCS 

Distance  between  source  image  and  field 
point 

Sine  of  PHIR 
Sine  of  THIR 

Theta  component  of  incident  ray  propagation 
direction  in  RCS 

Theta  component  of  ray  propagation  direc¬ 
tion  after  reflection  in  RCS 

2tt 

X,Y,Z  components  defining  unit  vectors  of 
the  source  image  coordinate  system  axes  in 
RCS 

Array  of  components  defining  unit  vectors 
of  the  source  image  coordinate  system  axes 
in  RCS 

Triply  dimensioned  array  of  image  locations 
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XQS 

KS 

I/O  VARIABLES: 

A.  INPUT 
D 

FLDPT 

FX 

FY 

FZ 

LDEBUG 

LNRFLD 

LUPRNT 

MP 

PHSR 

THSR 

TPI 

VXI 

XI 

XS 

B.  OUTPUT 
ERP 
ERT 

FX 


X,Y,Z  components  of  source  image  location 
(single  reflection  from  plate  MP) 

X,Y,Z  components  of  source  image  location 

Source  location  in  x,y,z  RCS 

LOCATION 

/DIR/ 

/NEAR/ 

/FLDXY Z/ 

/FLDXYZ/ 

/FLDXYZ/ 

/TEST / 

/NEAR/ 

/ADEBUG/ 

F.P. 

/DIR/ 

/DIR/ 

/PIS/ 

/IMAINF/ 

/IMAINF/ 

/SORINF/ 

LOCATION 

F.P. 

F.P. 

/FLDXYZ/ 


REFPLA 


FY 

/FLDXYZ/ 

FZ 

/FLDXYZ/ 

CALLING  ROUTINE: 
GTDORV 

CALLED  ROUTINES: 

ASSIGN 

BEXP 

CYLINT 

NFD 


PLAINT 


REFPLA 


( GTD) 
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SPECIFY  SOURCE 
IMAGE  LOCATIONS 


IF  NEAR  FIELD 
REQUESTED,  COMPUTE 
OBSERVATION  DIREC¬ 
TION  D  AND  ITS  ANGLES 
THSR,  PHSR 
NFD 


IF  NEAR  FIELD 
REQUESTED,  COMPUTE 
DISTANCE  BETWEEN 
SOURCE  IMAGE  AND 
FIELD  POINT 
SMAGNF 


/DOES\ 

./  RAY  FROM^ — 
^  SOURCE  IMAGE  ^ 
PASS  THRU  PLATE  MP 
\(OID  REFLECTION/ 
P^OCCUR^/ 


/DOESP 
""REFLECTED^ 
RAY  HIT 
ANOTHER  PLATJ 


/DOESP 
REFLECTED 
RAY  HIT  A 
^CYLINDER  . 


KNOWING  RADIATION 
DIRECTION,  COMPUTE 
INCIDENT  RAY 
DIRECTION 

THETA  AND  PHI  ANGLES 
PEFBP 


COMPUTE  INCIDENT  RAY 
x,  y,  z  COMPONENTS 


/  DOES  P 
RAY  FROM 
SOURCE  HIT  A 
PLATE  BEFORE 
^REFLECTION/ 


./  DOES  P^ 
/  RAY  FROM  P 
SOURCE  HIT 
A  CYLINDER  BEFORE 

Preelection?^/ 


SPECIFY  SOURCE  IMAGE 
AXES 


COMPUTE  SOURCE  FIELD 
FACTOR 

SOURCE 


MULTIPLY  SOURCE  FIELD 
FACTOR  BY  REFLECTION 
COEFFICIENT 


I  COMPUTE  PHASE  FACTOR  I 


COMPUTE  THETA  AND 
PHI  FIELD  COMPONENTS 


PURPOSE:  Read  common  blocks  and  data  sets  from  the  checkpoint  file. 


METHOD:  RESTRT  searches  for  the  desired  checkpoint  on  the  specified 
logical  unit.  The  commons  are  read,  the  peripheral  files  are 
opened,  and  the  data  sets  are  restored.  The  search  continues  until 
the  desired  checkpoint  has  been  read.  If  the  desired  checkpoint  is 
not  found,  the  routine  writes  an  error  message  and  terminates 
execution. 

INTERNAL  VARIABLES: 


VARIABLE 


DEFINITION 


ICHKPT 


The  checkpoint  number  that  will  be  used  for 
restarting  (default  =  1) 


ICKLOP 


Loop  index  for  reading  commons  and  periph¬ 
eral  files 


Flag  indicating  that  an  end-of-file  occur¬ 
red  while  reading  the  checkpoint  file 


ISDBON 

ISTRDT 


Flag  for  debug  command 

Flag  to  tell  when  to  store  data  on  periph¬ 
eral  files 


LOCTPO 

LUFILE 


Pointer  for  argument  list 

Checkpoint  logical  unit  number  specified  by 
RSTART  command 


NAMCPF 


Temporary  Hollerith  location 

Name  of  module  asked  for  by  RESTRT  command 
(default  =  INPUT) 


File  name  or  data  set  name 


NDXNCD 


Index  of  a  keyword  in  the  NCODES  array 


NREAD 


Flag  to  tell  RWCOMS  to  read  common  areas 
from  IOCKPT 


NUMCPF 


Pointer  to  keyword  of  module  for  which 
RESTRT  is  requested 


I/O  VARIABLES: 
INPUT 
IOCKPT 
ISOFF 
ISON 
KOLCOL 
KOLNAM 
KWCHKP 
KWNAME 
LUPRNT 
NAMSEG 
NARGTB 
NCOOES 
NDATBL 
NDEBUG 
NDXBLK 
NOPCOD 
NPDATA 
NPTASK 
NTSKTB 
NUMCHK 
RSTART 


LOCATION 

/SYSFIL/ 

/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/SEGMNT/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/SCNPAR/ 

/SEGMNT/ 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/SYSFIL/ 

/SEGMNT/ 


SEGTBL 


RESTRT 


(INPUT) 


« 


B.  OUTPUT 
CHKWRT 
DBGPRT 
IERRF 
IMDCHK 
INTARG 
IRSTRT 
LSTSYS 
LUOBUG 
NDEBUG 
NFINCD 
NPTASK 
NTSKTB 
RSTART 
RSTRTA 

CALLING  ROUTINE: 
INPDRV 

CALLED  ROUTINES: 

ASSIGN 

CONVRT 

ERROR 

GETSYM 

POST IP 


PUTSYH 


LOCATION 


/SYSFIL/ 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 

/ARGCOM/ 

/ADEBUG/ 

/SYSFIL/ 

/ADEBUG/ 

/SCNPAR/ 

/IOFLES/ 

/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/SYSFIL/ 


ROEFIL 

RWCOMS 

RWFILS 

STATIN 

STATOT 

WLKBCK 


1.  NAME:  RFDFIN  (GTO) 

2.  PURPOSE:  To  determine  the  reflection  point  on  an  elliptic  cylinder 
for  a  given  source  and  observation  location  in  the  near  field  of  the 
cylinder. 

3.  METHOD:  This  subroutine  solves  a  polynomial  equation  which  is  based 
on  the  geometry  and  satisfies  Snell's  law.  The  roots  define 
possible  reflection  point  locations.  The  true  point  is  singled  out 
using  the  laws  of  reflection.  This  procedure  is  explained  in 
chapter  IV  of  reference  A.  The  coefficients  are  given  on  page  100 
of  the  reference. 

INTERNAL  VARIABLES: 


VARIABLE  DEFINITION 


A 

B 

CA 

RT 

S 

SM 


UR 


Cylinder  radius  along  the  x  axis 

Cylinder  radius  along  the  y  axis 

Complex  coefficients  of  sixth  order  poly¬ 
nomial  equation 

Roots  of  polynomial  equation 

Smallest  distance  from  source  to  reflection 
point  to  observation  point 

Distance  from  source  to  reflection  point 
plus  the  distance  from  the  reflection  point 
to  the  observation  point 

Z  component  of  reflection  point  on  cylinder 
in  RCS 


VI 


X,Y,Z  components  of  reflected  ray  propaga¬ 
tion  direction 


VIM 


Normalization  constant  for  VI 


VM  Elliptical  angle  defining  possible  reflec¬ 

tion  point  on  cylinder 

VR  Phi  angle  defining  x  and  y  components  of 

reflection  point 


XC 


X,Y,Z  components  of  the  observation  point 
in  the  near  field  of  the  cylinder 


RFDFIN 


(STD) 


XR 

XS 

I/O  VARIABLES: 

A.  INPUT 
A 
B 

XC 
XS  * 

B.  OUTPUT 
UR 
VI 
VR 

CALLING  ROUTINES: 

GEOMPC 
RCLRPL 
REFCYL 
RPLRCL 

CALLED  ROUTINES: 

BTAN2 
POLYRT 
REFERENCE: 

A.  R.  J.  Marhefka,  "Analysis  of  Aircraft  Wing-Mounted  Antenna 
Patterns,"  Report  2902-25,  June  1976,  the  Ohio  State  University 
ElectroSclence  Laboratory,  Department  of  Electrical  Engineer¬ 
ing;  prepared  under  Grant  No.  NGL  36-008-138  for  National 
Aeronautics  and  Space  Administration. 


X,Y,Z  components  of  reflection  point  loca¬ 
tion  on  cylinder 

Source  location 

LOCATION 

/GEOMEL/ 

/GEOMEL/ 

F.P. 

/SORINF/ 

LOCATION 

F.P. 

F.P. 

F.P. 


DETERMINE 
COEFFICIENTS  OF 
POLYNOMIAL 
EQUATION 


SOLVE  POLYNOMIAL 
EQUATION  TO  OBTAIN 
POSSIBLE  REFLEC¬ 
TION  POINTS 

POLYRT 


DETERMINE  PHYSICALLY 
CORRECT  REFLECTION 
POINTS  (SATISFY 
LAW  OF  REFLECTION) 


1.  NAME:  RFDFPT  (GTD) 

2.  PURPOSE:  To  compute  the  ray  path  for  a  source  ray  which  is 
reflected  by  the  cylinder  and  then  diffracted  by  a  given  edge  on  a 
given  plate. 

3.  METHOD:  The  reflection-  point  on  an  elliptic  cylinder  and  the 
diffraction  point  on  a  plate  edge  for  the  ref lected-diffracted  ray 
in  a  given  observation  direction  is  calculated  via  an  iterative 
process.  The  equations  are  based  on  a  first  order  Taylor  series 
approximation  to  the  equations  governing  the  laws  of  reflection  and 
diffraction.  The  details  of  the  analysis  are  given  on  pages  141-148 
of  reference  A.  The  iteration  process  follows  the  same  basic  scheme 
outlined  in  the  description  for  subroutines  RFPTCL  and  DFRFPT.  The 
initial  start-up  procedure  for  this  subroutine  is  composed  of 
locating  the  reflection  point  on  the  cylinder  for  a  known  diffrac¬ 
tion  point  which  is  taken  to  be  on  the  corners  of  the  plate  edge 
under  consideration.  The  details  of  this  procedure  are  discussed  on 
pages  149-154  of  reference  A.  Pertinent  geometry  is  shown  in 
figure  1.  To  avoid  the  2-n-to-O  transition  in  4>  (a  numerical  jump  in 
the  variable  representing  the  angle),  the  reference  $  value  is  ro¬ 
tated  to  place  this  branch  cut  behind  the  cylinder  (shadowed  from 
the  plate  edge). 


Figure  1.  Illustration  of  Ray  Reflected  from  Cylinder  and 
then  Diffracted  by  a  Plate  Edge 


RFDFPT  (GTD) 


The  iteration  begins  with  an  initial  ref lected-diffracted  ray  which 
satisfies  the  laws  of  diffraction  and  reflection.  Starting  data  are 
obtained  in  one  of  two  ways.  If  a  previous  call  to  this  routine 
(for  the  same  plate  edge  and  source)  has  successfully  found  the 
ref lected-diffracted  ray  path,  this  previous  path  is  used  as 
starting  data.  Otherwise  the  starting  diffraction  point  is  defined 
on  the  corner  of  the  edge  closest  to  the  cylinder.  Then  the 
corresponding  reflection  point  is  found  by  enforcing  Snell's  law. 
The  first  method  is  preferred  to  the  second  since,  in  general,  far- 
field  ray  directions  (6)  in  subsequent  calls  to  RFDFPT  will  not 
differ  greatly.  For  example,  in  calculating  a  far-field  pattern 
cut,  the  far-field  e  and  $  angles  will  differ  by  only  a  few  degrees, 
and  the  closeness  of  the  starting  point  will  lead  to  fewer  itera¬ 
tions  in  order  to  obtain  convergence. 

The  path  of  the  starting  ray  defines  the  initi_al  cylinder  reflection 
point  (XFO  and  the  edge  diffraction  point  (XD).  In  almost  every 
instance,  the  resulting  radiation  direction  of  this  ray  will  not  be 
the  desired  radiation  direction.  The  angular  difference  between 
(0,<t>)  of  the  starting  ray  (THOR,  PHORP)  and  the  (0,$)  of  the  desired 
direction  (THSR,  PHSRP)  is  divided  into  a  number  of  small  angular 
steps  (A6,A<t>)  =  (DTSR,  DPSR).  The  purpose  of  the  iteration  is  to 
move  the  reflection  and  diffraction  points  from  their  initial  posi¬ 
tions  in  small  steps  corresponding  to  angular  changes  ( A 0 ,_4<j> )  50 
that  when  the  iteration  is  complete  the  resulting  XR  and  XD  will 
define  the  reflection  and  diffraction  points  that  give  the  desired 
D-directed  ray.  The  number  of  steps  to  be  taken  (IVD)  is  deter¬ 
mined  from  the  starting  data.  Should  convergence  not  be  reached  in 
IVD  steps,  the  number  of  steps  is  doubled  (up  to  32  steps)  and  the 
iteration  repeated.  The  doubling  process  is  the  outer  loop  of  the 
flowchart.  Should  convergence  be  reached  with  IVD  steps  and  the 
Snell's  law  error  be  significantly  smaller  than  required,  IVD  for 
the  plate  and  edge  under  consideration  is  halved  prior  to  exiting 
the  routine. 

The  iterations  which  step  through  the  ( 0 , 4> )  angles  by  (A0,A4>) 
correspond  to  the  inner  loop  of  the  flowchart.  Each  iteration  has 
three  steps: 

(1)  Compute  the  d if fraction  point  (XD)  from  known  reflection  point 

(XR),  source  point  (XS)  and  edge  unit  vector  (9).  This  is  done 
by  a  simple  application  of  Snell's  law.  All  the  far-field 

calculations  are  contained  in  this  subroutine,  but,  to  deter¬ 
mine  the  near-field  diffraction  point,  subroutine  DFPTWD  is 
called. 

(2)  The  change  in  cylinder  elliptic  angle  (DV)  and  z  coordinate 
(DU)  are  computed  from  a  Taylor  series  expansion.  The  expan¬ 
sion  requires  the  calculation  of  functions  and  partial 
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RFDFPT 


(STD) 


derivatives  of  equations  defining  elliptic  angle  (VR)  and 
2  coordinate  (UR)  in  terms  of  the  angles  (0,4)).  The  equations 
are  given  in  reference  A. 

(3)  The  coordinates  of  XR  are  computed  from  the  new  values  of  UR 
and  VR. 

At  the  end  of  the  prescribed  number  of  iterations,  the  initial 
observation  direction  has  been  stepped  slowly  to  the  desired  direc¬ 
tion  and  the  initial  reflection-diffraction  points  have  been  stepped 
from  their  initial  values  to  candidate  reflection-diffraction 
points.  Snell's  law  is  then  applied  to  the  final  reflection  and 
final  diffraction  points  to  see  if  they  qualify  as  the  bona  fide  ray 
path.  If  the  error  is  sufficiently  small,  the  outer  loop  is  exited. 
Otherwise,  the  number  of  steps  is  doubled,  as  described  above. 
Should  the  routine  not  converge  with  32  steps  (the  maximum  number), 
a  warning  message  is  printed  on  LUPRNT. 

This  routine  is  called  by  RCLDPL  (reflection  from  a  cylinder, 
diffraction  from  a  plate)  and  also  by  DPLRCL  (diffraction  from  a 
plate,  reflection  from  a  cylinder).  DPLRCL  only  calls  this  routine 
for  the  near-field  case  since,  if  the  observation  point  and  the 
source  point  are  reversed,  the  ray  path  would  be  the  same. 

INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

A  The  x  axis  radius  of  the  cylinder 

B  The  y  axis  radius  of  the  cylinder 

BCD  Diffraction  limit  for  ray  reflected  by  the 

cylinder  and  diffracted  from  the  plate 

D  The  unit  vector  of  the  observation 

direction 

DC  X,Y,Z  components  of  diffracted  ray  propaga¬ 

tion  used  in  iteration 

DCP  X,Y  components  of  phi  polarization  unit 

vector  for  diffracted  ray  used  in  iteration 

DCT  X,Y,Z  components  of  theta  polarization  unit 

vector  for  diffracted  ray  used  in  iteration 

DE  Dot  product  of  diffracted  ray  direction  and 

edge  vector  of  edge  ME 


RFOFPT 


(GTD) 


•-.'V 

Sis 

>:• 


DOTP 

DPSR 

DR 

DRM 

DRP 

DRT 

DRU 

DRV 

DTSR 

Co 

DV 

ERC 

FI 

FLDPT 

FP 

FT 

FU 

FV 

GI 


Test  parameter  used  to  determine  if  reflec¬ 
tion  is  legal 

Phi  angle  increment  size 

X,Y.,I  components  of  ray  direction  between 
refle'.tion  and  diffraction  points 

Distance  from  reflection  point  to  the 
diffraction  point 

Partial  derivative  of  DR  with  respect  to 

phi 

Partial  derivative  of  DR  with  respect  to 

theta 

Partial  derivative  of  DR  with  respect  to  UR 

Partial  derivative  of  DR  with  respect  to  VR 

Theta  angle  increment  size 

Change  in  UR  for  one  iteration  using  Taylor 
series  expansion 

Change  in  VR  for  one  iteration  using  Taylor 
series  expansion 

Error  detection  variable 

Equation  governing  the  law  of  reflection 

The  x,y,z  components  of  the  near-field 
observation  point 

Partial  derivative  of  FI  with  respect  to 
phi 

Partial  derivative  of  FI  with  respect  to 
theta 

Partial  derivative  of  FI  with  respect  to  UR 
Partial  derivative  of  FI  with  respect  to  VR 
Equation  governing  the  law  of  reflection 
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GP 

GT 

GU 

GV 

IVD 

LNRFLD 

LRDC 

ME 

MEP 

MP 

PHCR 

PHOR 

PHORP.PHSPR 

PHSR 

PHWR 

PI 

SNM 

SNPX 


RFOFPT  (GTD) 


Partial  derivative  of  GI  with  respect  to 

phi 

Partial  derivative  of  GI  with  respect  to 

theta 

Partial  derivative  of  GI  with  respect  to  UR 

Partial  derivative  of  GI  with  respect  to  VR 

Stored  number  of  steps  used  in  iteration 

Flag  to  indicate  near-field  (LNRFLD  =  1)  or 
far-field  (LNRFLD  =  0)  calculations  were 
requested 

Set  true  if  starting  point  data  are  avail¬ 
able  from  previous  pattern  angle 

Edge  on  plate  MP  where  diffraction  occurs 

Array  which  contains  the  number  of  edges  on 
each  plate 

Plate  where  diffraction  occurs 

Phi  component  of  diffracted  ray  direction 
used  in  iteration 

Phi  component  of  diffracted  ray  direction 
from  previous  time  RFDFPT  was  called  (or 
present  value  for  next  time  routine  is 
called) 

Phi  angle  of  diffracted  ray  direction  in 
rotated  RCS  system  (branch  cut  placed 
behind  cylinder) 

The  phi  angle  of  the  observation  direction 

Branch  cut  displacement  angle  for  the 
diffraction  point  along  edge  ME  of  plate  MP 

IT 

Normalization  constant  for  cylinder  tangent 

Partial  derivative  of  SNX  with  respect  to 
angle  VR 
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RFDFPT  (GTD) 


SNPY 

SNX 

SNY 

STP 

THCR 


Partial  derivative  of  SNY  with  respect  to 
angle  VR 

X  component  of  normal  to  cylinder 

Y  component  of  normal  to  cylinder 

Number  of  steps  used  in  iteration 

Theta  component  of  diffracted  ray  direction 
used  in  iteration 


THOR  Theta  component  of  diffracted  ray  direction 

from  previous  time  RFDFPT  was  called  (or 
for  next  time  routine  is  called) 

THSR  The  theta  angle  of  the  observation 

direction 


TPI  2  it 

UCD  Z  component  of  reflection  point  location  on 

cylinder  for  cyl inder-ref lected  ray 
diffracted  by  a  corner  on  edge  ME  of  plate 
MP  used  as  the  starting  point  if  previous 
data  do  not  exist 

UR  Z  component  of  reflection  point  location  on 

cylinder 

URO  Stored  components  defining  z  component  of 

starting  reflection  point  locations  on 
cyl inder 

V  Matrix  of  edge  unit  vectors  for  all  edges 

of  all  plates 

VCD  Elliptical  angle  defining  reflection  point 

on  cylinder  (2-D)  for  ray  which  is 
reflected  by  cylinder  and  diffracted  by  a 
corner  on  edge  ME  of  plate  MP  used  for  the 
starting  point  location  if  previous  data  do 
not  exist 

VI  X,Y,Z  components  of  unit  vector  of  ray 

incident  on  cylinder 

VIM  Distance  from  source  to  reflection  point 
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RFDFPT  (GTD) 

VIU 

Partial  derivative  of  VI  with  respect  to  UR 

VIV 

Partial  derivative  of  VI  with  respect  to 
angle  VR 

VR 

Elliptical  angle  defining  reflection  point 
on  cylinder  (2-D) 

VRO 

Stored  elliptical  angles  defining  starting 
reflection  point  locations  on  cylinder 

X 

Matrix  which  contains  all  the  corner  loca¬ 
tions  for  all  the  plates 

XD 

X,Y,Z  components  of  diffraction  point 
location 

XR 

X,Y,Z  components  of  reflection  point 
location  on  cylinder 

XS 

The  x,y,2  components  of  the  source  location 
in  RCS 

I/O  VARIABLES: 

A. 

INPUT 

LOCATION 

A 

/GEOMEL/ 

B 

/GEOMEL/ 

BCD 

/BNDRCL/ 

D 

/DIR/ 

DE 

F.P. 

FLDPT 

/NEAR/ 

LNRFLD 

/NEAR/ 

LRDC 

F.P. 

LUPRNT 

/ADEBUG/ 

ME 

F.P. 

MEP 

/GEOPLA/ 

rTv 


PHWR 


PI 

THSR 
TP  I 
UCD 

V 

VCD 

X 

XS 

e.  OUTPUT 
DOTP 
DR 
DRM 
LRDC 
SNH 

VI 
VIM 
VR 
XD 
XR 


/BRNPHW/ 

/PIS/ 

/DIR/ 

/PIS/ 

/8NDRCL/ 

/GEOPLA/ 

/BNDRCL/ 

/GEOPLA/ 

/SORINF/ 

LOCATION 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 

F.P. 


CALLING  ROUTINES 


CALLED  ROUTINES: 


BTAN2 

DFPTWD 

NFD 

REFERENCE: 


A.  R.  J.  Marhefka,  "Analysis  of  Aircraft  Wing-Mounted  Antenna 
Patterns,"  Report  2902-25,  June  1976,  The  Ohio  State  University 
ElectroScience  Laboratory,  Department  of  Electrical  Engi¬ 
neering;  prepared  under  Grant  No.  NGL  36-008-138  for  National 
Aeronautics  and  Space  Administration. 
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(COMPUTE  STARTING  POINT  (pp.  ISO  152 
REFERENCE  A; 


I  PERFORM  IVOR  STEP 
ITERATIONS  TO 
NUMERICALLY  COMPUTE 
REFLECTION  ANO 
DIFFRACTION  POINTS  i 

I 

STEP  THROUGH  ANGLE! 

defining  oiffracted 

RAY  PROPAGATION 
OiRECTlON 


_ * 

COMPUTE  REFLECTION 
POINT  •  y  i 
COMPONENTS  ON 
CYLINDER 


m 
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WRITE  WARNING 
MESSAGE.  SET 
LRDC=FALSE 


RETURN 


NAME:  RFPTCL 


(GTD) 


PURPOSE:  To  calculate  the  reflection  point  on  the  elliptic  cylinder 
for  a  source  ray  reflected  in  a  given  direction.  The  routine  also 
computes  cylinder  reflection  points  for  source  rays  that  are 
reflected  by  a  given  plate  and  then  reflected  by  the  cylinder. 

METHOD:  Figures  1,  2  and  3  show  the  geometry  involved.  The  reflec¬ 
tion  point  for  a  ray  reflected  in  a  direction  defined  by  the  pM 
angle  PHSR  is  calculated  via  an  iterative  process.  The  routine 
starts  with  the  tangent  ray  nearest  to  the  reflected  ray  direction 
(or  other  nearby  reflected  ray  whose  reflection  point  is  known)  and 
steps  along  the  cylinder  surface,  calculating  the  approximate 
reflection  point  for  each  reflected  ray  phi  angle  PHPR  (which  is 
stepped  from  PHOR  to  PHSR  in  evenly  spaced  steps).  Each  reflection 
point  calculation  uses  the  previous  reflection  point  as  a  reference. 
As  long  as  the  steps  are  sufficiently  small,  the  approximation  is 
accurate.  The  equations  are  based  on  a  first  order  Taylor  series 
approximation  of  the  equation  governing  the  laws  of  reflection. 
Further  details  are  given  on  pages  102-104  of  reference  A.  The 
point  obtained  at  the  end  of  the  process  is  the  estimated  reflection 
point.  The  routine  then  takes  the  sum  of  dot  products  of  the 
cylinder  normal  and  the  Incident  and  reflected  rays  (which  should  be 
zero  in  order  to  satisfy  the  law  of  reflection).  If  it  is  larger 
than  some  minimal  amount,  the  number  of  iteration  steps  for  angle 
PHPR  is  doubled  and  the  calculation  is  redone.  If  the  error  is  much 
smaller  than  necessary,  the  number  of  steps  used  in  the  next  calcu¬ 
lation  is  divided  by  two. 

Once  a  reflection  point  is  calculated  for  a  particular  geometry,  the 
elliptical  angle  defining  the  reflection  point  (VRO(MR))  is  saved, 
along  with  the  number  of  steps  used  to  calculate  it  (IVD(MR))  for 
the  next  time  RFPTCL  is  called  for  the  same  geometry.  Since  the 
next  pattern  angle  Is  likely  to  be  quite  close  to  the  previous  one, 
this  gives  the  computer  a  good  starting  point  in  defining  the  next 
reflection  point,  hence  minimizing  computer  time.  LRFC  is  a  logical 
variable  which  if  true  tells  the  user  that  there  are  data  from  the 
previous  pattern  angle  available  to  compute  the  next  reflection 
point.  If  a  reflection  does  not  occur,  LRFC  is  set  false,  and  the 
next  time  the  routine  Is  called,  it  will  start  at  the  nearest 
tangent  point. 

An  Important  note  to  users  Is  that  if  the  same  problem  is  run  with  a 
different  starting  angle  or  with  a  different  angle  increment,  the 
answers  may  not  be  precisely  the  same  due  to  the  iterative  approach. 
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Figure  1.  Illustration  of  Cylinder  Reflection  Point 
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Figure  2.  Geometry  for  Calculating  Reflection  Point 


Figure  3.  Illustration  of  Iterative  Method  Used  in  Computing 
the  Cylinder  Reflection  Point 


4.  INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

A  Radius  of  cylinder  along  the  x  axis 

B  Radius  of  cylinder  along  the  y  axis 

BTI  This  defines  unit  vectors  for  the  two  rays 

reflected  by  each  plate  and  tangent  to  the 

cylinder.  The  unit  vector  for  the  source 
ray  reflected  from  plate  MP  tangent  to 
tangent  point  1  is  given  by: 

T1  =  x  *  BTI (MP , 1 )  +  y  *  BTI(MP,2) 
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v  v  v  V 


>  /  /vVv' 
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V 


RFPTCL 


(STD) 


DPX.DPY 


The  unit  vector  for  the  source  ray 
reflected  from  plate  MP  tangent  to  tangent 
point  2  is  given  by: 

A  a 

T2  =  x  *  BTI ( MP , 3 )  +  y  *  BTI(MP,4) 

This  defines  unit  vectors  of  the  two  source 
rays  tangent  to  the  cylinder.  The  unit 
vector  for  the  source  ray  tangent  to 
tangent  point  1  is  given  by: 

A  "  ~ 

T1  =  x  *  BTS(l)  +  y  *  BTS ( 2 ) 

The  unit  vector  for  the  source  ray  tangent 
to  tangent  point  2  is  given  by: 

T2  =  x  *  BTS ( 3 )  +  y  *  BTS(4) 

Cosine  of  PHPR 
Cosine  of  PHSR 
Cosine  of  VR 

Normalization  constant  for  reflection  point 
normal  vector 

One  half  the  difference  between  the  dot 
products  of  the  reflected  ray  direction  and 
cylinder  unit  normal  and  the  incident  ray 
direction  and  cylinder  unit  normal 

Size  of  angle  step  used  in  iteration 

X  and  Y  components  of  partial  derivative  of 
reflected  ray  direction  with  respect  to  phi 
observation  angle 

Dot  product  of  incident  ray  unit  vector  and 
cylinder  unit  normal 

Dot  product  of  reflected  ray  propagation 
direction  unit  vector  and  cylinder  unit 
normal 

Change  in  angle  VR 


RFPTCL  (GTD) 


DVB 

DVT 

DX,DY 

ERC 

ERCA 

FI 

IVD 

IVDM 

LRFC 


MP 


MPXR 

MR 


Partial  derivative  of  the  reflection  law 

equation  (FI)  with  respect  to  elliptical 
angle  VR 

Partial  derivative  of  the  reflection  law 

equation  (FI)  with  respect  to  the  phi  angle 
of  the  observation  direction 

X  and  Y  components  of  unit  vector  of 
reflected  ray  (direction  defined  by  angle 
PHPR)  in  RCS 

Error  parameter  (sum  of  DS  and  DR) 

Absolute  value  of  ERC 

Equation  satisfying  the  law  of  reflection 

Number  of  iterations  used  to  find  reflec¬ 
tion  point  the  last  time  RFPTCL  was  called 
for  plate  HP 

Number  of  steps  used  in  iteration 

(Entering  routine)  set  true  if  reflection 
occurred  last  time  REFCYL  was  called. 
(LRFC  set  true  when  leaving  routine  if 
reflection  occurred  this  time) 

Used  to  specify  whether  source  or  source 

image  is  used 

MP=0  designates  source 

MP>0  designates  source  image  for  reflection 
from  plate  MP 

Maximum  number  of  plates  present 

Index  variable  (MP+MPXR+I)  for  storing  data 
for  next  call  to  RFPTCL 


PHE  Phi  angle  between  reflected  ray  direction 

and  tangent  point  2 

PHEP  Phi  angle  between  reflected  ray  direction 

and  tangent  point  1 


PHIR 


Phi  component  of  source  location  in  RCS 


PHOR 

PHORB 

PHORP 

PHPR 

PHSPR 

PHSR 

PI 

S 

SIPX.SIPY 

SIX.SIY 

SHPX.SNPY 

SNV 

SNX.SNY 

SPP 

SPS 


RFPTCL  (GTD) 


Reflected  ray  phi  angle  (stored  as  starting 
point  parameter  for  next  time  routine  is 
called) 

Phi  angle  defining  ray  tangent  to  tangent 
point  1 

Phi  angle  of  cylinder  reflected  ray  direc¬ 
tion  In  rotated  RCS 

Reflected  ray  phi  angle  (iterated  from  PHOR 
to  PHSR) 

Phi  angle  defining  reflected  ray  direction 
In  rotated  RCS 

Phi  component  of  reflected  ray  propagation 
direction  In  RCS 

it 

Distance  from  source  to  reflection  point  in 
x-y  plane 

X  and  Y  components  of  partial  derivative  of 
Incident  ray  vector  with  respect  to  ellip¬ 
tical  angle  VR 

X  and  Y  components  of  Incident  ray  propaga¬ 
tion  vector  In  RCS  (not  always  normalized) 

X  and  Y  components  of  partial  derivative  of 
cylinder  normal  at  reflection  point  with 
respect  to  elliptical  angle  VR 

Sine  of  VR 

X  and  Y  components  of  ray  normal  to 
cylinder  reflection  point  In  RCS  (not 
always  normalized) 


STP 


Sine  of  PHPR 
Sine  of  PHSR 

Number  of  steps  used  In  Iteration 


VR 

Elliptical  angle  defining  reflection  point 
In  RCS  x-y  plane 

VRO 

Elliptical  angles  defining  tangent  points 
for  source  ray  (or  source  ray  reflected 
from  plate)  tangent  to  cylinder 

VTI 

Array  of  elliptical  angles  defining  for 
each  plate  the  two  tangent  points  on  the 
cylinder  for  rays  which  are  reflected  from 
that  plate 

VTS 

• 

Consists  of  two  elliptical  angles  defining 
the  two  tangent  points  on  the  cylinder  from 
the  source 

XI 

Array  which  contains  the  source  Image  loca¬ 
tions  In  wavelengths  for  all  plate  single 
and  double  reflections 

XIS 

Source  location 

xs 

X,Y  and  Z  components  of  source  location 

I/O  VARIABLES: 

A. 

INPUT 

LOCATION 

A 

/GEOMEL/ 

B 

/GEOKiL/ 

BTI 

/BNDICL/ 

BTS 

/BNDSCL/ 

LRFC 

F.P. 

LUPRNT 

/ADEBUG/ 

HP 

F.P. 

MPXR 

/GROUND/ 

•i- 


RFPTCL 


(GTD) 


VTI 

/BNDICL/ 

VTS 

/BNDSCL/ 

XI 

/IMAINF/ 

XS 

/SORINF/ 

OUTPUT 

LOCATION 

DO 

F.P. 

DOTP 

F.P. 

LRFC 

F.P. 

S 

F.P. 

VR 

F.P. 

CALLING  ROUTINES: 

RCLRPL 

REFCYL 

RPLRCL 

CALLED  ROUTINE: 

BTAN2 

REFERENCE: 

A.  R.  J.  Marhefka,  "Analysis  of  Aircraft  Wing-Mounted  Antenna 
Patterns,"  Report  2902-25,  June  1976,  The  Ohio  State  University 
ElectroSclence  Laboratory,  Department  of  Electrical  Engi¬ 
neering;  prepared  under  Grant  No.  NGL  36-008-138  for  National 
Aeronautics  and  Space  Administration. 
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SPECIFY  SOURCE  LOCATION 


E  F  LECT  lON'V. 
PRESENT  LAST  TlME^V. 
^RFPTCL  WAS  CALLED  FOR  THIS*' 
V.  GEOMETRY?  .X* 

^^ILRFC  •  TRUE?]^^ 


(COMPUTE  STARTING  POINT) 

_ Jl _ 

SPECIFY  SOURCE  VECTORS  TANGENT 
TO  CYLINDER  AND  ELLIPTICAL 
ANGLES  DEFINING  TANGENT  POINTS 


COMPUTE  ANGLES  ANO  SPECIFY  WHICH 
TANGENT  VECTOR  IS  CLOSEST  TO  THE 
REFLECTED  RAY  PROPAGATION  DIREC 
TION  PHOR  (MR)  DEFINES  NEAREST 
TANGENT  VECTOR.  VRO  (MR)  DEFINES 
CORRESPONDING  TANGENT  POINT  IN 
RCS  *-y  PLANE  IVO  (MR)  -  1 


SPECIFY  STARTING  POINT 


PERFORM  IVD(MR)  STEP  NUMERICAL 
CALCULATIONS  TO  APPROXIMATE 
ELLIPTICAL  ANGLE  VR  DEFINING  THE 
REFLECTION  POINT  IN  RCS  »  y  PLANE 
(SEE  FIGURE  3) 


STEP  THROUGH  ANGLES 
(VARIABLE  •  PHPR) 


APPROXIMATE  NEW  ANGLE  VR.  USING 
VALUE  OBTAINED  FOR  PREVIOUS  VALUE 
OF  PHPR  CALCULATED  AS  A  REFERENCE 
(THE  LAST  TIME  THROUGH  THIS  LOOP 
PHPR  -  PHSR  ANO  ELLIPTICAL  ANGLE 
VR  THEN  DEFINES  THE  APPROXIMATE 
REFLECTION  POINT  FOR  RAY  PR  OP  A 
CATION  IN  DIRECTION  6) 


(CHECK  ACCURACY  OF  REFLECTION  POINT  ANGLE  CALCULATED) 


COMPUTE  ERROR  FUNCTION  ERC,  THE  SUM  OF 
DOT  PROOUCTS  OF  REFLECTION  POINT  UNIT 
NORMAL  ANO  THE  INCIDENT  AND  REFLECTED 
RAY  UNIT  VECTORS 
IRC-  0»n 


NAME:  ROMBNT 


(GTO,  MOM) 


PURPOSE:  To  numerically  compute  the  Integral  of  the  function  exp 
(-jkr)/kr  for  a  wire  segment  or  I(l/kr)3  +  j(l/kr)2J  exp  (-jkr)  for 
a  patch  observation  point. 

METHOD:  To  evaluate  the  field  due  to  a  segment,  a  local  cylindrical 
coordinate  system  Is  defined  with  origin  at  the  center  of  the  seg¬ 
ment  and  z  axis  in  the  segment  direction.  This  geometry  is  illus¬ 
trated  In  the  discussion  of  subroutine  TNEFLD.  Subroutine  ROMBNT  is 
called  by  subroutines  TNEFLD  and  SOURCE  to  evaluate  the  Integral 


■/ 


kA/2 


►  “jkr 


G  *  /  s-j^-  d(kz) 

-kA/2 


where 


I 


r  -  p'2  +  <z  -  z')2 


1/2 


and  other  symbols  are  defined  In  the  discussion  of  subroutine 
TNEFLD.  To  evaluate  the  magnetic  field  at  a  patch  observation 
point,  this  subroutine  Is  called  by  the  subroutine  TNHFLD  to  eval¬ 
uate  the  Integral 


e"jkr  d(kz) 


The  numerical  integration  technique  of  Romberg  Integration  with 
variable  interval  width  Is  used.  The  Romberg  Integration  formula  Is 
obtained  from  the  trapezoidal  formula  by  an  Iterative  procedure  (see 
reference  A).  The  trapezoidal  rule  for  Integration  of  the  function 
f(x)  over  an  Interval  (a,  b)  using  2k  subintervals  Is 

T0k  «  (b  -  a)/N  £l/2  fQ  +  t1  +  ...  +  fN-1  +  1/2  fJ 
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where 


N  =  2k 
fi  -  f(xi) 

x^  =  a  +  1  (b  -  a)/N 


These  trapezoidal  rule  answers  are  then  used  in  the  Iterative 
formula 


T  n  =  x/<4  “ 

m ,  n 


4®*p  _ 

m-l,n+l  m-l,n 


The  results  Tffln  may  be  arranged  in  a  triangular  matrix  of  the  form 


T01  T10 
T02  TU  T20 


where  the  elements  in  the  first  column,  Tgk*  represent  the  trape¬ 
zoidal  rule  results  and  the  elements  In  the  diagonal,  T|<g,  are  the 
Romberg  Integration  results  for  2k  sub intervals. 

Convergence  to  Increasingly  accurate  answers  takes  place  down 
the  first  column  and  the  diagonal  as  well  as  toward  the  right  along 
the  rows.  The  row  convergence  generally  provides  a  more  realistic 
Indication  of  error  magnitude  than  two  successive  trapezoidal  rule 
or  Romberg  answers. 

This  convergence  along  the  rows  Is  used  to  determine  the 
Interval  width  In  the  variable  Interval  width  scheme.  The  complete 
Integration  Interval  Is  first  divided  Into  a  minimum  number  of  sub¬ 
intervals  (presently  set  to  one)  and  Tgg,  Tgi,  and  Tig  are  computed 
on  the  first  subinterval.  The  relative  difference  of  Tgi  and  Tig  Is 
then  computed  and  If  less  than  the  error  criterion,  Rx,  Tig  Is 
accepted  as  the  Integral  over  that  Interval  and  Integration  proceeds 
to  the  next  Interval.  If  the  difference  of  Tgi  and  Tig  Is  too 
great,  Tn2»  Tn,  and  T20  are  computed.  The  relative  difference  of 
Tn  and  120  Is  then  computed  and  If  less  than  Rx,  T2g  Is  accepted  as 
the  Integral  over  the  subinterval.  If  the  difference  of  Tn  and  T2g 
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Is  too  great,  the  subinterval  is  divided  in  half  and  the  process 
repeated  starting  with  Tqq  for  the  new  left-hand  subinterval.  The 
subinterval  Is  repeatedly  halved  until  convergence  to  less  than  Rx 
Is  found.  The  process  is  repeated  for  successive  subintervals  until 
the  right-hand  side  of  the  integration  interval  is  reached.  When 
convergence  has  been  obtained  with  a  given  subinterval  size,  the 
routine  attempts  doubling  the  subinterval  size  for  a  few  times  to 
maintain  the  largest  subinterval  size  that  will  give  the  required 
accuracy.  Thus,  the  routine  will  use  many  points  in  a  rapidly 
changing  region  of  a  function  and  few  points  where  the  function  is 
smoothly  varying. 

Since  the  function  to  be  Integrated  is  complex,  the  convergence 
of  both  real  and  imaginary  parts  is  tested  and  both  must  be  less 
than  Rx.  The  same  subinterval  sizes  are  used  for  real  and  imaginary 
parts. 

When  the  field  of  a  segment  is  being  computed  at  the  segment's 
own  center  the  length  r  becomes 


(z  -  z' ) 


where  b  Is  the  wire  radius.  For  small  values  of  b,  the  real  part  of 
the  Integrand  Is  sharply  peaked  and  hence  difficult  to  Integrate 
numerically.  Therefore,  the  Integral  Is  divided  into  the  components 


kA/2 


-kA/2 


d(kz) 


ZA/2 

i* 

2 


d(kz) 


G'  must  be  computed  numerically,  however  the  integrand  Is  no  longer 
peaked.  G",  which  contains  the  sharp  peak,  can  be  computed  as 
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G**  =  21og  (J  b2  +  A2  +  A)/b 


To  further  reduce  Integration  time  for  the  self  term,  the  integral 
of  G'  Is  computed  from  -kA/2  to  0  and  the  result  doubled  to  obtain 
G\ 


INTERNAL  VARIABLES: 

VARIABLE 

BK 

OZ 

DZOT 

ELI 

EL2 

EP 

FNM 

FNS 

G1I 

G1R 

G2I 

G2R 

G3I 

G3R 

G4I 

G4R 

G5I 


DEFINITION 
Wire  radius,  b 

Subinterval  size  on  which  TQQ,  TQ.,  ...are 
computed  uu  Ui 

0.5  OZ 

-kA/2 

kA/2 

Tolerance  for  ending  the  integration 

Interval  kA/NM*10 

Real  number  equivalent  of  NM 

Real  number  equivalent  of  NS 

Imaginary  part  of  f ^ 

Real  part  of  f. 

Imaginary  part  of  f2 

Real  part  of  f^ 

Imaginary  part  of  f^ 

Real  part  of  f3 

Imaginary  part  of  f^ 

Real  part  of  f^ 

Imaginary  part  of  fg 


IT  (GTO,  MOM) 


Indicates  self-term  in:egration  when  equal 
to  zero,  or  patch  observation  point  when 
equal  to  one 

65536  =  Z16 

Present  subinterval  size  Is  kA/NS 

Counter  to  control  increasing  of  sub¬ 
interval  size 

Larger  values  retard  increasing  of  sub¬ 
interval  size 

Maximum  allowed  subinterval  size  Is  kA/NX 
Rx 

Imaginary  part  of  G 
Real  part  of  G 
A 

Relative  difference  of  Tnl  and  Tln  for 

imaginary  part 

Relative  difference  of  TQ,  and  T.Q  for  real 
part  4 

Relative  difference  of  T,,  and  T-q  for 

imaginary  part 

Relative  difference  of  T,.  and  T2Q  for  real 
part 

Imaginary  part  of  Tqq 
Real  part  of  Tqq 
Imaginary  part  of  TQ1 
Real  part  of  TQ1 
Imaginary  part  of  Tq2 
Real  part  of  TQ2 
Imaginary  part  of  T^Q 
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T10R 

Real  part  of  T^q 

Till 

Imaginary  part  of  T^ 

T11R 

Real  part  of  T^ 

T20I 

Imaginary  part  of  T^q 

T20R 

Real  part  of  T£q 

Z 

Integration  variable  at  left-hand  side  of 
subi’iterval 

ZE 

kA/2 

ZEND 

kA/2  -  EP  (EP  =  tolerance  term) 

ZP 

Integration  variable 

5.  I/O  VARIABLES: 

A.  INPUT 

LOCATION 

BK 

F.P. 

ELI 

F.P. 

EL2 

F.P. 

IJ 

F.P. 

B.  OUTPUT 

LOCATION 

ICALL 

/ADEBUG/ 

NUMWRD 

/ADEBUG/ 

SGI 

F.P. 

SGR 

F.P. 

6.  CALLING  ROUTINES:* 

SOURCE  (2) 

TNEFLD  (3) 

TNHFLD  (3) 

♦2-GTD 

3-MOM 

916 

Ov>VV'.*'V,^V'  V. 

-  «  -  .  t*  •  9»  r  f  ,•  .*  -  •  '  «t'  • .  ’ 

1.  NAME:  ROTATE 


(GTD,  INPUT) 


2.  PURPOSE:  To  rotate  a  point  to  or  from  the  origin  of  a  given  coordi¬ 
nate  system. 

3.  METHOO:  The  coordinates  of  a  point  rotated  in  a  given  coordinate 
system  are  given  by 


X' 

Y* 

z' 


=  I  R_  I  •  R  •  R 
Z  y  x 


X 

Y 

Z 


where 


r  i 

10  0 

1  Rxl  * 

0  cosY  -sinY 

l  *J 

0  sinY  cosY 

m  «■ 

r  ' 

cosd  0  sin0 

I  R  |  * 

0  10 

L  yJ 

-sin0  0  cos0 

[■■] 


cos$ 

sin$ 

0 


-sin$  0 
cos$  0 
0  1 


and  Y,  0,  4>  are  the  rotation  angles  about  the  x,  yt  and  z  axis 
respectively. 

4.  INTERNAL  VARIABLES: 


VARIABLE 

DEFINITION 

A 

nra 

The  components  of  the  rotation  matrix 

CP 

cos$ 

CS 

cosY 

CT 

COS0 

NOP 

Operation  code  to  designate  rotation  from 

or  to  the  origin 


PHI 


the  rotation  about  the  z  axis 
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PHISV 

Saved  value  of  PHI 

PSI 

¥,  the  rotation  about  the  x  axis 

PSISV 

Saved  value  of  PSI 

RX,RY,RZ 

Rotation  angle  about  the  x,  y,  and  z  axis, 
respectively 

SP 

s1n<j» 

SS 

sin¥ 

ST 

sine 

THETA 

9,  the  rotation  about  the  y  axis 

THTSV 

Saved  value  of  THETA 

X.Y.Z 

Coordinates  of  Input/output  variables  to  be 
changed 

I/O  VARIABLES 

A.  INPUT 

LOCATION 

ISOFF 

/AOEBUG/ 

NOP 

F.P. 

RX.RY.RZ 

F.P. 

X.Y.Z 

F.P. 

B.  OUTPUT 

LOCATION 

X.Y.Z 

F.P. 

6.  CALLING  ROUTINES:* 
COORDS  (1) 

CYAXIS  (2) 

GTDCS  (1) 

PATCH  (1) 

WYRDRV  (1) 


*1- INPUT 
2-GTD 


ROTATE 


(GTD,  INPUT) 


CALLED  ROUTINES: 

ASSIGN 

STATIN 

STATOT 

WLKBCK 


COMPUTE  MATRIX 
ELEMENTS 


CALCULATE 

TRANSFORMED 

COORDINATES 


1ALCULATE  INVERSE 
TRANSFORMED 
COORDINATES 


NAME:  ROTRAN  (GTD) 

PURPOSE:  To  transform  and  rotate  a  point  or  vector  defined  In  the 
global  coordinate  system  (as  stored  In  the  segment  table  SEGTBL)  to 
the  cylinder-centered  reference  coordinate  system  (RCS)  used  for  the 
GTO  calculations. 

METHOD:  The  point  Yv  defined  In  the  global  coordinate  system  may  be 
represented  by  point  Xr^ In  the  cylinder-centered  RCS  where  (refer  to 
figure  1): 


Xrt  '  [Vcl  ]  V  *  X*  ‘  Xo 

XO(l) 
XO  ( 2 ) 
XO(3) j 


XRT ( 1 )  XCL( 1 
XRT( 2 )  *  YCL( 1 
XRT(3)J  [zCL(l 


)  XCL( 2)  XCL( 3 ) 
YCL(l)  YCL( 2 )  YCL( 3 ) 
ZCL(l)  ZCL( 2)  ZCL(3) 


where  Y0  Is  the  location  of  the  cylinder -centered  RCS  origin  defined 
in  the  global  coordinate  system  and  x,  y,  2  are  unit  vectors 
defining  the  cylinder  centered  RCS  axes  in  global  coordinate  system 
(g)  components: 


x  *  *g  XCL(l)  +  yg  XCL( 2 )  +  zg  XCL(3) 

A  A  A  A 

Y  *  *g  YCL(l)  +  yg  YCL( 2 )  +  zg  YCL(3) 

A  A  A  A 

z  -  *g  ZCL(l)  +  yg  ZCL( 2 )  +  zg  ZCL{3). 
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POINT  TO  BE 
TRANSFORMED 


NEW  COORDINATE 
SYSTEM  ORIGIN  AND 
AXES 

- 


OLD  COORDINATE  SYSTEM 
ORIGIN  AND  AXES 


*x  -  *old  +  fold  XX<2>  +  5old  XX<3> 

=  x  XRT(l)  +  y  XRT(2)  +  z  XRT(3) 

Figure  1.  Illustration  of  Old  and  New  Reference 
Coordinate  Systems 


INTERNAL  VARIABLES: 
VARIABLE 


DEFINITION 

This  defines  the  cylinder-centered  refer¬ 
ence  coordinate  system  x  axis  unit  vector 
in  global  system  components 

X,Y,  and  Z  components  of  the  cylinder- 
centered  reference  coordinate  system  origin 
location  defined  In  the  global  coordinate 
system 

X,Y,  and  Z  components  of  point  location  in 
RCS 

X,Y,  and  Z  components  of  point  location 
after  translating  global  coordinate  system 
origin  to  point  XO 

X,Y,  and  Z  components  of  point  location  in 
global  coordinate  system 

This  defines  the  cylinder-centered  refer¬ 
ence  coordinate  system  y  axis  in  global 
reference  system  components 
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This  defines  the  cylinder-centered  refer¬ 
ence  coordinate  system  z  axis  In  global 
reference  system  components 


I/O  VARIABLES: 

A.  INPUT 

LOCATION 

XCL 

/ROTRDT/ 

XO 

F.P. 

XX 

F.P. 

YCL 

/ROTRDT/ 

ZCL 

/ROTRDT/ 

B.  OUTPUT 

LOCATION 

XRT 

F.P. 

CALLING  ROUTINE: 

GTDORV 

CALLED  ROUTINE: 

NAME:  RPLDPL 


(6TD) 


PURPOSE:  To  calculate  the  unobstructed  electric  field  for  a  unit 
source  ray  that  is  reflected  off  plate  MR  and  diffracted  by  edge  ME 
on  plate  MP  into  a  given  far-field  observation  direction  or  to  a 
given  near-field  observation  point. 

METHOD:  RPLDPL  is  the  driver  routine  to  compute  the  plate-reflected 
and  then  edge-diffracted  fields.  Pertinent  geometry  is  shown  in 
figure  1  and  computation  details  are  given  in  references  A-C.  The 
fields  are  first  Initialized  to  zero.  Then  the  diffraction  point  on 
edge  ME  of  plate  MP  is  found  for  the  given  observation  direction  or 
point  based  on  the  location  of  the  source  Imaged  through  plate  MR. 
The  diffraction  point  is  computed  in  subroutine  DFPTWD.  If  diffrac¬ 
tion  did  not  occur,  debug  Information  (if  requested)  is  printed  on 
file  LUBRNT  and  control  Is  returned  to  the  calling  routine.  If 
diffraction  did  occur  on  the  vector  tangent  to  the  edge,  the  point 
location  is  checked  to  see  if  it  Is  between  the  corners.  If  it  is 
not,  the  diffraction  point  Is  set  at  the  closest  corner  and  a  flag, 
LDIF,  Is  set  to  indicate  only  corner-diffracted  fields  exist. 


/ 

f 

/ 

t 

•''SOUWCE  IMAGE 

location  xis 

Figure  1.  Illustration  of  a  Ray  Reflected  by  a  Plate 
and  Then  Diffracted  by  a  Plate  Edge 


Now  the  ray  path  Is  checked  for  any  obstructions.  If  the  path  is 
shadowed,  debug  information  (If  requested)  is  printed  and  control 
returns  to  the  calling  routine.  If  the  path  is  clear  the  field 
computations  begin. 
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First,  necessary  diffraction  angles  and  geometry  are  calculated. 
Then  the  source  field  pattern  factor  is  found  by  calling  subroutine 
SOURCE  and  is  multiplied  by  the  reflection  coefficient.  Then  the 
incident  field  perpendicular  and  parallel  to  the  edge  can  be 
computed. 

If  slope  diffraction  is  requested,  the  incident  slope  field  pattern 
factor  is  computed  by  calling  subroutine  SOURCP.  This  factor  also 
must  be  multiplied  by  the  reflection  coefficient  to  account  for  axis 
direction  changes  due  to  single  reflection. 

The  phase  term  is  now  computed.  The  edge  diffraction  coefficient  is 
determined  in  subroutine  DW.  Now  the  edge-diffracted  fields  are 
computed,  first  in  terms  of  a  parallel  and  perpendicular  orientation 
and  then  in  theta  and  phi  components.  By  calling  subroutine  XYZFLD 
the  x,  y,  z  components  of  the  edge-diffracted  field  are  computed  and 
accumulated  with  all  other  fields  computed  by  other  reflection  and 
diffraction  interactions. 

If  corner  diffraction  was  requested,  the  far-field  corner-diffracted 
fields  are  computed  for  each  corner  on  edge  ME  in  the  same  manner  as 
for  the  edge  diffraction. 

After  all  fields  have  been  computed  and  the  x,  y,  z  components 
accumulated,  debug  information  (if  requested)  is  printed  on  file 
LUPRNT.  This  consists  of  the  field  magnitude  and  theta  and  phi 
components  of  the  total  field.  Control  is  then  returned  to  the 
calling  routine. 


INTERNAL  VARIABLES: 

VARIABLE 

AON 

AFN 

BETN 

BETP 

BO 


DEFINITION 

Dot  product  of  vector  from  plate  MP  to  the 
source  image  and  the  plate  unit  normal 

Wedge  angle  number 

Difference  in  diffracted  and  incident  phi 
angles 

Sum  of  diffracted  and  incident  phi  angles 

Diffracted  field  beta  polarization  unit 
vector  in  diffraction  edge-fixed  coordinate 
system  (in  x,y,z  RCS  components) 


BOP 


Incident  field  beta  polarization  unit 
vector  In  diffraction  edge-fixed  coordinate 
system  (in  x,y,z  RCS  components) 
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BRD 

Lower  and  upper  limit  for  edge  diffraction 
angle 

CHP 

Cosine  of  half  wedge  angle 

CORN 

Corner  diffraction  coefficient 

CPH 

Cosine  of  PSR 

CPHJ 

Cosine  of  PH JR 

CPHO 

Cosine  of  PSOR 

CTH 

Cosine  of  THR 

CTHJ 

Cosine  of  THJR 

CTHP 

Cosine  of  THPR 

DEL 

Parameter  used  in  transition  function 

DH 

Diffraction  coefficient  for  hard  boundary 
condition 

OHIR 

Distance  from  reflection  point  to  diffrac¬ 
tion  point 

DHIT 

Distance  from  source  to  reflection  point 
(from  PLAINT) 

OHT 

Distance  from  source  to  hit  point  (from 
PLAINT  and  CYLINT) 

DIN 

Edge  diffraction  coefficient  (from  subrou¬ 
tine  DICOEF)  for  incident  diffracted  field 

DIP 

Edge  diffraction  coefficient  (from  subrou¬ 
tine  DICOEF)  for  reflected  diffracted  field 

DPH 

Slope  diffraction  coefficient  for  hard 
boundary  condition 

OPS 

Slope  diffraction  coefficient  for  soft 
boundary  condition 

DS 

Diffraction  coefficient  for  soft  boundary 
condition 
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EIPLP 


EIPRP 


EIX,EIY,EIZ 


Dot  product  of  edge  unit  vector  and 
diffracted  ray  propagation  direction  unit 
vector 

Phi  component  of  corner-diffracted  E-field 

Theta  component  of  corner-diffracted 
E-field 

Phi  component  of  edge-diffracted  E-field 

Component  of  diffracted  field  parallel  to 
the  edge 

Component  of  diffracted  field  perpendicular 
to  the  edge 

Theta  component  of  edge-diffracted  E-field 

Theta  component  of  corner-diffracted 
E-field  in  RCS 

Phi  component  of  corner-diffracted  E-field 
in  RCS 

Component  of  incident  field  parallel  to  the 
edge 

Pattern  factor  for  component  of  incident 
slope  field  parallel  to  the  edge 

Component  of  incident  field  perpendicular 
to  the  edge 

Pattern  factor  for  component  of  incident 
slope  field  perpendicular  to  the  edge 

Source  pattern  factor  for  x,y,  and  2  compo¬ 
nents  of  incident  E-field 

Complex  phase  term 

Wedge  angle  number 

Wedge  angle  indicator 

Angle  exterior  to  wedge  in  degrees 
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Dot  product  of  the  diffracted  ray  direction 
and  the  vector  from  the  origin  to  the 
diffraction  point 

Sign  change  variable 

Index  variable 

Logical  variable  set  true  if  diffraction 
point  is  on  edge  tangent  but  not  within 
corners.  (The  diffraction  is  set  to  the 
closest  corner) 

Logical  variable  set  true  if  edge  diffrac¬ 
tion  exists  (from  subroutine  DFPTWD) 

Set  true  if  ray  hits  a  plate  or  cylinder 
(from  PLAINT  or  CYLINT) 

Index  variable  used  to  step  through  corners 

Edge  on  plate  MP  where  diffraction  occurs 

Corner  at  end  of  edge  ME 

Plate  for  which  diffraction  occurs 

Plate  where  reflection  occurs 

DO  loop  variable 

DO  loop  variable 

DO  loop  variable 

Dot  product  of  diffraction  edge  bi normal 
and  diffracted  ray  propagation  direction 

Diffracted  field  phi  polarization  unit 
vector  in  diffraction  edge-fixed  coordi¬ 
nate  system  (in  x,y,z  RCS  components) 

Phi  component  of  reflected  ray  propagation 
direction  in  RCS 

Phi  component  of  incident  (source)  ray 
propagation  direction 
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PHO 

Incident  field  phi  polarization  unit  vector 
in  diffraction  edge-fixed  coordinate  system 
(in  x,y,z  RCS  components) 

PHSR 

Phi  component  of  ray  propagation  direction 
after  diffraction  in  RCS 

PP 

Negative  dot  product  of  diffraction  edge 
binormal  and  incident  ray  unit  vector 

PS 

PSR*DPR 

PSD 

Diffracted  ray  phi  angle  in  edge-fixed 
coordinate  system 

PSO 

PS0R*DPR 

PSOD 

Incident  ray  phi  angle  in  edge-fixed  coor¬ 
dinate  system 

PSOR 

Phi  component  of  incident  ray  direction  in 
edge-fixed  coordinate  system 

PSR 

Phi  component  of  diffracted  ray  propagation 
direction  in  edge-fixed  coordinate  system 

QO 

Dot  product  of  diffraction  plate  normal  and 
diffracted  ray  propagation  direction 

QI 

Negative  of  dot  product  of  diffracted  plate 
normal  and  incident  ray  propagation  iirec- 
tion 

SBO 

Sine  of  BO,  the  angle  the  diffracted  ray 
makes  with  the  edge 

SNF 

Distance  between  diffraction  point  and 
near-field  observation  point 

SNP 

Sine  of  half  wedge  angle 

SP 

Distance  from  source  image  to  diffraction 
point  (from  subroutine  DFPTWD) 

SPH 

Sine  of  PSR 

SPHJ 

Sine  of  PH JR 

SPHO 

Sine  of  PSOR 
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SPP 

STHJ 

STHR 

TERM 

THIR 

THJR 

THPR 

THR 

TPP 

VAX 

VC 

VCM 

VECT 

VI 

VIP 

VJ 

VMG 

XI 


^-Distance  from  source  image  to  modified 
diffraction  point 

Sine  of  THJR 

Sine  of  THR 

Coefficient  of  corner-diffracted  fields 

Theta  component  of  reflected  ray  direction 
in  RCS 

Theta  component  of  incident  (source)  ray 
propagation  direction 

Angle  diffracted  ray  makes  with  edge 

Angle  between  edge  unit  vector  and  ray  from 
source  image  location  to  corner  MC 

Distance  parameter  used  in  calculating 
diffraction  coefficients 

3X3  matrix  defining  the  source  image 
coordinate  system  axes 

Unit  vector  from  source  image  to  corner  1 
or  2  of  edge  ME 

Distance  from  source  image  to  corner  1  or  2 
of  edge  ME 

Vector  used  to  move  diffraction  point  off 
edge  for  shadowing  tests 

Unit  vector  of  ray  Incident  on  edge  from 
plate  reflection  (from  subroutine  DFPTWD) 

Unit  vector  of  ray  from  source  Image  to 
modified  diffraction  point 

X,Y,  and  Z  components  of  source  ray  propa¬ 
gation  direction 

Distance  along  the  edge  from  first  corner 
of  edge  to  diffraction  point 

Single  reflection  source  image  location 
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XC 

XO 

XOl 

XDP 


XS 

XSS 

ZP 

I/O  VARIABLES: 
A.  INPUT 
D 

DP 

DPR 

DT 

FLDPT 

FNN 

LCORNR 

LDEBUG 

LNRFLD 

LSLOPE 

LSURF 


Corner  location 

Diffraction  point  (calculated  in  subroutine 
DFPTWD)  in  RCS 

Diffraction  point  location 

Modified  diffraction  point  used  for 
shadowing  tests 

Source  image  location  (for  reflection  from 
plate  MR) 

Source  location  in  RCS 

Single  reflection  source  image  location 

Dot  product  of  propagation  direction  unit 
vector  and  vector  from  diffraction  point  to 
corner  MC 

LOCATION 

/DIR/ 

/THPHUV/ 

/PIS/ 

/THPHUV/ 

/NEAR/ 

F.P. 

/LOGDIF/ 

/TEST/ 

/NEAR/ 

/LOGDIF/ 

/SURFAC/ 

/ADEBUG/ 


LUPRNT 


RPLDPL 


HE 

MEP 

HP 

HR 

PHSR 

PI 

THSR 
TP  I 
V 

VHAG 

VN 

VP 

VXI 

X 

XI 
XS 

B.  OUTPUT 
ECPH 
ECTH 
EOPH 
EDTH 

CALLING  ROUTINE: 
GTDDRV 


»'  -V 


F.p. 

/GEOPLA/ 

F.P. 

F.P. 

/OIR/ 

/PIS/ 

/DIR/ 

/PIS/ 

/GEOPLA/ 

/EDHAG/ 

/GEOPLA/ 

/GEOPLA/ 

/IHAINF/ 

/GEOPLA/ 

/IHAINF/ 

/SORINF/ 

LOCATION 

F.P. 

F.P. 

F.P. 

F.P. 
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CALLED  ROUTINES: 

ASSIGN 

BEXP 

BTAN2 

CYLINT 

DFPTWD 

DICOEF 

DW 

FFCT 

NFD 

PLAINT 

REFBP 

SMAGNF 

SOURCE 

SOURCP 

STATIN 

STATOT 

TPNFLD 

WLKBCK 

XYZFLD 

REFERENCES: 
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Surface,"  Proc.  IEEE,  Vol.  62,  November  1974,  pp.  1448-1461. 


RPLDPL 


(GTD) 


W.  0.  Burnside  and  P.  H.  Pathak,  "A  Corner  Diffraction  Coeffi¬ 
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Y.  M.  Hwang  and  R.  G.  Kouyoumjlan,  "A  Dyadic  Diffraction 
Coefficient  for  an  Electromagnetic  Wave  Which  Is  Rapidly 
Varying  at  an  Edge,"  USNC-URSI  1974  Annual  Meeting,  Boulder, 
CO.,  Oct.  1974. 
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ANGLES 


COMPUTE  PHASE 


COMPUTE  OIFFRACTION 
POLARIZATION  UNIT 
VECTORS 


COMPUTE  EDGE  DIF¬ 
FRACTION  COEFFICIENTS 


COMPUTE  SINE  OF 
OIFFRACTION  ANGLE 


COMPUTE  SOURCE  FIELD 
PATTERN  FACTOR 

SOURCE 


MULTIPLY  SOURCE 
FACTOR  SY  REFLEC¬ 
TION  COEFFICIENT 


COMPUTE  COMPONENTS 
OF  INCIDENT  FIELD 
PERPENDICULAR  ANO 
PARALLEL  TO  THE 
EDGE 


IF  SLOPE  DIFFRACTION  IS 
OESIRED.  COMPUTE 
INCIDENT  SLOPE  FIELO 
PATTERN  FACTORS 
SOURCE 


MULTIPLY  SLOPE  FAC¬ 
TORS  BY  REFLECTION 
COEFFICIENT 
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FIELD 
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NAME:  RPLRCL  (GTO) 


PURPOSE:  To  compute  the  unobstructed  electric  field  from  a  unit 
source  reflected  by  plate  MP  and  then  reflected  by  the  cylinder  into 
the  far-field  observation  direction  or  to  the  near-field  observation 
point. 

METHOD:  RPLRCL  functions  as  a  service  routine  for  subroutine 
RPLSCL,  where  the  actual  plate-cylinder  fields  are  computed.  The 
geometrical  optics  reflected  field  components  ETH  and  EPH  computed 
in  RPLRCL  are  used  only  for  reference  purposes.  The  field  compo¬ 
nents  calculated  in  RPLRCL,  which  are  used  in  RPLSCL,  are  the  hard 
and  soft  components  of  the  plate-reflected  field  incident  on  the 
cylinder  at  the  reflection  point.  These  components,  along  with 
several  other  useful  parameters,  are  passed  to  subroutine  RPLSCL 
through  common  block  FUDGI.  Pertinent  geometry  is  shown  in 
figure  1. 


SOURCE  IMAGE 
LOCATION 


Figure  1.  Illustration  of  Plate-Reflected, 
Cylinder-Reflected  Ray 
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The  code  first  makes  two  checks  to  determine  if  it  is  possible  for 
reflection  to  occur  off  the  cylinder.  One  check  determines  if  the 
observation  direction  or  point  is  in  the  lit  or  dark  geometrical 
optic  region  of  the  cylinder.  The  other  check  is  to  determine  if 
the  observation  direction  or  point  is  in  the  paraxial  region  beyond 
the  end  caps.  If  reflection  cannot  occur,  a  flag  is  set  to  indicate 
starting  data  are  not  available  for  the  next  time  subroutine  RPLRCL 
Is  called.  Another  flag  is  also  set  to  indicate  the  field  does  not 
exist.  The  fields  are  set  to  zero  and  control  is  returned  to  the 
calling  routine. 

If  cylinder  reflection  could  occur,  the  reflection  point  is  deter¬ 
mined  by  calling  subroutine  RFDFIN  for  near-field  calculations  and 
subroutine  RFPTCL  for  far-field  calculations.  The  code  then  checks 
to  see  if  the  cylinder  reflection  point  is  beyond  the  cylinder  end 
caps.  If  It  is,  a  flag  is  set  to  indicate  that  the  geometrical 
optics  field  does  not  exist.  The  fields  are  set  to  zero  and  control 
is  returned  to  the  calling  routine.  If  the  reflection  point  is  on 
the  curved  surface  of  the  cylinder,  the  ray  path  from  the  cylinder 
reflection  point  in  the  far-field  observation  direction  or  to  the 
near-field  observation  point  is  checked  for  obstructions.  Then,  the 
code  checks  to  see  if  reflection  occurs  from  plate  MP.  Next,  the 
remainder  of  the  complete  ray  path  is  checked  for  obstructions.  If 
at  any  point  the  ray  path  is  shadowed  or  reflection  does  not  occur 
from  plate  MP,  the  code  sets  a  flag  to  indicate  that  the  geometrical 
optics  field  does  not  exist.  Those  fields  are  set  to  zero  and 
control  is  returned  to  the  calling  routine. 

If  the  plate  and  cylinder  reflections  do  occur,  and  the  ray  path  is 
not  obstructed,  the  field  computations  can  begin.  First,  the  source 
field  pattern  factor  is  found  by  calling  subroutine  SOURCE.  The 
source  factor  is  then  multiplied  by  the  reflection  coefficient. 
Next,  the  polarization  unit  vectors  perpendicular  and  parallel  to 
the  plane  of  incidence  are  computed.  These  are  used  to  compute  the 
Incident  field  components  parallel  and  perpendicular  to  the  plane  of 
incidence.  The  plate-reflected  field  components  parallel  and 
perpendicular  to  the  plane  of  incidence  are  then  computed.  Next, 
the  theta  and  phi  components  of  the  reflected  field  are  computed. 
The  subroutine  then  ends  by  computing  the  theta  and  phi  components 
of  the  hard  and  soft  field  components  incident  upon  the  cylinder. 


INTERNAL  VARIABLES: 
VARIABLE 


DEFINITION 

Dot  product  of  cylinder  normal  and  reflec¬ 
tion  propagation  direction  unit  vector 

Cosine  of  WR 


>r*  r 


Propagation  direction  after  cylinder 
reflection  in  x,  y,  z  RCS  components 

Direction  unit  vector  used  for  determining 
if  cylinder  reflection  is  possible 

Direction  unit  vector  used  for  determining 
if  cylinder  reflection  is  possible 

X-Y  distance  from  the  z  axis  to  the 
cylinder  reflection  point 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  1  (2-D) 

Dot  product  of  unit  vector  of  propagation 
direction  and  cylinder  tangent  unit  vector 
through  tangent  point  2  (2-D) 

Distance  from  reflection  point  on  plate  to 
reflection  point  on  the  cylinder 

Distance  from  source  to  hit  point  (from 
PLAINT) 

Distance  to  hit  point  (from  subroutine 
PLAINT) 

X,Y,  and  Z  components  of  incident  ray 
direction  on  cylinder  in  RCS 

Variable  used  to  indicate  if  reflection 
occurred  and  satisfied  Snell's  Law 

X,Y,Z  components  of  propagation  direction 
of  ray  incident  on  plate 

Dot  product  of  vector  from  origin  to  source 
image  location  and  propagation  direction 
(2-D) 

Pattern  factor  of  theta  component  of  inci¬ 
dent  field  in  RCS 

Pattern  factor  of  phi  component  of  incident 
field  in  RCS 


RPLRCL 


(GTD) 


ehphi 


EHTHI 


ERX,ERY,ERZ 


ESPHI 


ESTHI 
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FPTXY 


Phi  component  of  the  hard  component  of 

field  Incident  on  cylinder  (parallel  to 

plane  of  incidence) 

Theta  component  of  the  hard  component  of 

field  incident  on  cylinder  (parallel  to 

plane  of  incidence) 

Incident  cylinder  field  component  parallel 
to  plane  of  incidence 

Incident  cylinder  field  component  perpen¬ 
dicular  to  plane  of  incidence 

Phi  component  of  cylinder-reflected  E-field 

Cylinder-reflected  field  component  parallel 
to  plane  of  incidence 

Cylinder-reflected  field  component  perpen¬ 
dicular  to  plane  of  incidence 

X,Y,Z  components  of  field  incident  on  (or 
reflected  from)  cylinder  in  RCS 

Phi  component  of  the  soft  component  of 
field  incident  on  cylinder  (perpendicular 
to  plane  of  incidence) 

Theta  component  of  the  soft  component  of 
field  Incident  on  cylinder  (perpendicular 
to  plane  of  Incidence) 

Theta  component  of  cylinder-reflected  E- 
field 

Pattern  factor  of  x,y,z  components  of 
source  field  incident  on  cylinder  in  RCS 

Location  of  field  point  in  z=0,  x-y  plane 

Set  true  if  ray  hits  plate  (from  PLAINT) 

Set  true  if  reflection  data  are  available 
from  previous  pattern  angle  (or  for  next 
pattern  angle  (when  leaving  routine)) 
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LTRFI 

MP 

ORIGIN 

PH 

PHIR 

PHJR 

PHSR1 

PHSR2 

RHOII 

RH02 

S 

SMAGI 

SNF 

SQRH 

SXN.SYN.SZN 

THIR 

THJR 


Set  true  if  geometrical  optics  reflected- 
reflected  field  does  not  exist 

Plate  on  which  reflection  occurs 

The  origin  of  the  reference  coordinate 
system  (RCS)  (0.,  0.,  0.) 

Complex  phase  and  ray  spreading  coefficient 

Phi  component  of  incident  ray  direction  on 
cylinder 

Phi  angle  for  direction  of  ray  incident  on 
plate 

Phi  angle  of  D1 
Phi  angle  of  D2 

Ray  spreading  radius  in  plane  of  cylinder 
curvature  at  reflection  point 

Ray  spreading  radius  normal  to  plane  of 
incidence  at  reflection  point 

Oi stance  from  source  image  to  cylinder 
reflection  point 

Length  of  ray  from  reflection  point  on 
cylinder  to  source  image 

Distance  between  field  point  and  cylinder 
reflection  point 

Part  of  spr  :ading  factor 

X,Y,  and  Z  components  of  unit  vector  of 
ray  from  reflection  point  on  cylinder  to 
source  image  location  in  RCS 

Theta  component  of  incident  ray  direction 
on  cylinder 

Theta  angle  which  defines  direction  of  ray 
incident  on  plate 


THSR1 


Theta  angle  of  D1 
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THSR2 


UIPPX,UIPPY,UIPPZ 


UIPRX,UIPRY,UIPRZ 


URPPX ,URPPY , URPPZ 


Theta  angle  of  D 2 

Unit  vector  of  binonnal  to  cylinder  at 

reflection  point  (2-D) 

Unit  vector  of  normal  to  cylinder  at 

reflection  point  (2-0) 

X,Y,Z  components  of  incident  field  polar¬ 
ization  unit  vector  parallel  to  plane  of 

incidence 

X,Y,Z  components  of  incident  reflected 
field  polarization  unit  vector  perpendic¬ 
ular  to  plane  of  incidence 

X,Y,Z  components  of  reflected  field 
polarization  unit  vector  parallel  to  plane 
of  incidence 

Matrix  defining  source  coordinate  system 
axes  in  RCS  components 

Phi  angle  at  which  cylinder  reflection 
occurs 

Vector  in  the  direction  of  the  more  posi¬ 
tive  end  cap  used  to  determine  if  cylinder 
reflection  can  occur 

Vector  in  the  direction  of  the  more  nega¬ 
tive  end  cap  used  to  determine  If  cylinder 
reflection  can  occur 

X,Y,Z  components  of  source  image  location 
also  reflection  point  on  plate 

Location  of  reflection  point  on  cylinder  In 
x.y ,z  RCS 

Source  image  location 

Cylinder  tangent  point  one  location  in  x-y 
plane  of  vector  from  source  image 

Cylinder  tangent  point  two  location  in  x-y 
plane  of  vector  from  source  image 
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I/O  VARIABLES: 

A.  INPUT 

LOCATION 

A 

/GEOMEL/ 

B 

/GEOMEL/ 

BTI 

/BNDICL/ 

CPS 

/DIR/ 

CTC 

/GEOMEL/ 

CTHS 

/DIR/ 

0 

/DIR/ 

DP 

/THPHUV/ 

DT 

/THPHUV/ 

DTI 

/BNDICL/ 

FLDPT 

/NEAR/ 

LDEBUG 

/TEST/ 

LNRFLD 

/NEAR/ 

LRFI 

/CLRFI/ 

LUPRNT 

/ADEBUG/ 

HP 

F.  P. 

PHSR 

/DIR/ 

PI 

/PIS/ 

SPS 

/DIR/ 

STHS 

/DIR/ 

THSR 

/DIR/ 

TP  I 

/PIS/ 

VTI 

/BNDICL/ 

CALLEO  ROUTINES: 

ASSIGN 

RFOFIN 

BEXP 

RFPTCL 

BTAN2 

SMAGNF 

CYLINT 

SOURCE 

NANDB 

STATIN 

NFD 

STATOT 

PLAINT 

TPNFLD 

REFBP 

WLKBCK 
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PURPOSE:  To  compute  the  unobstructed  electric  field  from  a  unit 
source  due  to  double  reflection  from  specified  plates  (reflection 
off  plate  MP  and  then  plate  MPP). 

METHOD:  RPLRPL  Is  the  driver  routine  which  directs  the  ray  tracing, 
physics  and  field  calculations  for  double  reflection  off  specified 
plates  In  a  given  far-field  direction  or  to  a  near-field  observation 
point  from  a  unit  source.  The  pertinent  geometry  Is  shown  in 
figure  1. 


PLATE  MPP 


SOURCE  XS 


A.  owls'! 

DHIT  // 


/  \ 


n  DHIJ 

W 

\ 


9 

V 


plate  mp 


SINGLE  REFLECTION  DOUBLE  REFLECTION 
SOURCE  IMAGE  SOURCE  IMAGE 

Xis  XIJ 


Figure  1.  Geometry  for  Double  Reflected  Ray 


First  the  ray  path  from  the  double  reflection  source  image  is 
checked  to  see  If  It  passes  through  plate  MPP.  If  It  does,  the  ray 
path  from  the  reflection  point  on  plate  MPP  In  the  desired  far-field 
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direction  or  to  the  near-field  observation  point  is  checked  for 
obstructions.  If  this  ray  did  not  hit  another  plate  or  a  cylinder, 
the  rays  between  the  first  and  second  reflection  points  and  between 
the  source  and  first  reflection  point  are  also  checked  to  see  if 
they  are  blocked.  If  none  of  these  separate  paths  is  blocked,  then 
it  is  known  that  reflections  on  the  two  plates  specified  did  occur 
in  the  correct  order  and  that  the  complete  ray  path  is  unobstructed. 
If  at  any  check  the  ray  was  blocked  by  another  plate  or  a  cylinder, 
or  the  ray  did  not  pass  through  plates  MP  and  MPP  as  required  for 
double  reflection,  the  code  immediately  sets  the  theta  and  phi 
components  of  the  electric  field  to  zero,  and  no  other  computations 
except  debug  functions  (if  requested)  are  performed  in  this  routine. 

If  the  ray  path  is  unobstructed  and  the  reflections  occurred,  the 
sourca  field  pattern  factor  at  the  double  reflection  source  image 
location  is  computed  in  subroutine  SOURCE.  The  phase  factor  is  then 
computed.  For  far  field  this  will  refer  the  electric  field  back  to 
the  origin  of  the  reference  coordinate  system  (RCS).  For  near  field 
the  phase  factor  includes  the  spherical  wave  spread  factor.  Now  the 
theta  and  phi  components  of  the  electric  field  are  computed.  The 
electric  field  is  given  as: 


(EF  9 


EG  $) 


,j2*(XU  *  D>  for  far 


field 


from  subroutine  from  sub-  EX  -  phase  factor 
SOURCE  -  routine 

SOURCE  - 

theta  component  phi  component 
of  source  of  source 

factor  factor 


and 


E  =  (EF  8 


+  EG  <j>) 


- j2nSNF 
e  J _ 

SNF 


for  near  field. 


from  subroutine  from  subroutine 
SOURCE  -  SOURCE  - 

theta  component  phi  component 
of  source  of  source 

factor  factor 


EX  -  phase  factor 
where  SNF  =  I FlDPT  -  XIJ I 
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The  x,  y,  z  components  of  the  electric  field  are  computed  and 
accumulated  by  calling  subroutine  XYZFLD. 

If  the  debug  capabilities  have  been  requested,  the  doubly  reflected 
field  magnitude  is  computed.  The  magnitude,  theta  and  phi  complex 
components  of  the  field  are  printed  on  file  LUPRNT. 

4.  INTERNAL  VARIABLES: 


VARIABLE 

CPHI 

CPHJ 

CTHI 

CTHJ 


DEFINITION 
Cosine  of  PHIR 
Cosine  of  PHJR 
Cosine  of  THIR 
Cosine  of  THJR 


D 


X,Y,Z  components  of  ray  propagation  direc¬ 
tion  after  second  reflection  in  RCS 


DHIJ  Distance  from  double  reflection  image  to 

hit  point  on  plate  HPP 

DHIS  Distance  between  reflection  points 


DHIT  Distance  from  source  to  reflection  point 

(from  PLAINT) 

DHT  Distance  (calculated  in  PLAINT  or  CYLINT) 

from  source  or  point  from  which  ray 
originates  to  hit  point 

OICOEF  X,Y,Z  components  of  incident  ray  propaga¬ 

tion  direction  in  RCS 


DJ  X,Y,Z  components  of  propagation  direction 

of  ray  incident  on  plate  MPP 


EF 

Theta  component  of 
factor 

source 

field 

pattern 

EG 

Phi  component  of 
factor 

source 

field 

pattern 

EIX,EIY,EIZ 

X,Y,Z  components  of 
factor 

source 

field 

pattern 
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Phi- component  of  electric  field 

Theta  component  of  electric  field 

Complex  phase  factor 

The  electric  field  magnitude 

The  x,y,z  location  of  the  near-field  obser¬ 
vation  point 

Phase  distance  to  origin  (dot  product  of 
double  reflection  image  location  and 
reflected  ray  propagation  direction) 

Logical  variable  set  true  if  debug  option 
requested 

Set  true  if  ray  intersects  a  plate  or 
cylinder  (from  PLAINT  or  CYLINT) 

Flag  to  indicate  if  far-field  (LNRFLD=0)  or 
near-field  (LNRFLD=I)  calculations  are 
requested 

Plate  from  which  first  reflection  occurs 

Plate  from  which  second  reflection  occurs 

Phi  angle  of  incident  ray  propagation 
direction  in  RCS 

Phi  angle  of  ray  direction  between  reflec¬ 
tions  in  RCS 

Phi  angle  of  ray  propagation  direction 
after  reflection  in  RCS 


Distance  from  double  reflection  image  loca¬ 
tion  to  observation  field  point 

Sine  of  PHIR 

Sine  of  PH JR 

Sine  of  THIR 


Sine  of  THJR 
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THIR 

Theta  angle  of  incident  ray  propagation 
direction  in  RCS 

THJR 

Theta  angle  of  ray  direction  between 
reflections  in  RCS 

THSR 

Theta  angle  of  ray  propagation  direction 
after  reflections  in  RCS 

VAX 

X,Y,Z  components  defining  unit  vectors  of 
the  source  image  coordinate  system  axes  in 
RCS  components 

VAXP 

X,Y,Z  components  defining  unit  vectors  of 
the  source  image  coordinate  system  axes  in 
RCS  for  double  reflection 

XI 

Triply  dimensioned  array  of  image  locations 

XIJ 

X,Y,Z  components  of  double  reflection  image 
location 

XIS 

X,Y,Z  components  of  single  reflert'on 
source  image  location  (single  reflection 
from  plate  MP) 

XQ 

XtY,Z  components  of  double  reflection  image 
location 

xs 

Source  location  in  x,y,z  RCS 

I/O  VARIABLES: 

A.  INPUT  LOCATION 


0  /DIR/ 
FLDPT  /NEAR/ 
LDEBUG  /TEST/ 
LNRFLD  /NEAR/ 

/ADEBUG/ 


LUPRNT 


COMPUTE  00U8LE 
REFLECTED 
SOURCE  IMAGE  AXES 
IMDIR 


COMPUTE  SOURCE  FIELD 
PATTERN  FACTOR 


COMPUTE  PHASE  TERM 


COMPUTE  AND 
ACCUMULATE 
x,y,z  COMPONENTS  OF 


COMPUTE  FIELD 
.MAGNITUDE 


PRINT  MAGNITUDE, 
E-THETA,  E-PHI  ON 
FILE  LUPRNT 
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2.  PURPOSE:  To  compute  the  unobstructed  electric  field  from  a  unit 
source  reflected  by  plate  MP  and  then  scattered  by  the  cylinder  in 
the  given  far-field  observation  direction  or  to  a  given  near-field 
observation  point. 

3.  METHOD:  RPLSCL  is  the  driver  routine  which  directs  all  the  ray 
tracing,  physics  and  field  calculations  for  determining  the  electric 
field  reflected  by  a  plate  and  then  scattered  by  the  cylinder  in  the 
given  far-field  observation  direction  or  to  a  given  near-field 
observation  point.  The  geometry  is  shown  in  figure  1. 


Figure  1.  Illustration  of  Ray  Reflected  by  a  Plate  and 
Then  Scattered  by  the  Clyinder 


The  code  begins  by  initializing  the  fields  to  zero.  The  code  then 
makes  a  check  to  see  if  the  plate-reflected  rays  can  illuminate  the 
cylinder  curved  surface.  If  they  cannot,  debug  information  (if 
requested)  is  printed  on  file  LUPRNT.  Control  is  then  returned  to 
the  calling  routine.  If  the  plate-reflected  rays  can  illuminate  the 
curved  surface  of  the  cylinder,  the  code  steps  through  the  tangent 
vectors  to  calculate  fields  based  on  a  value  ALR.  ALR  is  the 
reflected  ray  phi  angle  in  the  tangent  point  coordinate  system.  The 
tangents  are  found  from  the  image  of  the  source  through  plate  MP 
tangent  to  the  cylinder  curved  surface.  If  ALR  is  less  than  it,  a 
plate-reflected  cylinder-reflected  ray  can  be  determined.  If  ALR  is 
greater  tha.  t,  a  creeping  wave  exists  for  this  tangent  on  the 
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cylinder.  If  ALR  is  approximately  tt,  grazing  incidence  occurs  on 
the  cylinder.  After  the  fields  associated  with  one  tangent  have 
been  found,  the  code  proceeds  to  the  next  tangent  and  calculates  the 
new  value  of  ALR.  The  fields  associated  with  this  tangent  are  then 
determined.  All  the  fields  are  accumulated  in  subroutine  XYZFLD  in 
common  block  FLDXYZ.  After  the  total  field  is  found,  debug  informa¬ 
tion  (if  requested)  is  printed  on  file  LUPRNT.  The  debug  informa¬ 
tion  consists  of  field  magnitude,  theta  and  phi  components  of  the 
total  field.  Control  is  then  returned  to  the  calling  routine. 

If  ALR  is  less  than  tt,  a  reflected  ray  path  will  possibly  be  found 
on  the  cylinder.  Subroutine  RPLRCL  calculates  the  ray  path  for  a 
plate-reflected  ray  followed  by  a  cylinder-reflected  ray.  It  also 
calculates  parameters  associated  with  this  ray  path  and  field. 
These  parameters  and  the  field  incident  upon  the  cylinder  are  passed 
to  subroutine  RPLSCL  through  common  block  FUDGI.  Once  control  is 
returned  to  RPLSCL,  this  subroutine  checks  to  make  sure  reflected 
fields  are  present.  If  they  are  not,  the  code  will  proceed  to  the 
next  tangent.  If  reflected  fields  are  present,  the  code  checks  to 
see  if  reflection  should  be  handled  by  the  second  tangent.  If  it 
should,  the  tangent  index  is  set  for  the  second  tangent.  Field 
computations  can  now  begin.  The  hard  and  soft  components  of  the 
field  incident  on  the  cylinder,  obtained  from  common  block  FUDGI, 
are  converted  into  the  cylinder-reflected  field.  This  is  the  total 
plate-reflected  cylinder-reflected  field.  The  x,  y,  z  components  of 
this  field  are  computed  in  subroutine  XYZFLD.  The  field  is  also 
accumulated  in  this  subroutine  with  fields  from  other  interactions. 

If  ALR  is  greater  than  n,  a  creeping  wave  can  occur  along  the 
cylinder.  The  code  to  determine  the  incident  point  and  the  point  it 
which  the  creeping  wave  leaves  the  cylinder  is  different  for  near¬ 
field  and  far-field  calculations.  Please  refer  to  the  accompanying 
flowchart  for  the  specifics  of  this  procedure.  While  in  the  differ¬ 
ent  near-field  and  far-field  paths,  the  code  does  compute  the  same 
items.  The  x,  y,  z  components  of  the  point  at  which  the  creeping 
wave  leaves  the  cylinder  and  the  x,y,z  components  of  the  point  at 
which  the  creeping  wave  begins  on  the  cylinder  are  computed.  These 
points  are  checked  to  make  sure  that  they  exist  on  the  curved  sur¬ 
face  of  the  cylinder  and  not  beyond  the  end  caps.  Also,  the  code 
checks  to  see  if  reflection  from  plate  MP  can  occur.  The  ray  path 
between  the  reflection  and  initial  diffraction  point  are  checked  for 
shadowing.  The  ray  path  from  the  source  to  the  plate  reflection 
point  is  checked  for  shadowing  also.  While  still  in  the  separate 
near-f ield/far-f ield  paths,  the  incident  source  field  pattern  factor 
is  computed  by  calling  subroutine  SOURCE.  The  source  factor  is 
multiplied  by  the  reflection  coefficient.  Various  other  parameters 
associated  with  the  path  lengths  needed  for  phase  factors  and  ray 
spreading  radii  are  computed.  The  code  comes  back  together  and 
checks  if  the  cylinder-diffracted  ray  is  obstructed.  If  not,  the 
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total  phase  factor  is  computed,  and  the  hard  and  soft  components  of 
the  creeping  wave  are  determined.  From  this,  the  total  plate- 
reflected  cylinder-diffracted  field  can  be  found.  It  is  converted 
into  x,y,z  components  and  accumulated  in  subroutine  XYZFLD. 

If  ALR  is  approximately  equal  to  tt,  grazing  incidence  can  occur  on 
the  cylinder.  In  this  section,  the  code  first  checks  to  see  if 
reflection  from  plate  MP  can  occur.  If  it  can,  the  ray  between  the 
reflection  point  and  the  far-field  observation  direction  or  near¬ 
field  observation  point  is  checked  for  obstructions.  Then  the 
source  ray  is  checked  for  shadowing.  If  the  ray  paths  are  clear, 
the  grazing  incident  point  is  computed.  This  point  is  checked  to 
make  sure  it  is  on  the  curved  surface  of  the  cylinder  and  not  beyond 
the  end  caps.  If  the  ray  path  is  legitimate,  then  the  source  field 
pattern  factor  is  computed  by  calling  subroutine  SOURCE.  The  source 
factor  is  multiplied  by  the  reflection  coefficient.  A  phase  factor 
is  determined.  The  hard  and  soft  components  of  the  field  incident 
at  the  grazing  point  are  determined.  By  combining  this  field  with 
the  grazing  incident  transition  field,  the  total  plate-reflected 
cylinder-scattered  field  can  be  determined.  The  x,y,z  components  of 
the  grazing  incidence  field  are  computed  and  accumulated  in  subrou¬ 
tine  XYZFLD. 

4.  INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

ALR  Cylinder-reflected  ray  phi  angle  in  tangent 

point  coordinate  system  (2-D) 

ALS  Phi  angle  defining  direction  of  ray  from 

RCS  origin  to  source  image  in  tangent  point 
coordinate  system 


v- 

% 


BX,BY,BZ 


Km 

CCC 

s,‘ 

i 

CFH 

i 

CFS 

DEPH 

DETH 

X,Y,Z  components  of  polarization  unit 
vector  of  soft  component  of  field  incident 
on  cylinder  (parallel  to  cylinder  surface 
and  normal  to  incident  ray  propagation 
direction) 

Real  part  of  the  Fresnel  integral 

Hard  transition  field  coefficient 

Soft  transition  field  coefficient 

Phi  component  of  transition  field  in  RCS 

Theta  component  of  transition  field  in  RCS 


DHIT 


Oi stance  from  source  image  to  plate  reflec¬ 
tion  point  (from  PLAINT) 


DHIV 

Distance  from  plate  reflection  point  to 
cylinder 

DHT 

Distance  from  source  to  hit  point  (from 
PLAINT) 

DICOEF 

Unit  vector  of  ray  incident  on  cylinder 

OIJ 

X,Y  plane  vector  from  a  source  image 
tangent  ray  to  the  point  the  creeping  wave 
leaves  the  cylinder. 

DIJXDJ 

Cross  product  of  DIJ  and  DJT 

DIT 

Cylinder  incident  ray  vector 

DIXDIJ 

Cross  product  of  DIT  and  DIJ 

DJ 

X,Y,Z  components  of  unit  vector  of  propaga¬ 
tion  direction  of  source  ray  incident  on 
plate 

DJT 

X-Y  plane  components  of  observation 
direction 

DMAG 

Distance  between  plate  reflection  point  and 
the  near-field  observation  point  for 
grazing  incidence  calculations 

EF 

Pattern  factor  for  theta  component  of 
Incident  field  in  RCS 

EG 

Pattern  factor  for  phi  component  of 
Incident  field  in  RCS 

EHP 

Phi  component  of  hard  component  of  field 
incident  on  cylinder  in  RCS 

EHT 

Theta  component  of  hard  component  of  field 
incident  on  cylinder  in  RCS 

EIX,EIY,EIZ 

Pattern  factor  for  x,y,z  components  of 
incident  field  in  RCS 

EP 

Phi  component  of  cylinder-scattered  E-field 
with  phase  referred  to  RCS  origin 
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Dot  product  of  unit  vector  tangent  to 
cylinder  and  the  propagation  direction  unit 
vector 

Phi  component  of  field  reflected  by  plate 
then  reflected  by  cylinder 

Theta  component  of  field  reflected  by  plate 
then  reflected  by  cylinder 

Phi  component  of  soft  component  of  field 
incident  on  cylinder  in  RCS 

Theta  component  of  soft  component  of  field 
incident  on  cylinder  in  RCS 

Theta  component  of  cylinder-scattered 
E-field  with  phase  referred  to  RCS  origin 

The  x-y  plane  components  of  the  near-field 
observation  field  point 

Variable  used  to  step  through  tangent 
points 

Set  true  if  ray  hits  a  plate  (from  PLAINT) 

(Returned  from  RPLRCL)  set  true  if  geo¬ 
metrical  optics  cylinder-reflected  field 
does  not  exist 

Logical  variable  set  true  first  time 
creeping  wave  computations  begin 

Plate  where  reflection  occurs 


RCS  origin  (0. ,0. ,0.) 

Phi  component  of  propagation  direction  of 
ray  incident  on  cylinder 

Phi  component  of  propagation  direction  of 
source  ray  incident  on  plate 

Length  of  vector  from  source  image  to 
tangent  point  (2  or  3-0).  Also,  the  total 
distance  between  the  source  image  point  and 
the  cylinder  incidence  point 


X-Y  plane  distance  between  the  source  image 
and  incident  point 

X-Y  plane  distance  between  the  field  point 
and  the  point  on  the  cylinder  at  which  the 
creeping  wave  leaves 

Distance  between  near-field  observation 
point  and  point  at  which  ray  path  leaves 
cylinder  (from  reflection  or  creeping  wave) 

Distance  of  path  along  the  cylinder 

Imaginary  part  of  the  Fresnel  integral 

Elliptical  angle  defining  the  source 
tangent  point  x-y  location 

Theta  component  of  propagation  direction  of 
ray  incident  on  cylinder 

Theta  component  of  propagation  direction  of 
source  ray  incident  on  plate 


Unit  binormal  at  reflection  point  phi  angle 
(2-0)  in  x-y  plane 

Unit  normal  at  reflection  point  phi  angle 
(2-0)  in  x-y  plane 

Source  image  axes 

Elliptical  angle  defining  point  where 
creeping  wave  leaves  cylinder 

Elliptical  angle  used  to  define  tangent 
points  (2-D) 

Elliptical  angles  defining  the  two  tangent 
points  on  the  cylinder  for  the  vector  from 
the  field  point  tangent  to  the  cylinder 

The  elliptical  angle  defining  the  x-y  plane 
point  on  the  cylinder  at  which  the  creeping 
wave  leaves  the  cylinder 

Elliptical  angle  defining  lower  range  of 
creeping  wave  travel  on  cylinder  (2-D) 
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RPLSCL 


(GTD) 


XD,YD,ZD 
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I/O  VARIABLES: 
l\.  INPUT 


Elliptical  angle  defining  upper  range  of 
creeping  wave  travel  on  cylinder  (2-D) 

X,Y,Z  components  of  direction  of  ray  from 

source  to  cylinder  tangent  point  (incident 
ray  for  creeping  and  grazing  incidence 
cases) 

X,Y,Z  components  of  point  where  incident 

creeping  wave  (or  grazing  wave)  meets 
cylinder 

X,Y,Z  components  of  image  source  location 

(for  reflection  from  plate  MP) 

X,Y,Z  components  of  point  where  ray  leaves 
cylinder 

Incident  point  on  cylinder 

X,Y,Z  components  of  point  where  creeping 

wave  leaves  cylinder 

Source  image  location  through  plate  MP 

Argument  of  the  Fresnel  integral 


LOCATION 

/GEOMEL/ 

/GTD/ 

/GEOMEL/ 

/BNDICL/ 

/COMP/ 

/COMP/ 

/DIR/ 

/GEOMEL/ 


'iSS 


RPLSCL 


(GTD) 


VXI 

/ I M A INI/ 

XI 

/IMAINI/ 

XRI 

/FUDGI/ 

XS 

/SORINF/ 

zc 

/GEOMEL/ 

B.  OUTPUT 

LOCATION 

EP 

F.P. 

ERP 

F.P. 

ERT 

F.P. 

ET 

F.P. 

LRFI 

/CLRFI/ 

CALLING  ROUTINE: 

GTDDRV 

CALLEO  ROUTINES: 

ASSIGN 

QFUN 

BEXP 

RADCV 

BTAN2 

REFBP 

CYLINT 

RPLRCL 

DQG32 

SMAGNF 

FCT 

SOURCE 

FKARG 

STATIN 

FRNELS 

STATOT 

NANDB 

TANG 

NFD 

WLKBCK 

PFUN 

XYZFLO 

PLAINT 


RPLSCL 


ACCUMULATE  TOTAL 
FIELDS 


RPLSCL  (GTD) 
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SWITCH  TO  LAST 
TANGENT 

6*  1 

f 

DETERMINE  x.y 
COMPONENTS  OF 

POINT  WHERE  CREEPING 
WAVE  LEAVES  CYLINOER 

1 

r 

DETERMINE 
PONENTS  OF 
CREEPING  W 

x.yCOM- 

WHERE 

AVE  BEGINS 

COMPUTE  t  LOCATIONS 
OF  WHERE  CREEP'NG 
WAVE  BEGINS  AND  OF 
WHERE  IT  LEAVES  THE 
CYLINDER 
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COMPUTE  PHASE 
FACTOR 


COMPUTE  HARO  AND 
SOFT  CREEPING  WAVE 
COMPONENTS 
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/‘<^source'v'sv 

RAY  SHADOWED 
BEFORE  HITTING 
^\PLAT£  MP^ 


SPECIFY  SOURCE  IMAGE 
AXES 


(COMPUTE  GRAZING 
INCIDENCE  FIELD) 
ALR  S  7T 


>^DOES\^ 
'REFLECTION^' 
FROM  PLATE 
\MP  OCCURS 


RAYS,. 

"SHADOWED^ 

BETWEEN 

REFLECTIONS 


COMPUTE  SOURCE  FIELD 
PATTERN  FACTOR 


/  IS  \ 

^  SOURCE  S 
RAY  SHADOWED. 


MULTIPLY  SOURCE  FAC¬ 
TOR  BY  REFLECTION 
COEFFICIENT 


CALCULATE  GRAZING 
INCIDENCE  POINT 


COMPUTE  POINT  WHERE 
CREEPING  WAVE  LEAVES 
CYLINDER 


.X' POINT 
'ON  CYLINDER  ^ 
CURVED  SURFACE 


-v 

v  > 


W-  ■ 
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NAME:  RWCOMS  (GTD,  INPUT,  MOM,  OUTPUT) 

PURPOSE:  To  read  or  write  common  blocks  on  the  checkpoint  file. 


METHOD:  RWCOMS  uses  arrays  the  same  length  as  each  common  area. 
These  common  arrays  can  be  read  from,  or  written  to,  the  checkpoint 
file. 


INTERNAL  VARIABLES: 

VARIABLES 

ADEBG 

AMPZJ 

ARGCM 

CSTM 

DFDT 

FLDCM 

GEODT 

GTDDT 

I  COM 

ICOMSV 

I  EOF 

IFILE 

INDXP1 

INDXWB 

INTM 

IOFLS 

IWBSAV 

JNCN 

MOLE 


DEFINITION 

Array  the  same  length  as  ADEBUG  common 

Array  the  same  length  as  AMPZIJ  common 

Array  the  same  length  as  ARGCOM  common 

Array  the  same  length  as  CSYSTM  common 

Array  the  same  length  as  DEFDAT  common 

Array  the  same  length  as  FLDCOM  common 

Array  the  same  length  as  GEODAT  common 

Array  the  same  length  as  GTDDAT  common 

Index  to  commons 
Save  the  index 
End-of-file  indicator 

File  number  for  input  or  output  (CHKPNT  or 
MODCHK) 

INDXWB+1 

Size  of  array  NAMOLD 
Array  the  same  length  as  INTMAT  common 

Array  the  same  length  as  IOFLES  common 

Saved  value  of  walk  back  table  index 
Array  the  same  length  as  JUNCOM  common 

Array  the  same  length  as  MODULE  common 
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RWCOMS 


(GTO,  INPUT,  MOM,  OUTPUT) 


NAMCOM  Contains  the  names  of  commons  written  to 

the  checkpoint  file 

NAME  Name  of  common  block  being  read  in 

NAMOLO  Array  to  store  first  four  locations  of 

walkback  table 

NCOMSZ  Array  containing  the  length  of  each  common 

NMWRDS  Number  of  words  in  common  being  read 

NREAD  Input  variable  that  will  indicate  a  read  or 

write 


NUMCOM 

PARTB 

SCNPR 

SGMNT 

SMSTR 

SYSFL 

TEMP 

I/O  VARIABLES: 
A.  INPUT 
ADEBG 
AMPZJ 
ARGCM 
CSTM 
OFDT 
FLDCM 
GEODT 
GTDDT 
IFILE 


Number  of  commons 

Array  the  same  length  as  PARTAB  common 
Array  the  same  length  as  SCNPAR  common 
Array  the  same  length  as  SEGMNT  common 
Array  the  same  length  as  SYMSTR  common 
Array  the  same  length  as  SYSFIL  common 
Array  the  same  length  as  TEMP01  common 


LOCATION 

/ADEBUG/ 

/AMPZIJ/ 

/ARGCOM/ 

/CSYSTM/ 

/DEFDAT/ 

/FLDCOM/ 

/GEOOAT/ 

/GFDOAT/ 

F.P. 
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(GTO,  INPUT,  MOM,  OUTPUT) 


IOFLS 

ISOFF 

ISON 

JNCN 

LUPRNT 

MDLE 

MX WALK 

NAMRTN 

NREAD 

NRSUBS 

PARTB 

SCNPR 

SGMNT 

SMSTR 

SYSFL 

TEMP 

OUTPUT 

ADEBG 

AMPZJ 

ARGCM 

CSTM 

DFDT 

FLOCM 


/IOFLES/ 

/ADEBUG/ 

/ADEBUG/ 

/JUNCOM/ 

/ADEBUG/ 

/MODULE/ 

/ADEBUG/ 

/ADEBUG/ 

F.P. 

/ADEBUG/ 

/PARTAB/ 

/SCNPAR/ 

/SEGMNT/ 

/SYMSTR/ 

/SYSFIL/ 

/TEMPOl/ 

LOCATION 

/ADEBUG/ 

/AMPZ1J/ 

/ARGCOM/ 

/CSYSTM/ 

/DEFDAT/ 

/FLDCOM/ 

/GEODAT/ 


GEODT 


RWCOMS 


(GTD,  INPUT,  MOM,  OUTPUT) 


1.  NAME:  RWFILS 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  Read  or  write  all  data  files  to/from  a  checkpoint  file. 


METHOD:  The  symbol  table  is  searched  for  a  data  file  that  is  not 
null.  The  data  set  name  is  read/written  and  file  attributes 
calculated.  PUTSYM  is  called  to  handle  the  file-to-file  transfer 
via  the  TEMP  or  SEGTBL  array.  Should  TEMP  not  be  large  enough  to 
handle  one  record  of  the  file,  a  fatal  error  is  generated. 


INTERNAL  VARIABLES: 


VARIABLE 


I  FI  LE 


KGEOM 


NBITWD 


NFILE 


NPRBUF 


NPRELM 


NPRREC 


NREAD 


NRECS 


NUMREC 


DEFINITION 


Pointer  to  symbol  table  entry 


Checkpoint  file  number,  either  IOCKPT  or 
MODCHK 


Internal  variable  equal  to  zero 


Flag  indicating  geometry  data  set 


User-assigned  name  of  data  set 


Data  set  bit  attribute  word 


Internal  variable  to  save  file  length  while 
reinitializing  file 


Logical  unit  number  for  data  set 


Number  of  records  which  will  fit  into 
buffer  array 


Number  of  words  per  data  set  element 


Number  of  words  on  data  set  record 


Flag  Indicating  whether  to  read  (ISOFF)  or 
write  (ISON)  files 


Internal  variable  equal  to  NPRREC 


Hollerith  equivalent  of  NAME 


Number  of  records  contained  in  data  set 


*  »'•  , 


RWFILS 


(GTD,  INPUT,  MOM,  OUTPUT) 


I/O  VARIABLES: 

A.  INPUT 

LOCATION 

IFILE 

F.P. 

ISOFF 

/ADEBUG/ 

ISON 

/ADEBUG/ 

KBCPLX 

/PARTAB/ 

KBGEOM 

/PARTAB/ 

KOLB IT 

/PARTAB/ 

KOLCOL 

/PARTAB/ 

KOLLOC 

/PARTAB/ 

KOLNAM 

/PARTAB/ 

KOLROW 

/PARTAB/ 

LUPRNT 

/ADEBUG/ 

NOATBL 

/PARTAB/ 

NDFILE 

/IOFLES/ 

NPDATA 

/PARTAB/ 

NREAD 

F.P. 

NTEMPS 

/TEMPOl/ 

SEGTBL 

/SEGMNT/ 

TEMP 

/TEMPOl/ 

J.  OUTPUT 

LOCATION 

IERRF 

/ADEBUG/ 
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RWFILS  (GTO,  INPUT,  MOM,  OUTPUT) 

CALLING  ROUTINES:* 

RESTRT  (1) 

STRTUP  (2,3,4) 

WRTCHK  (1,2, 3, 4) 

CALLED  ROUTINES: 

ASSIGN 

CLSFIL 

CONVRT 

IBITCK 

OPNFIL 

PUTSYM 

RDEFIL 

STATIN 

STATOT 

WLKBCK 

WRTFIL 


*  1  -  INPUT 

2  -  GTO 

3  -  MOM 

4  -  OUTPUT 


(GTD,  INPUT,  MOM,  OUTPUT) 


NAME:  SCALE2  (MOM,  OUTPUT) 

PURPOSE:  To  scale  a  linear  axis. 

METHOO:  Given  a  minimum  value  XMIN,  a  maximum  value  XMAX,  and  the 
number  of  intervals  N,  SCALE2  finds  a  new  maximum  XMAXP,  a  new  mini¬ 
mum  XMINP,  and  the  size  of  the  intervals  DIST. 

INTERNAL  VARIABLES: 


VARIABLE 

A 

AL 

B 

DEL 

DIST 

FM1 

FM2 

Ml 

M2 

N 

NAL 

NP 

VINT 

XMAX 

XMAXP 

XMIN 


DEFINITION 

Approximate  interval  size 
Log10  of  A 
A  scaled 
Round-off  error 

Distance  between  each  scale  mark 

Minimum  value  divided  by  the  interval  size 

Maximum  value  divided  by  the  interval  size 

Variable  used  to  keep  points  within  the 

minimum  value 

Variable  used  to  keep  points  within  the 

maximum  value 

Input  number  of  Intervals 

Variable  used  to  make  the  minimum/maximum 
interval  large  enough 

Check  for  the  number  of  intervals 

Array  containing  number  of  interval  sizes 

Input  maximum  value 

Output  maximum  value 

Input  minimum  value 

Output  minimum  value 


XMINP 


SCALE2  (MOM,  OUTPUT) 


5.  I/O  VARIABLES: 


INPUT 

LOCATION 

LUPRNT 

/ADEBUG/ 

N 

F.P. 

XMAX 

F.P. 

XMIN 

F.P. 

OUTPUT 

LOCATION 

DIST 

F.P. 

XMAXP 

F.P 

XMINP 

F.P. 

6.  CALLING  ROUTINE: 
PAGPLT 

7.  CALLED  ROUTINES: 
ASSIGN 

STATIN 

STATOT 

WLKBCK 
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NAME:  SCALE3 


(MOM,  OUTPUT) 


PURPOSE:  To  scale  a  log  axis. 

METHOD:  Given  a  minimum  value  XMIN,  a  maximum  value  XMAX,  and  N 

intervals  where  N  is  greater  than  one,  SCALE3  finds  a  new  range 
XMINP  and  XMAXP  divided  into  N  equal  logarithmic  intervals.  The 
ratio  of  adjacent  uniformly  spaced  scale  values  is  DIST. 

INTERNAL  VARIABLES: 


VARIABLE 

A 

AL 

B 

DEL 

DIST 

DISTL 

FM1 

FM2 

FN 

Ml 

M2 

N 

NAL 

NP 

NX 

VINT 

XMAX 


DEFINITION 

Approximate  interval  size 
Log10  of  A 
A  scaled 
Round-off  error 

Output  distance  between  scale  marks 
Approximate  interval  size 
Approximate  minimum  limits 
Approximate  maximum  limits 
Number  of  intervals 
Minimum  limit  factor 
Maximum  limit  factor 
Input  number  of  decades 
Integer  power  of  the  full  interval 
Approximate  number  of  decades 
Minimum  power 

Integer  containing  number  of  decades 
Input  maximum  value 
Logarithmic  value  of  XMAX 


XMAXL 


SCALE3 


(MOM,  OUTPUT) 


7. 


XMAXP 


ASSIGN 

STATIN 

STATOT 

WLKBCK 


Output  maximum  value 


XHIN 

Input  minimum  value 

a 

XMINL 

Logarithmic  value  of  XMIN 

XMINP 

Output  minimum  value 

- - - 

I/O  VARIABLES: 

...  ; 

A.  INPUT 

LOCATION 

*•  ;• 

LUPRNT 

/ADEBUG/ 

>'■ 

N 

F.P. 

XMAX 

F.P. 

XMIN 

F.P. 

-  -V 

B.  OUTPUT 

LOCATION 

DIST 

F.P. 

•;  >1 

XMAXP 

F.P. 

XMINP 

F.P. 

CALLING  ROUTINE: 

PAGPLT 

* - H 

CALLED  ROUTINES: 

,  ^ 

■  .-1 

SSI 
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1.  NAME:  SCAN 


(INPUT) 


2.  PURPOSE:  Scan  the  input  command  and  generate  the  scan  tables  for 

the  parsing  routine. 

3.  METHOD:  SCAN  breaks  up  the  input  command  text  into  fields.  These 

are  numeric,  alpha,  keyword,  and  operator  fields.  The  numeric  field 
is  converted  into  an  integer  or  floating  point  number  and  put  in  the 
NVAL  or  VAL  array.  The  alpha  field  is  packed  into  the  NVAL  array. 
Keyword  fields  have,  the  keyword  number  put  in  the  NVAL  array.  The 
operator  number  is  put  in  the  NVAL  array.  Then  a  field  code  for 
each  will  be  put  in  the  NCODE  array. 

When  the  SCAN  table  is  completed,  the  next  command  is  read  to  look 
for  a  possible  end  of  file  before  the  END  command. 

If  the  IGNORE  flag  is  on,  the  scan  table  is  scanned  for  monadic 

signs.  These  signs  are  added  to  the  number  that  follows  it  in  the 

SCAN  table. 

4.  INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

FIRST  This  flag  tells  SCAN  that  the  input  data 

are  the  first  of  this  record 

FRAC  Fractional  portion  of  a  floating  point 

number 

IFOUND  This  flag  tells  SCAN  that  the  exponent  was 

found  for  this  number 


INTOVR 

Flag  that  gets  set  when  a 
many  digits 

number 

has  too 

IP 

Exponent  sign 

IPROVR 

Flag  for  too  many  digits  in 

the  exponent 

IPWR 

Exponent  for  a  number 

IS 

Sign  of  numbers  in  SCAN 
compressing  tables 

table. 

used  In 

ISHIFT 

Byte  shifter 

NCOMCD 

Flag  for  comment  card 

993 


SCAN 


(INPUT) 


NOCARD 

NDXENO 

NFLDS 

NPl 

NTABSV 

NTAB1 

NUMCHR 

NUMDEC 

NUMFLD 

NUMREC 

NXTCHR 

5.  I/O  VARIABLES: 
A.  INPUT 
IBLANK 
ICOMMA 
IDIG 
IGNORE 
IMINUS 
IPER 
IPLUS 
IRIGHT 
ISLASH 
ISOFF 
ISON 


Variable  containing  the  code  for  the  END 
command 

Index  to  the  NCODES  array  for  the  keyword 
END 

Number  of  entries  in  the  SCAN  table 
Index  N  plus  one 
Saves  the  value  of  NTAB 
NTAB  plus  one 

Number  of  characters  in  an  alpha  field 

Number  of  decimals  in  a  number 

Field  counter 

Record  counter 

Next  character  in  a  field 

LOCATION 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/SCNPAR/ 

/ADEBUG/ 

/ADEBUG/ 
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SCAN 


(INPUT) 


i  ■ 
£ 


|  ISTAR 

/SCNPAR/ 

r*>* 

ISYMBL 

/SCNPAR/ 

Sv 

vV 

H 

JOIG 

/SCNPAR/ 

V  - 

i 

•  KWENO 

/PARTAB/ 

FS  KWNAME 

|%Y' 

/PARTAB/ 

K 

LETR 

/SCNPAR/ 

| 

LSTCOL 

/SCNPAR/ 

LUPRNT 

/ADEBUG/ 

LUTASK 

/AOEBUG/ 

£ 

\A 

•>s 

a \\‘ 

MXANCT 

/SCNPAR/ 

MXEXFP 

/SCNPAR/ 

| 

MXFPCT 

/SCNPAR/ 

l 

MXINCT 

/SCNPAR/ 

>: 

MXSYMB 

/SCNPAR/ 

1  ..  *• 

■  -"- 
•  m 

NBYTSZ 

/AOEBUG/ 

*'•  * 

« 

NCARD 

/SCNPAR/ 

§ 

NCOOES 

/PARTAB/ 

$ 

NCOMCH 

/SCNPAR/ 

h£i 

Ifl 

NCONCH 

/SCNPAR/ 

NDFILE 

/ IOFLES/ 

NTALPH 

/ADEBUG/ 

NTEND 

/ADEBUG/ 

i 

NTFIPT 

/ADEBUG/ 

NTINT 

/ADEBUG/ 

1 

NTKEYW 

/ADEBUG/ 

r«r.  V.v/.- v.v.* ■•/•-.v.v.v 

HifejMaiaaaisai 


*.  •*.  a 

V  '.-••• 


995 


J  J- JL  AV  >-  W-  ■■-  v.  .'■ 


i 

NTSYMB 

/ADEBUG/ 

NVALMX 

/SCNPAR/ 

B.  OUTPUT 

LOCATION 

NARGS 

/SCNPAR/ 

i  ■ 

NCARD 

/SCNPAR/ 

NCARDS 

/SCNPAR/ 

i 

i 

NCCARD 

/SCNPAR/ 

j 

NCHAR 

/SCNPAR/ 

vl 

$ 

> 

NCODE 

/SCNPAR/ 

j 

NDFILE 

/IOFLES/ 

1 

wl 

I,- 

NOGOFG 

/ADEBUG/ 

4 

m 

,  4 

NSCNER 

/SCNPAR/ 

(  4 

j 

NSCOL 

/SCNPAR/ 

0 

NTAB 

/SCNPAR/ 

% 

« 

t 

NVAL 

/SCNPAR/ 

VAL 

/SCNPAR/ 

6.  CALLING  ROUTINES: 


INPORV 

WYRDRV 

7.  CALLED  ROUTINES: 
ASSIGN 
GETKWQ 
STATIN 
STATOT 
WLKBCK 


i 


i 


NAME:  SCLRPL 


(GTD) 


PURPOSE:  To  compute  the  unobstructed  electric  field  from  a  unit 
source  scattered  by  an  elliptic  cylinder  and  then  reflected  by  plate 
MP  in  the  given  far-field  observation  direction  or  to  the  given 
near-field  observation  point.  The  cylinder  scattering  can  occur  by 
a  creeping  wave  and  a  reflected  ray  or  by  two  creeping  waves. 

METHOD:  SCLRPL  is  the  driver  routine  which  directs  all  the  ray 
tracing  physics  and  field  calculations  for  determining  the  electric 
field  scattered  by  a  cylinder  and  then  reflected  by  plate  MP  in  a 
given  far-field  direction  or  to  the  given  near-field  observation 
point.  The  ray  paths  are  shown  in  figures  1  and  2. 

The  code  begins  by  making  two  special  checks.  One  check  is  to 
determine  if  the  source  location  is  in  the  paraxial  region  beyond 
the  end  caps  such  that  it  cannot  illuminate  the  curved  surface  of 
the  cylinder.  The  other  check  is  to  make  sure  the  ray  path  which 
leaves  the  cylinder  can  cause  reflection  to  occur  from  plate  MP  in 
the  correct  direction.  If  the  result  of  either  of  these  checks  is 
such  that  cylinder-scattered  plate-reflected  fields  cannot  be  found, 
the  code  will  print  debug  information,  if  requested,  on  file  LUPRNT. 
Control  is  then  returned  to  the  calling  routine.  If  the  correct 
interaction  can  be  found,  the  tangent  vectors  from  the  source  to  the 
cylinder  are  specified.  Ray  tracing  and  field  computations  are  per¬ 
formed  for  both  of  these  tangents.  The  type  of  field  calculation  is 
based  on  the  value  ALR.  ALR  is  the  reflected  ray  phi  angle  in  the 
tangent  coordinate  system.  If  it  is  approximately  it,  grazing  inci¬ 
dence  occurs.  If  ALR  is  greater  than  n,  a  cylinder  creeping  wave 
occurs  for  this  tangent.  If  ALR  is  less  than  tt,  a  reflected  wave 
occurs.  After  performing  the  ray  path  calculations  and  field  calcu¬ 
lations,  the  code  will  check  to  see  if  the  second  tangent  remains  to 
be  addressed.  If  it  does,  a  new  value  of  ALR  will  be  computed  and 
the  appropriate  field  calculations  will  be  performed.  If  both 
tangents  have  been  addressed,  debug  information  (if  requested)  will 
be  printed  on  file  LUPRNT.  The  debug  information  consists  of  the 
total  field  magnitude  and  the  theta  and  phi  components  of  the  field. 
Control  Is  then  returned  to  the  calling  routine. 

If  ALR  Is  less  than  it,  it  is  possible  that  a  reflected  ray  can  occur 
from  the  cylinder.  Subroutine  RCLRPL  calculates  the  ray  path  for  a 
ray  reflected  by  a  cylinder  and  then  reflected  by  a  plate.  It  also 
computes  the  incident  field  upon  the  cylinder  at  the  cylinder 
reflection  point  and  other  geometry-specific  terms.  These  terms  and 
field  are  passed  to  SCLRPL  through  common  block  FUDGJ.  The  ray 
spreading  radii  and  phase  factor  are  computed  in  this  subroutine 
based  upon  parameters  from  common  block  FUDGJ.  Then  the  field 
computations  can  begin.  First,  the  cylinder-reflected  field  Is 
computed.  Then  the  field  reflected  from  plate  MP  is  computed  in 
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Figure  1.  Illustration  of  Ray  Scattered  by  the  Cylinder  and 
Reflected  by  a  Plate  into  the  Far  Field 
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Figure  2.  Illustration  of  Ray  Scattered  by  Cylinder  and 

Reflected  by  Plate  to  the  Near-Field  Observation 
Point 
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theta  and  phi  components.  This  is  the  total  cylinder-relected 
plate-reflected  field.  Subroutine  XYZFLO  is  called  to  compute  the 
x,  y,  z  components  of  the  total  field  and  to  accumulate  them  with 
other  fields  from  Drevious  interactions. 

If  ALR  is  approximately  equal  to  n,  grazing  incidence  occurs.  The 
grazing  incidence  portion  of  the  code  first  checks  to  make  su^e  that 
reflection  from  plate  MP  can  occur.  If  it  cannot,  the  cooe  will 
check  to  see  if  the  second  tangent  remains  and  will  address  it 
appropriately.  If  reflection  from  plate  MP  can  occur,  the  ray  path 
is  checked  to  see  if  it  is  shadowed  anywhere.  If  it  is,  the  code 
will  check  for  another  tangent.  If  shadowing  did  not  occur  the 
grazing  incidence  point  is  computed.  This  point  is  checked  to  make 
sure  that  it  is  on  the  curved  surface  of  the  cylinder.  The  source 
field  pattern  factor  is  found  by  calling  subroutine  SOURCE.  Next, 
the  phase  factor  is  computed.  The  hard  and  soft  theta  and  phi  com¬ 
ponents  of  the  cylinder  field  are  now  found.  The  grazing  incidence- 
scattered  field  is  now  computed.  From  this  the  total  field  being 
made  up  of  the  grazing  incidence-scattered  field  and  the  field 
reflected  from  the  plate  is  computed  in  theta  and  phi  components. 
The  x,  y,  z  components  are  computed  by  calling  subroutine  XYZFLO. 

If  ALR  is  greater  than  it,  a  creeping  wave  can  occur  from  the 
cylinder.  The  location  of  the  point  at  which  the  creeping  wave 
leaves  the  cylinder  and  the  point  at  which  the  creeping  wave  begins 
on  the  cylinder  are  computed  in  two  different  manners  depending  upon 
whether  near  field  or  far  field  was  requested.  For  the  computation 
sequence,  see  the  accompanying  flowchart.  While  in  the  separate 
field  computation  paths,  the  code  checks  the  ray  from  the  source  to 
the  cylinder  for  obstructions.  If  the  path  is  blocked,  the  code 
proceeds  to  the  next  tangent.  If  the  path  is  clear,  calculations 
continue  for  this  tangent.  The  incident  source  field  pattern  factor 
is  computed  by  calling  subroutine  SOURCE.  For  both  near  field  and 
far  field,  the  code  checks  to  see  if  reflection  can  occur  from  plate 
MP.  If  it  cannot,  the  code  will  look  to  see  if  a  tangent  remains  to 
be  addressed.  If  reflection  does  occur  on  plate  MP,  the  ray  from 
the  cylinder  to  the  plate  and  into  the  observation  direction  or  to 
the  observation  point  is  checked  for  any  obstruction.  If  it  is 
obstructed,  the  code  will  proceed  to  the  next  tangent.  If  the  ray 
in  unobstructed,  the  phase  factor  is  computed,  the  hard  and  soft 
creeping  waves  are  computed,  and  then  the  total  cylinder  creeping 
wave  field  is  computed.  From  this,  the  total  field  scattered  by  the 
cylinder  and  reflected  by  the  plate  is  computed  in  theta  and  phi 
components.  Subroutine  XYZFLO  is  called  to  compute  and  accumulate 
the  x,  y,  z  components. 
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INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

A1,A2 

Field  components  of  ray  incident  on  plate 
normal  and  tangent  to  plate 

A3 

Determinant  of  polarization  transformation 

ALR 

Phi  angle  defining  propagation  direction  in 
tangent  point  coordinate  system  (2-D) 

ALRS 

Difference  between  ALS  and  ALR 

ALS 

Phi  angle  defining  direction  of  ray  from 
RCS  origin  to  source  in  tangent  point  coor¬ 
dinate  system 

AN 

Distance  from  plate  MP  plane  to  the  obser¬ 
vation  field  point  (from  subroutine  IMAGE) 

ANR 

Distance  from  plate  plane  to  XRF,  the  point 
at  which  the  creeping  wave  leaves  the 
cylinder  (from  subroutine  IMAGE) 

BX.BY.BZ 

X,Y,Z  components  of  polarization  unit 
vector  of  soft  component  field  incident  on 
cylinder  (parallel  to  cylinder  surface  and 
normal  to  incident  field  propagation 
direction) 

C11,C12,C2I,C22 

Coefficients  used  to  convert  polarization 
from  theta  and  phi  components  in  RCS  to 
components  normal  and  tangent  to  plate  (and 
vice-versa) 

CCC 

Real  part  of  Fresnel  integral  (from  subrou¬ 
tine  FRNELS) 

CF 

Phase  term  and  ray  spreading  factor 

CFH 

Hard  transition  field  coefficient 

CFS 

Soft  transition  field  coefficient 

DEPH 

Phi  component  of  transition  field  in  RCS 

DETH 

Theta  component  of  transitior  field  in  RCS 

DHIT 

Distance  from  source  to  plate  hit  point 
(from  subroutine  PLAINT) 

DHJT 

Distance  from  where  ray  leaves  cylinder  and 
hits  plate  MP  (from  subroutine  PLAINT) 

DHT 

Distance  to  hit  point  (from  subroutine 
PLAINT  and  CYLINT) 

DICOEF 

X,Y,  and  Z  components  of  incident  ray 
direction  on  cylinder  in  RCS 

OIJ 

X-Y  plane  vector  from  a  source  tangent  to 
the  point  the  creeping  wave  leaves  the 
cylinder 

DIJXDJ 

Cross  product  of  DIJ  and  DJT 

DIT 

Incident  ray  vector 

DIXDIJ 

Cross  product  of  DIT  and  DIJ 

OJ 

X,Y,Z  components  of  propagation  direction 
of  ray  between  cylinder  and  plate  in  RCS 

DJI 

Direction  unit  vector  towards  XE1 

DJ2 

Direction  unit  vector  towards  XE2 

DJT 

X-Y  plane  components  of  observation 
direction 

DMAG 

Distance  between  the  observation  field 
point  and  the  plate  reflection  point 

EF 

Theta  component  of  source  field  pattern 
factor  in  RCS 

EG 

Phi  component  of  source  field  pattern 
factor  in  RCS 

EHP 

Phi  component  of  hard  component  of  geomet¬ 
rical  optics  field  incident  on  cylinder  in 
RCS 

EHT 

Theta  component  of  hard  component  of  geo¬ 
metrical  optics  field  incident  on  cylinder 

SCLRPL  (GTO) 


X  component 
factor 

of 

the 

source 

field 

pattern 

Y  component 
factor 

of 

the 

source 

field 

pattern 

Z  component 
factor 

of 

the 

source 

field 

pattern 

Phi  component  of  scattered-ref lected  field 
in  RCS 

Oot  product  of  cylinder  tangent  unit  vector 
and  reflected  ray  propagation  direction 
(2-D) 

Phi  component  of  geometrical  optics 
ref lected-ref lected  field  in  RCS 

Theta  component  of  geometrical  optics 
ref lected-ref lected  field  in  RCS 

Phi  component  of  soft  component  of  geomet¬ 
rical  optics  field  incident  on  cylinder  in 
RCS 

Theta  component  of  soft  component  of  geo¬ 
metrical  optics  field  incident  on  cylinder 
in  RCS 

Theta  component  of  scattered-ref lected 
field  in  RCS 

Near-field  observation  field  point  image 
location  (imaged  through  plate  MP) 

The  x-y  plane  location  of  the  field  point 
image 

Intermediate  variable  for  transition 
function 

Logical  variable  set  true  if  a  plate  is  hit 
(from  subroutine  PLAINT  and  CYLINT) 


Logical  variable  set  true  if  cylinder- 
reflected  field  is  not  present 


SCLRPL 


(GTD) 


LVJ 

MP 

ORIGIN 

PHIR 

PHJR 

PHJR1 

PHJR2 

RGF 

RGI 

RT 

S 

51 

52 

SKWIG 

SNF 

SPHJ 


Logical  variable  set  true  the  first  time 
creeping  wave  computations  begin 

Plate  where  reflection  occurs 

Origin  of  the  reference  coordinate  system 

(0.,  0.,  0.) 

Phi  angle  of  incident  ray  direction  on 
cylinder 

Phi  angle  of  ray  propagation  direction 
between  cylinder  and  plate 

Phi  angle  of  DJI 

Phi  arsle  of  DJ2 

Longitudinal  radius  of  curvature  of  the 
cylinder  at  the  point  the  creeping  wave 
leaves  the  cylinder 

Longitudinal  radius  of  curvature  of  the 
cylinder  at  the  creeping  wave  incidence 
point 

Transverse  radius  of  curvature  of  the 
cylinder 

The  total  distance  between  the  source  point 
and  the  incident  point 

X-Y  plane  distance  between  the  source  and 
incident  point 

X-Y  plane  distance  between  the  field  point 
image  and  the  point  on  the  cylinder  at 
which  the  creeping  wave  leaves  the  cylinder 

Parameter  used  in  transition  function 

Distance  between  plate  reflection  point  and 
field  point 

Sine  of  PHJR 


SS 


Distance  of  path  along  the  cylinder 
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NAME:  SOURCP 


(GTO) 


PURPOSE:  To  compute  the  tangential  components  of  the  normal  deriva¬ 
tive,  of  the  Incident  field  pattern  factor  for  a  source  ray 
incident  on  a  given  edge. 

METHOD:  A  source  Is  located  and  oriented  according  to  figure  1  and 
emits  a  ray  Incident  on  a  diffracting  edge.  The  slope  field  Is 
given  by 


3n  '  s'sin  8„  3$^ 
o  o 


where 


E  =  Eo(e'^') 


-jks' 


0'  and  are  the  spherical  angles  of  the  source  ray  In  the  source 
coordinate  system,  s'  Is  the  distance  from  the  source  to  the  edge. 


DIFFRACTING 

EDGE 


EDGE  CENTERED 
SYSTEM 


SOURCE 

COORDINATE 

SYSTEM 


SOURCE 

CENTER 


Figure  1.  Pertlnenfc  Gegpietiry  for  Slope  Diffraction.  The  Unit 
Vectors  0',  8  and  $  Are  Defined  In  the 

Standard  Hay.  0  0 


PREVIOUS  PAGE 
IS  BLANK 


SOURCP 


(GTO) 


The  required  derivative  In  (1)  Is  computed  from  the  spherical 
components  of  the  source  field  and  their  derivatives  (In  the  source 
coordinate  system)  by  using  the  chain  rule  of  taking  derivatives: 


ai  3<E0f0’)  3<V*’> 

n  i  3  r\  i  k  l 


94) 


3$ 


H. 


(3) 


«  i  <\8i  8  |  jai 

=  e’  +  Ee*  +  ♦'  +  V  4^ 


ET 


a* 


3E«, 


'0' 


iil  4.  3E9'  Ml 

80'  80q  8$ '  30q 


ETT 


ETP 


EP 


(4) 

(5) 


8E^, 

3*o 


3V  801  ,  aV  Ml 
80'  8*c  +  80 '  8*q 

EPT  EPP 


(6) 


The  angular  derivatives  are  given  by 


801 


-sin  b; 


0' 


Hi 


sin  B '  ~  ~ 

sin  0'  vo  v 


TPHO 


PPHO 


H;  -  *in  e; 


<♦, 


0’)  s'  -  cot0 '  (A  •  ♦’)  ♦' 


(8) 


TPHO 


VI 


PPHO 


I 


I, 


sin  p  a  a  a 

SIST*  <♦„  •  «' 


PPHO 


p*  «  sin  0'  s'  +  cos  0*  0’ 


STHP  VI 


0*  =  (XTH,  YTH,  ZTH) 

« 

♦'  =  (XPH,  YPH,  ZPH) 


CTHP 


The  above  expressions  may  be  combined  to  yield: 


A  A 

oil  A  1 


••)  «•  -  sr-  <♦„  •  9’>  ♦ 


JB  1  r  3Ea|  a  A  A  3E  .  |  A  A  A  1 

IS  *  <♦.  •  9'>  91  •  **-<*»■  9,>  *'J 

1  r  3Ea  |  A  A  A  3E  .  ,  A  A  A  "I 

♦  [-  tt~  «♦.  ■  ♦’>  9'  -  w*-  «♦.  •  *'J 

[<*  A  A  A  A  "  I  (  a"!* 

-  Ee*  (♦o  *  *’}  ♦'  +  V  <*o  •  ♦’>  e'J  (  “7T2 


Note  that  sin  Bq  Is  eliminated  from  the  expression  In  (13).  Also, 
only  the  tangential  components  are  retained. 

e-jks' 

SOURCP  returns  the  quantity  within  the  brackets,  as  the  term  — — 

s'* 

Is  added  elsewhere  In  the  code.  The  two  field  components  EIPRP  and 
EIPLP  are  given  by: 


ft#  .-jks' 

g  -  [EIPRP  *  +  EIPLP  B' J  - — 5 — 
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EIPRP  Is  computed  by  taking  the  dot  product  of  (13)  with  L, 
Ignoring  the  exponential  term.  EIPLP  Is  obtained  by  dotting  (13) 
with  §o* 

Since  Ea  1  and  E^'  for  GTD  and  HOH  sources  are  analytic  In  the  far 
field,  it  Is  straightforward  to  obtain  the  necessary  partial  deriva¬ 
tives.  For  example,  consider  a  wire  segment  with  pulse  current 


ET  -  E01 


In 

2  'X  ’ 


sin  0 ' 


=  0 


ETT 


ETP 


3E0, 

30' 

3Ee. 

34' 


(|^)  cos  0'  EPT 
0  EPP 


(15) 


These  terms  are  then  substituted  Into  the  formulas  for  EIPRP  and 
EIPLP  at  statement  100. 


4.  INTERNAL  VARIABLES: 
VARIABLE 
ACTHP 

BOP 

CCDK2 

CDK2 

CNDK2 

CONST 

CPHP 

crop 


DEFINITION 

Absolute  value  of  the  cosine  of  0'  (source 
coordinate  system) 

Rectangular  components  of  the  beta- 
component  of  the  Incident  field  In  RCS 
system 

Cosine  of  CDK2 

Wave  number  times  wire  half-length  times 
cosine  of  0' 

Cosine  of  0*2 

jn(%AA)*X  used  In  computing  patch  slope 
fields 

Cosine  of  4 '(source  coordinate  system) 

Cosine  of  0'  (source  coordinate  system) 
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DK2 

E1,E2 

EA.EB 

EFA 

EFB 

EIPLP 

EIPRP 

EP 

EPP 

EPT 

ET 

ETP 

ETT 

IM 

PHO 


Wave  number  times  wire  half-length 

Temporary  storage  for  terms  used  to  compute 
slope  fields  of  stiff  dipole  source 

Temporary  storage  for  terms  used  to  compute 
slope  fields  of  stiff  dipole  source 

Partial  derivative  of  dipole  pattern  factor 
with  respect  to  theta 

Dipole  pattern  factor  divided  by  sin(e') 

Parallel  polarized  component  (soft)  of 
normal  derivative  of  Incident  fields 
(parallel  to  edge) 

Perpendicularly  polarized  (hard)  component 
of  normal  derivative  of  incident  fields 
(perpendicular  to  edge) 

Phi  polarized  component  of  incident  field 
(source  coordinate  system) 

Partial  derivative  of  EP  with  respect  to 

phi 

Partial  derivative  of  EP  with  respect  to 

theta 

Theta  polarized  component  of  incident  field 
(source  coordinate  system) 

Partial  derivative  of  ET  with  respect  to 

phi 

Partial  derivative  of  ET  with  respect  to 

theta 

GTD  source  type 

Rectangular  components  of  the  phi -component 
of  the  Incident  field  In  RCS 

Dot  product  of  phi  polarization  unit  vector 
of  source  coordinate  system  and  beta  polar¬ 
ization  unit  vector  of  edge-centered  coor¬ 
dinate  system 


PPBO 
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PPHO 


RDX 

ROY 

SCDK2 

SN 

SNDK2 

SP1 

SP2 

SPHP 

STHP 

TPBO 


TPHO 


VAX 

VI 


Oot  product  of  phi  polarization  unit  vector 
of  source  coordinate  system  and  phi  polari¬ 
zation  unit  vector  of  edge-centered  coor¬ 
dinate  system 

Projection  of  incident  ray  direction  onto 
source  x  axis 

Projection  of  incident  ray  direction  onto 
source  y  axis 

Sine  of  CDK2 

Sign  of  CTHP 

Sine  of  0K2 

Source  parameter  one:  wire  radius  (wave¬ 
lengths)  or  patch  area  (square  wavelengths) 

Source  parameter  two:  wire  length  (wave¬ 
lengths) 

Sine  of  (source  coordinate  system) 

Sine  of  0'  (source  coordinate  system) 

Oot  product  of  theta  polarization  unit  vec¬ 
tor  of  source  coordinate  system  and  the 
beta  polarization  unit  vector  of  edge-cen¬ 
tered  coordinate  system 

Dot  product  of  theta  polarization  unit  vec¬ 
tor  of  source  coordinate  system  and  the  phi 
polarization  unit  vector  of  edge-centered 
coordinate  system 

Source  axes  direction  cosines  (RCS  rectan¬ 
gular  component  projections) 

Direction  cosines  of  Incident  ray  propaga¬ 
tion  direction 


XPH,YPH,2PH  Rectangular  components  of  the  phi  unit 

polarization  unit  vector  In  the  source 
coordinate  system  (RCS  components) 


XTH.YTH.ZTH 


Rectangular  components  of  the  theta  unit 
polarization  unit  vector  In  the  source 
coordinate  system  (RCS  components) 


I/O  VARIABLES: 
A.  INPUT 


LOCATION 


/COMP/ 

/AMPZIJ/ 

/SRC/ 


/PIS/ 

/SRC/ 

/SRC/ 

/PIS/ 


B.  OUTPUT 
EIPLP 
EIPRP 

CALLING  ROUTINES: 

DIFPLT 

DPLRPL 

RPLDPL 

CALLED  ROUTINES: 

ASSIGN 

STATIN 

STATOT 

rfLKBCK 


LOCATION 


SOUR CP 


\ / 


CALCULATE  COMPO¬ 
NENTS  OF  PROPAGA¬ 
TION  DIRECTION 
ALONG  SOURCE  AXES 


CALCULATE  POLARIZA¬ 
TION  UNIT  VECTORS 
FOR  RAY  IN  SOURCE 
COORDINATE  SYSTEM 


PROJECT  POLARIZA¬ 
TION  VECTORS 
PARALLEL  &  PERPEN¬ 
DICULAR  TO  PLATE 
EDGE 

~r~ 


INITIALIZE  FIELDS 
TO  ZERO 


I - 

COMPUTE  INCIDENT 
FIELDS  AND  DERIVA¬ 
TIVES  FOR  SPECIFIED 
SOURCE  TYPE 


100 


COMPUTE  THE  PAR¬ 
ALLEL  ANO  PERPEN¬ 
DICULAR  SLOPE  FIELD 
COMPONENTS 


NAME:  SPWDRV 


(MOM) 


PURPOSE:  Generate  the  plane  or  spherical  wave  excitation  on  the 
structure. 

METHOD:  The  coordinate  system  and  parameters  are  illustrated  In 
figure  1,  where  E  is  the  linear  component  and  EP  is  the  polarization 
component  of  the  incident  electric  field.  The  eccentricity  e 
specifies  |fP|/_|E|_and  for  positive  e,  the  wave  vector  k  has  direc¬ 
tion  given  by_E  x  EP  for  negative  e,  the  wave  vector  k  has  direction 
glvenby  E  x  EP  and  for  negative  £»  _  k  has  direction  given  by 
EP  x  E.  For  plane  wave  excitation,  k  is  oriented  toward  the  origin; 


Figure  1.  Geometry  for  Plane  or  Spherical  Wave  Excitation 
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whereas  for  spherical  wave  excitation,  k  is  oriented  toward  the 
observation  point.  The  total  field  is  given  by 


Efc  =  E1  +  ER 


E1  =  E  +  jEP 

-R  -R 

where  E  is  the  reflected  field.  E  is  given  by: 


ER  =  (E1?. 


+  ER  )  R, 


-Ei  R  i 


-n  -P 

where  E  u  and  E  y  are  the  components  of  the  total  field  incident  normal  to 
(n)  and  parallel  to  (P)  the  reflective  surface.  Elis  the  total  incident 
field  perpendicular  to  the  plane  of  reflection.  R  ||  and  Ri  are  the 
modified  Fresnel  reflection  coefficients  for  the  in-plane  and  out-of-plane 
components  as  described  in  the  Engineering  Manual. 


The  excitation  for  a  segment  located  at  R.  oriented  parallel  to 
with  length  I |  is 


E.  =  -l.  •  (R. ) 
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INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

ARGI 

ARGR 

COSARG 

—  k  •  R 

X  Ki  Ki 

in  Z  .  D 

X  1  l 

-jk  •  R. 

Real  component  of  e 

COSETA 

(E  •  $)/IEI 

COSP 

Cosine  of  $ 

COST 

Cosine  of  0 

DXSN 

X  component  of  source  with  respect  to 
specular  point 

DYSVI 

Y  component  of  source  with  respect  to 
specular  point 

DZSW 

Z  component  of  source  with  respect  to 
specular  point 

ECCN  «  ECCEN 

Eccentricity  (IEPI/IEI) 

EN 

Magnitude  IEI 

EPRX.EPRY.EPRZ 

X,Y,Z  components  of  reflected  polarization 
component  of  incident  wave 

EPX.EPY.EPZ 

X.Y,Z  components  of  reflected  linear 
component  of  Incident  wave 

ERX,  ERY,  ERZ 

XtY,Z  components  of  reflected  linear 
component  of  incident  wave 

ESX,ESY,ESZ 

X,Y,Z  components  of  linear  component  of 
Incident  wave 

ETAE 

Angle  between  E  and  0 

ETAINV 

1/377  (mho) 

ETAP 

Polarization  angle 

EXI 

Total  x  component  of  Incident  wave 
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EXPARG 

EXR 

EXRI 

EXS 

EYI 

EYR 

EYR1 

EYS 

EZI 

EZR 

EZRI 

F 

GROUND 

HCNVRT 

HXI.HYI.HZI 

HXR.HYR.HZR 

KIX 

KIXSQ 

KIY 

KIYSQ 

KIZ 

KRX 

KRXSQ 


-jk  •  R, 
e 

Total  x  component  of  reflected  wave 

In-plane  reflected  x  component 

Linear  x  component 

Total  y  component  of  Incident  wave 

Total  y  component  of  reflected  wave 

In-plane  reflected  y  component 

Linear  y  component 

Total  z  component  of  Incident  wave 

Total  z  component  of  reflected  wave 

In-plane  reflected  z  component 

Logical  .FALSE. 

Logical  .TRUE.  If  ground  present 

Logical  .TRUE.  If  convert  to  H-fleld  for 
plane  wave 

X,Y,Z  components  of  Incident  H-fleld  for 
patches 

XtY,Z  components  of  reflected  H-fleld  for 
patches 

kt  •  X 

(KIX) 2 

*”  A 

k1  *  Y 

(KIY) 2 

kr  2 
k]  •  X 
(KRX)2 
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KRY 

ii  .  9 

r 

KRYSQ 

(KRY)2 

KRZ 

‘V* 

KSYMP 

1  -  no  ground 

2  *  ground 

KXKY 

KIX  *  KIY 

LBLK 

Block  number  containing  subsection  being  ' 

considered 

LINEAR 

.TRUE,  for  ECC  -  0  ! 

.FALSE,  for  ECC  >  0 

NAMEXC 

Symbolic  name  of  excitation  data  set 

NDXBLK 

Block  number  of  current  geometry  data  | 

NI 

Index  to  TEMP  array  for  Imaginary  component 

NR 

Index  to  TEMP  array  for  real  component 

NUMYRS 

Number  of  wire  segments  j 

NX.NY.NZ 

Components  of  patch  normal  vector 

PHI 

Spherical  angle  ♦  In  radians 

PHIS 

Spherical  angle  $  In  degrees  j 

PLNMAV 

.TRUE,  for  plane  wave  excitation 

RF 

RFI 

RHO 

RHOSQ 

RI 


.FALSE,  for  spherical  wave  excitation 

>  0  -  location  of  spherical  wave  source 
<  0  =  plane  wave  source 

Distance  from  source  to  specular  point 

I/RF 

Reflection  plane  component  of  RF 
(RHO)2 

Distance  from  source  to  field  point 
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In-plane  reflection  coefficient 

Out-of-plane  reflection  coefficient 

Location  of  wave  excitation  source 

R2  .  . 

-jk  •  R, 

Imaginary  component  of  e 
?  1/2 

(1-COSETA  ; 

Sine  (♦) 

Sine  (0) 

Logical  .TRUE. 

Internal  symbols  for  SEGTBL  parameters 

Internal  symbols  for  SEGTBL  parameters 

8  (radians) 

8  (degrees) 


Unused 
IE  +  je  EPI 
Total  excitation 
Unused 

Logical  .TRUE,  if  element  Is  a  wire 
X  coordinate  of  field  point 
X  coordinate  of  specular  point 
X  coordinate  of  source  point 
X  component  of  £ ^ 

Y  coordinate  of  field  point 

Y  coordinate  of  specular  point 

Y  coordinate  of  source  point 


RHH 
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^4 


^  . 


YU 

ZC 

ZR 

ZRSQRT 

ZS 

ZW 

5.  I/O  VARIABLES: 

A.  INPUT 
OGTORO 
ECCEN 
ETA 
ETAP 
IPERF 
ISGTBL 
ISOFF 
ISON 
KSYMP 
MAXSEG 
NAMEXC 
NDXBLK 
NPATCH 
NWIRE 
PHIS 
RS 


Y  component  of  4^ 

Z  coordinate  of  field  point 

Z  coordinate  of  specular  point 

Intermediate  value  In  calculation  of  RINP 
and  ROUT 

Z  coordinate  of  source  point 
Z  component  of  4j 

LOCATION 

/GEOOAT/ 

F.P. 

/AMPZIJ/ 

F.P. 

/AMPZIJ/ 

/SEGMNT/ 

/ADEBUG/ 

/AOEBUG/ 

/AMPZIJ/ 

/SEGMNT/ 

F.P. 

/SEGMNT/ 

/SEGMNT/ 

/SEGMNT/ 

F.P. 

F.P. 
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SEGTBL 

THETS 


/SEGMNT/ 


WAVNUM 

ZRATI 

B.  OUTPUT 


UPDBLK 

CALLING  ROUTINE: 
EXCDRV 

CALLED  ROUTINES: 

ASSIGN 

GETSEG 

STATIN 

STATOT 

WLKBCK 


/AMPZIJ/ 

/AMPZIJ/ 

LOCATION 

/TEMPOl/ 

/SEGMNT/ 


1.  NAME:  STATFN 


(GTO,  INPUT,  MOM,  OUTPUT) 


>■'551 


PURPOSE:  Subroutine  to  print  the  timing  statistics  compiled  during 
code  execution  and  to  write  the  end-of-module  checkpoint. 

METHOD:  The  number  of  times  a  subroutine  Is  entered,  the  total  time 
to  run  the  subroutine,  and  Its  percentage  of  the  total  code 
execution  time  are  compiled  and  printed. 


INTERNAL  VARIABLES: 

VARIABLE 

IMDCHK 


ITEMS 


MODCHK 

NITEMS 


RITEMS 

TOTAL 


I/O  VARIABLES: 
A.  INPUT 
CHKPNT 
IERRF 
IOCKPT 
ISOFF 


DEFINITION 

Flag  Indicating  that  an  end-of-module 
checkpoint  Is  being  written 

Dummy  array  used  to  store  the  subroutine 
times  for  sorting 

Pointer  to  the  next  subroutine  statistics 
to  be  printed 

An  order  array  used  to  Indicate  the  sorted 
subroutine  timing  statistics 

End-of-module  checkpoint  file  logical  unit 

The  number  of  entries  for  which  there  will 
be  timing  statistics 

The  percentage  of  time  spent  In  any  one 
subroutine 

Real  array  equlvalenced  to  ITEMS 

Total  amount  of  GEMACS  code  computer  time 
accounted  for  In  subroutines 


LOCATION 

/SYSFIL/ 

/ADEBUG/ 

/SYSFIL/ 

/ADEBUG/ 


*  “  (L*  -  «  "  •  -  •  •  %  •  «_  •.  V  %  \  • 
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STATFN 


(GTD,  INPUT,  MOM,  OUTPUT) 


ISON 
LUPRNT 
MODCHK 
HODNAM 
NRNAHS 
NRSUBS 
NRTIHS 
RSUHS 
B.  OUTPUT 
COHPLT 
IHOCHK 
LSTMOO 
MODLST 
RSTRTA 

6.  CALLING  ROUTINES: 
MAIN  PROGRAM 
ERROR 

7.  CALLED  ROUTINES: 
CLSFIL 

OPNFIL 

SHELL 

WRTCHK 


/ADEBUG/ 

/AOEBUG/ 

/SYSFIL/ 

/MODULE/ 

/AOEBUG/ 

/AOEBUG/ 

/AOEBUG/ 

/AOEBUG/ 

LOCATION 

/SYSFIL/ 

/AOEBUG/ 

/MODULE/ 

/MODULE/ 

/SYSFIL/ 
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(GTD,  INPUT,  MOM,  OUTPUT) 


SET  NITEMS  TO  NUMBER 
OF  SUBROUTINES. 

SET  TOTAL  TO  0. 

SET  I  TO  1. 


SET  LOCIII  •  I 
SET  RITEMS  III  TO 
-RSUMS  (II  FOR 
DESCENDING  SORT 


NRTIMSIII  \ 
I  THAN  OR  EQUAL 
TO  ZERO?  / 


TOTAL  *  TOTAL  ♦ 
RSUMSdl 


EQUAL  TO 
WTEME 

V  ?  > 


EXECUTE  AN 

ENOQE-MOOULE 

CHECKPOINT 


NAME:  STATIN  (GTO,  INPUT,  MOM,  OUTPUT) 

PURPOSE:  To  Initialize  timing  statistics  for  all  subroutines  which 
call  It. 

METHOD:  The  current  wall  clock  time  Is  entered  Into  the  RTINS 

array.  The  NRTIMS  array  Is  Incremented  by  1  and  the  total  time  In 
the  previous  subroutine  Is  loaded  In  the  RSUMS  array. 

INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

MSAVE  Input  array  containing  the  statistics  for 

the  previous  subroutine. 

N  Subroutine  number  of  the  previous  subrou¬ 

tine  to  call  this  routine. 

NAME  Input  argument,  current  subroutine  name. 

NRTIMS  Array  to  accumulate  the  number  of  times  a 

given  subroutine  Is  entered. 

NUMSB  Subroutine  number  of  calling  subroutine. 

RSUMS  Total  time  accumulated  In  the  previous  sub¬ 

routine  to  call  this  routine. 

RTINS  The  current  clock  time  for  the  current  sub¬ 

routine  calling  this  subroutine. 

TIMIN  Current  clock  time. 

I/O  VARIABLES: 

A.  INPUT 
ISON 
LTRACE 
LUPRNT 
MSAVE 
NAME 
NUMSB 


STATIN 


(GTO,  INPUT,  MOM,  OUTPUT) 


B.  OUTPUT 
NRTIMS 
RSUMS 
RTINS 


LOCATION 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 


CALLING  ROUTINES: 
All  major  routines. 
CALLED  ROUTINES: 
NONE 
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.v 


V.Y 


.v .  ■ 
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STATIN 


(GTD,  INPUT,  MOM,  OUTPUT) 


SUBTRACT  FROM  THIS 
TIME  THE  TIME  SUB¬ 
ROUTINE  N  WAS 
ENTERED,  AND  ADD 
THAT  TO  ITS  TOTAL 
TIME _  _ 


RETURN 


NAME:  STATOT 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  To  close  the  timing  statistic  upon  exit  from  a  subroutine. 

METHOD:  The  Index  of  the  calling  subroutine  Is  retrieved.  The  wall 
clock  time  is  then  determined,  and  the  accumulated  time  for  the  cur¬ 
rent  subroutine  Is  determined.  Then  the  clock  is  restarted  for  the 
subroutine  which  called  the  calling  subroutine. 

INTERNAL  VARIABLES: 


VARIABLE 

MSAVE 

N 

NAME 

NUMSB 

RSUMS 

RTINS 

TIMOUT 


DEFINITION 

Input  argument  array  containing  the  sub¬ 
routine  which  called  the  calling  subroutine 

The  subroutine  number  of  the  subroutine 
which  called  the  calling  subroutine 

Input  argument  containing  the  name  of  the 
subroutine  for  which  the  statistic  is  being 
accumulated 

The  subroutine  number  of  the  subroutine  for 
which  the  statistic  is  being  accumulated 

The  accumulated  time  spent  in  the  subrou 
tine  for  which  the  statistic  is  being 
accumulated 

The  current  time  reset  for  the  subroutine 
which  called  the  calling  subroutine 

Current  wall  clock  time 


I/O  VARIABLES: 


INPUT 

LOCATION 

ISON 

/ADEBUG/ 

LTRACE 

/ADEBUG/ 

LUPRNT 

/ADEBUG/ 

MSAVE 

F.P. 

NAME 

F.P. 

NUMSB 

F.P. 

PREVIOUS  PAGE 
IS  BLANK 
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STATOT 


(GTD,  INPUT,  MOM,  OUTPUT) 


RSUMS 

RTINS 

6.  CALLING  ROUTINES: 
All  major  routines. 

7.  CALLED  ROUTINES: 
None. 


LOCATION 

/ADEBUG/ 

/ADEBUG/ 


n 


STATOT 


(GTD,  INPUT,  MOM,  OUTPUT) 
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NAME:  STRTUP 


(GTO,  MOM,  OUTPUT) 


PURPOSE:  Inltallze  commons  and  reset  data  files  to  begin  module 
execution. 

METHOD:  The  last  checkpoint  on  module  checkpoint  file  MODCHK  Is 
read  to  Initialize  commons  and  read  In  data  files.  If  the  RSTART 
flag  Is  on,  the  module  name  requested  In  the  RESTRT  command  Is 
compared  to  the  name  of  this  module.  If  they  do  not  match,  execu¬ 
tion  stops  so  that  this  same  checkpoint  can  be  used  by  a  subsequent 
module  for  restarting.  If  a  match  occurs,  the  routine  returns  to 
GEMACS,  which  Invokes  TSKXQT  to  continue  execution. 

For  standard  module  start-up  (no  RSTART  flag),  the  error  flag  Is 
first  checked.  If  an  error  occurred  In  a  previous  module,  execution 
Is  terminated  with  an  appropriate  warning  message.  Otherwise, 
parameters  are  reset  to  their  default  values  and  data  sets  rewound 
to  just  before  the  first  word  of  their  first  editions. 


INTERNAL  VARIABLES: 
VARIABLE 


ICHKPT 

ICKLOP 


IOSTOR 

IRSAV 


KCODE 


LOCNOW 

LUFILE 


MODNOW 

NREAD 


DESCRIPTION 

Loop  Index  pointing  to  symbol  table  entry 

Maximum  number  of  checkpoints  on  tape 

Loop  Index  over  number  of  checkpoints  read 

Flag  Indicating  end-of-flle  on  LUFILE 

Logical  unit  number  of  data  set 

Internal  variable  to  save  value  of  IRSTRT 

Loop  Index  over  module  names 

Array  of  keyword  numbers  describing  module 
names 

Present  position  of  data  file 

Internal  variable  set  to  MODCHK 

Array  of  Hollerith  format  module  names 

Internal  variable  equal  to  MODNAM 

Flag  set  to  ISOFF  so  that  commons  and  files 
are  read  from  MODCHK 


PREVIOUS  PAGE 
IS  BLANK 
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I/O  VARIABLES 
A.  INPUT 


LOCATION 


IERRF 
IOFILE 
IRSTRT 
ISOFF 
ISON 
KOLAST 
KOLFST 
KOLLOC 
LSTHOO 
LSTSYS 
LUPRNT 
MOOCHK 
HOOLST 
MOONAM 
NOATBL 
NPDATA 
B.  OUTPUT 
CHKHRT 
0B6PRT 
FRQHHZ 
IMOCHK 
IRSTRT 


/AOEBUG/ 

/IOFLES/ 

/AOEBUG/ 

/AOEBUG/ 

/AOEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/MOOULE/ 

/SYSFIL/ 

/AOEBUG/ 

/SYSFIL/ 

/MODULE/ 

/MODULE/ 

/PARTAB/ 

/PARTAB/ 

LOCATION 

/SYSFIL/ 

/AOEBUG/ 

/AMPZIJ/ 

/AOEBUG/ 

/AOEBUG/ 

/INTMAT/ 


KJFLD 


STRTUP 


(GTO,  MOM,  OUTPUT) 


KJGTO  /INTMAT/ 

KJMOM  /INTMAT/ 

NDATBL  /PARTAB/ 

NOSTAT  /ADEBUG/ 

RSTART  /SYSFIl/ 

6.  CALLING  ROUTINE: 

GEMACS 

7.  CALLED  ROUTINES: 

ASSIGN 

GETSYM 

MOVFIL 

PUTSYM 

RDEFIL 

RWCOMS 

RHFILS 

STATIN 

STATOT 

VILKBCK 


STRTUP 


(GTD,  MOM,  OUTPUT) 
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STRTUP 


(GTD,  MOM,  OUTPUT) 


1.  NAME:  SUBPAT 


(INPUT) 


2.  PURPOSE:  Augments  the  segment  table  when  a  wire-patch  connection  is 
found. 

3.  METHOD:  First,  the  subroutine  SUBPAT  determines  the  maximum  number 
of  entries  and  the  maximum  number  of  data  b  ocks  that  will  result 
from  all  wire-patch  connections.  Then,  the  segment  table  Is 
searched  for  the  first  patch  that  is  connected  to  a  wire  segment. 
This  patch  is  divided  Into  four  smaller  patches  of  ejgual  area  and 
oriented  with  respect  to  the  surface  vectors  t^  and  t2  as  shown  In 
figure  1. 


Figure  1.  Orientation  of  Subpatches  Connected  to  a  Wire  Segment 


Also,  the  unit  vector  of  each  of  the  smaller  patches  Is  the  same  as 
that  of  the  original  patch.  The  geometry  data  of  the  four 
subpatches  are  stored  in  a  temporary  array  along  with  all  succeeding 
patch  data  and  any  patches  augmented  by  a  connection  to  a  wire. 
After  all  patches  have  been  searched,  the  data  In  the  temporary 
array  are  stored  In  the  segment  table  starting  at  the  point  where 
the  first  wire-patch  connection  was  found.  Finally,  all  wire 
segments  are  checked  for  a  connection  to  a  patch.  If  such  a  connec¬ 
tion  Is  found,  the  connection  data  are  corrected  to  reflect  the  new 
patch  number  which  has  resulted  from  the  Increase  In  the  number  of 
patches  dlie  to  a  wire-to-patch  connection. 

4.  INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

AREA  Surface  area  of  patch 

Index  for  wire  segment  data  blocks 


IBLK 


The  Index  to  the  data  block  which  locates 
the  first  patch  connected  to  a  wire  segment 

Number  of  columns  In  the  connection  array 

Saved  value  of  ICOL 

Connection  data  for  wire  segment 

Connection  data  for  end  1  of  wire  segment 

Connection  data  for  end  2  of  wire  segment 

Flag  Indicating  whether  a  wire  segment 
connected  to  patch  has  been  found 

Logical  unit  on  which  symbol  Is  stored 

The  number  of  segments  In  requested  data 
block 

The  data  block  with  the  first  wire  segment- 
to-patch  connection 

Index  to  data  block  containing  the  Initial 
patch  data 

An  array  containing  current  position 
pointer  for  the  search  file 

Scratch  file  for  temporary  storage  of 
SEGTBL  data 

Integer  Identifying  which  patch  Is 
connected  to  end  1  of  wire  segment 

Integer  identifying  which  patch  Is 
connected  to  end  2  of  wire  segment 

The  number  of  patches  In  requested  data 
block 

The  location  of  first  patch  In  requested 
data  block 

Number  of  rows  In  connection  data  array 

Location  within  data  block  of  first  patch 
with  connection  to  a  wire  segment 


I WORDS 
JBIAS3 


MAXBLK 

MAXSEG 


MXBLKW 


NCONT 


NDXBLK 

NELMNT 


NNCON 


NOGOFG 


NPATCH 

NPRSEG 


Tag  Identifier 

Equlvalenced  to  WORDS 

Integer  to  bias  connection  data  to  a  patch 

The  number  of  connections  In  each  column  of 
connection  array 

Total  number  of  data  blocks 

Maximum  number  of  segments  per  data  block 

The  total  number  of  data  blocks  after 
accounting  for  wire  segment-to-patch 
connections 

The  total  number  of  data  blocks  containing 
wire  segments 

An  array  containing  the  old  and  new  patch 
numbers  for  a  wire  segment-to-patch 
connection 

The  number  of  wire  segment-to-patch 
connections 

Index  to  current  data  block 

The  original  number  of  wires  and  patches 
before  augmenting  for  wire  segment-to-patch 
connections 

Counter  for  a  new  location  of  patches  In 
SEGTBL 

Maximum  length  of  array  containing  old  and 
new  patch  numbers  for  wire  segment-to-patch 
connections 

No  go  flag 

Original  patch  number  before  patches  were 
divided 

The  total  number  of  patches 

The  number  of  data  Items  for  each  SEGTBL 
entry 


NUMSEG 

NUMWP 

NUIRE 

NXTBLK 

RH 

SIDE 

TIX.TIY.TIZ 

T2X,T2Y,T2Z 

WORDS 

XNPA,YNPA,ZNPA 

XPC,YPC,ZPC 

XSUBPA 

YSUBPA 

I/O  VARIABLES: 
A.  INPUT 
IOFILE 
I0SCR2 
IP217 
ISGTBL 


The  total  number  of  wire  segments  and 
patches  adjusted  for  wlre-to-patch 
connections 

Number  of  wire  segments  and  patches  In 
geometry 

The  total  number  of  wire  segments 

Index  to  the  next  data  block 

v/(XNPA  *  XNPA  +  YNPA  *  YNPA) 

The  x  or  y  dimension  for  the  distance  from 
the  center  point  of  the  patch  to  the  center 
point  of  a  subpatch 

A 

The  x,y,  and  z  components  of  tj 

A 

The  x,y,  and  z  components  of  t? 

Temporary  storage  array 

X,Y,  and  Z  components  of  patch  normal 
vector 

XfY,  and  Z  components  of  the  center  point 
of  the  patch 

The  x  coordinate  of  the  subpatch  with 

respect  to  the  patch  center  (=  SIDE) 

The  y  coordinate  of  the  subpatch  with 

respect  to  the  patch  center  (=  SIDE) 


LOCATION 

/IOFLES/ 

/SYSFIL/ 

/GEOOAT/ 

/SEGMNT/ 

/ADEBUG/ 


ISOFF 


SU8PAT  (INPUT) 


s-v- 


ISON 

JBIAS3 

KBREAL 

KOLCOL 

KOLLOC 

KOLNAM 

LUPRNT 

MAXBLK 

MAXSEG 

NAMSEG 

NCONT 

NOATBL 

NDXBLK 

NOPCOD 

NPATCH 

NPDATA 

NPRSEG 

NUMGTD 

NUMSEG 

NWIRE 

SEGTBL 

ZERO 

OUTPUT 

ISGTBL 

MAXBLK 


/ADEBUG/ 

/SEGMNT/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/SEGMNT/ 

/SEGMNT/ 

/SEGMNT/ 

F.P. 

/PARTAB/ 

/SEGMNT/ 

/ADEBUG/ 

/SEGMNT/ 

/PARTAB/ 

/SEGMNT/ 

/GTDDAT/ 

/SEGMNT/ 

/SEGMNT/ 

/SEGMNT/ 

/ADEBUG/ 

LOCATION 

/SEGMNT/ 

/SEGMNT/ 
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SUBPAT 


(INPUT) 


NOATBL 

/PARTAB/ 

NOGOFG 

/ADEBUG/ 

NPATCH 

/SEGHNT/ 

NUHSEG 

/SEGHNT / 

SEGTBL 

/SEGHNT/ 

UPDBLK 

/SEGHNT/ 

CALLING  ROUTINE: 

GEOORV 

CALLED  ROUTINES: 

ASSIGN 

MOVFIL 

STATIN 

CLSFIL 

OPNFIL 

STATOT 

GETSEG 

PUTSYM 

SYHDEF 

SYMDEF 

WLKBCK 

WRTFIL 


GETSYM 


RDEFIL 


SYMUPD 


Page 


NAME:  SYMDEF  (GTD,  INPUT,  MOM,  OUTPUT) 

PURPOSE:  To  define  or  redefine  a  symbol  during  program  execution. 

METHOD:  The  symbol  name  is  searched  for  in  the  NDATBL  array  and  if 
located  the  warning  message  is  printed  out  to  the  user  that  the 
symbol  is  being  redefined.  If  the  attributes  of  the  symbol  as 
defined  match  the  attributes  in  the  call  to  SYMDEF,  only  a  new 
edition  of  the  symbol  is  created,  not  a  completely  new  file.  Other¬ 
wise,  the  present  symbol  is  purged  from  the  symbol  table,  its  data 
file  closed,  and  an  entry  made  in  the  symbol  table  with  new 
attributes  as  specified  in  the  call.  Data  on  the  file  are  lost. 

If  not  located,  the  symbol  name  is  added  to  the  end  of  NOATBL,  and 
the  next  file  available  for  storage  is  assigned  to  the  symbol.  The 
file  is  opened,  and  the  NOATBL  pointers  are  reset.  Should  there  not 
be  a  file  available,  a  fatal  error  is  generated.  If  the  addition  of 
this  symbol  would  overflow  the  NDATBL  array,  a  fatal  error  is 
generated. 

INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

IBIT1  Bit  set  attributes  of  present  edition  of 

symbol 

IB ITS  Input  attribute  containing  the  bit  set 

attributes  of  the  symbol  being  defined 

INEW  Flag  indicating  new  data  set  (INEW=1)  or 

new  edition  of  present  data  set  (INEW=0) 

IOSTOR  Logical  unit  designator  of  file  for  this 

symbol 

LOCFST  Index  to  first  data  entry  of  this  edition 

of  symbol,  either  core  storage  or  file 
storage. 

LOCLST  Index  to  last  data  entry  of  this  edition  of 

symbol,  either  core  storage  or  file 
storage. 

NAME  User-assigned  name  of  symbol  to  be  defined 

or  redefined 


SYMOEF 


(GTD,  INPUT,  MOM,  OUTPUT) 


NCOL1 

NCOLS 

NEED 

NEUSYM 

NPDASV 

NROW1 

NROWS 

NSYMBL 

I/O  VARIABLES: 

A.  INPUT 
DBGPRT 
IBITS 
IOFILE 
IOSCR1 
I0SCR2 
IOSYMB 
I  PASS 
ISON 
KBCPLX 
KOLAST 
KOLBIT 
KOLCOL 


Number  of  columns  defined  for  present 
edition  of  symbol 

Number  of  columns  required  for  symbol  being 
defined  or  redefined 

Amount  of  In-core  storage  required  to  store 
data  In  FLTSYM 

Internal  variable  for  NAME 

Saved  value  of  variable  NPDATA  upon  entry 
to  SYMDEF 

Number  of  rows  defined  for  present  edition 
of  symbol 

Number  of  rows  required  for  symbol  being 
defined  or  redefined 

Number  of  active  entries  In  the  NDATBL 
array 

LOCATION 
/ADEBUG/ 

F.P. 

/IOFLES/ 

/SYSFIL/ 

/SYSFIL/ 

/SYSFIL/ 

/ARGCOM/ 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 


KOLFST 


SYMDEF 

(GTO,  INPUT,  MOM, 

OUTPUT) 

KOLLOC 

/PARTAB/ 

KOLNAM 

/PARTAB/ 

KOLROW 

/PARTAB/ 

LUPRNT 

/ADEBUG/ 

MAXSTR 

/SYMSTR/ 

NAME 

F.P. 

NCOLS 

F.P. 

NDATBL 

/PARTAB/ 

NDATMX 

/PARTAB/ 

NFILES 

/IOFLES/ 

NPDATA 

/PARTAB/ 

NROWS 

F.P. 

NXTSYH 

/SYMSTR/ 

B.  OUTPUT 

LOCATION 

IERRF 

/A0E8UG/ 

NDATBL 

/PARTAB/ 

NPDATA 

/PARTAB/ 

NXTSYH 

/SYMSTR/ 

6.  CALLING  ROUTINES*: 

BANDIT  (3) 

GEOORV  (1) 

SOLDRV  (3) 

DMPDRV  (1,2, 3,4) 

LODDRV  (3) 

SUBPAT  (1) 

EGFMAT  (3) 

LUOORV  (3) 

TSKXQT  (2,3,4) 

EXCDRV  (2,3) 

PUTSYM  (1,2, 3,4) 

ZIJDRV  (2,3) 

FLDDRV  (2,3,4) 

SETDRV  (3) 

*1- INPUT 

2- GTO 

3- MOM 

4- OUTPUT 

YMDEF 


(GTD,  INPUT,  MOM,  OUTPUT) 
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RESET  SYMBOL  TABLE 
POINTER 


IE 


'Sl 


2. 

PURPOSE:  Search  the  symbol  and  literal  tables  for  the  next  entry  in 

the  scan  table  and. 

if  found,  load  the  index  to  the  proper  table 

into  the  next  argument  list  table  entry. 

3. 

METHOD:  If  the  next 

scan  table  entry  is  an  alpha  field,  the  subrou- 

tine  SYMSCH  is  called  to  find  the  next  scan  table  entry.  If  the 

entry  is  not  found. 

an  error  condition  is  set;  if  it  is  found,  the 

index  to  the  symbol 

table  is  loaded  into  the  argument  list  table. 

If  the  next  scan  table  entry  is  not  an  alpha,  subroutine  LITSCH  is 

called  and  the  index 
table. 

is  returned  and  loaded  into  the  segment  list 

4. 

INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

INC 

Inline  function  to  position  NPARGL 

INDEX 

Index  to  symbol  table 

IN0EX1 

Index  to  literal  table 

5. 

I/O  VARIABLES: 

A.  INPUT 

LOCATION 

ISOFF 

/ADEBUG/ 

NARGMX 

/PARTAB/ 

NCOOE 

/SCNPAR/ 

NPARGL 

/PARTAB/ 

NPEARG 

/INPERR/ 

NPRSER 

/SCNPAR/ 

NTAB 

/SCNPAR/ 

NTALPH 

/ADEBUG/ 

NVAL 

/SCNPAR/ 
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v.-v:. 

B.  OUTPUT  LOCATION 

NARGTB  /PARTAB/ 

NPARGL  /PARTAB/ 

NPRSER  /SCHPAR/ 

NUMWRD  /ADEBUG/ 

CALLING  ROUTINES: 

FNOARG 

PARSE 

CALLED  ROUTINES 

ASSIGN  STATIN  SYMSCH 


FABL02 


STATOT 


WLKBCK 


NAME:  SYMMOD 


(MOM) 


PURPOSE:  This  subroutine  forms  a  symmetrical  combination  of  an 

Input  matrix  by  using  an  input  symmetry  operator. 

METHOD:  If  Z  is  a  matrix  containing  N  submatrices  and  S  is  the 

matrix  representation  of  the- symmetry  operator,  then  the  operation 

[Z]  =  [S]  [Z] 

forms  the  symmetrical  combinations  of  the  submatrices  Z. 

INTERNAL  VARIABLES: 


VARIABLE 


DEFINITION 


SYMOP 


Scratch  array 

Scratch  accumulation 

Number  of  columns  in  the  input  matrix 

Number  of  rows  in  the  input  matrix 

Dimension  of  the  symmetry  operator 

Index  pointer 


SYMOP 


S  matrix 


Input  matrix 


I/O  VARIABLES: 
A,  INPUT 


LOCATION 


previous  PAoe 
is  blank 


'D 


SYMMOD 


LOOP  OVER 
COLUMNS 
OF  INPUT 
MATRIX 


LOAD  SUBMATRIX 
ELEMENTS  INTO 
SCRATCH  ARRAY 


RETURN 


NAME:  SYMSCH 


(INPUT) 


PURPOSE:  Search  the  NDATBL  array  for  the  occurrence  of  the  symbol 
name  specified  in  the  argument  list. 

METHOD:  The  symbol  table  is  searched  for  the  occurrence  of  the  name 
specified  in  the  subroutine  argument  list.  If  the  name  is  found, 
and  was  not  supposed  to  be  previously  entered,  an  error  flag  is  set. 
If  the  name  is  not  found,  and  was  supposed  to  be  previously  entered, 
an  error  flag  is  set.  If  the  name  was  not  found  and  was  not  sup¬ 
posed  to  be  previously  entered,  it  is  entered  into  the  symbol  table 
and  the  index  is  returned  through  the  subroutine  argument  call. 

INTERNAL  VARIABLES: 


VARIABLES 

IEND 

INDEX 

NAME 

I/O  VARIABLES: 
A.  INPUT 
ISOFF 
ISON 
KOLNAM 
MATCH 
NAME 
NDATBL 
NOATMX 
NOMTCH 
NPDATA 
NPENOM 


DEFINITION 

Last  entry  in  the  symbol  table 
Index  to  the  symbol  table 
Symbol  name 


LOCATION 

/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 

/SCNPAR/ 

F.P. 

/PARTAB/ 

/PARTAB/ 

/SCNPAR/ 

/PARTAB/ 

/INPERR/ 

/INPERR/ 


NPESEX 


SYHSCH 


(INPUT) 


NPESYM 

/INPERR/ 

M 

^  A  a 

L'.i« 

NTAB 

/SCNPAR/ 

B.  OUTPUT 

LOCATION 

INOEX 

F.P. 

i 

NCOOE 

/SCNPAR/ 

NDATBL 

/PARTAB/ 

|5 

NPDATA 

/PARTAB/ 

H 

NPRSER 

/SCNPAR/ 

jfe 

NTAB 

/SCNPAR/ 

NUMWRD 

/ADEBUG/ 

CALLING  ROUTINES: 

FNOARG 

PLIST 

SYMLIT 

CALLED  ROUTINES: 
ASSIGN 
FABL02 
STATIN 


STATOT 

WLKBCK 


-A  ^  >.'• 


NAME:  SYMUPD 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  Update  attributes  of  entries  in  the  NOATBL  array. 

METHOO:  The  attribute  of  the  symbol  to  be  updated  is  called  through 
the  argument  list.  The  column  to  be  changed  in  the  NOATBL  array  is 
checked  for  validity  and,  if  invalid,  a  fatal  error  is  generated. 
If  valid,  the  appropriate  column  is  updated. 

INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

IFILE 

Logical  unit  on  which  the  symbol  is  stored 

KLM 

Input  argument  designating  column  of  NOATBL 
array  to  be  changed 

KOL 

Internal  variable  for  KLM 

NAMSYM 

Symbolic  name  for  symbol  to  be  updated 

NEWDAT 

Argument  containing  new  data  to  be  placed 
in  NOATBL  array 

NEWNAM 

Internal  representation  for  argument  NAMSYM 

NS 

Saved  value  of  loop  indexed  while  searching 
NDATBL  array 

NSYMBL 

Number  of  entries  in  the  NDATBL  array 

I/O  VARIABLES: 

A.  INPUT 

LOCATION 

DBGPRT 

/ADEBUG/ 

ISON 

/ADEBUG/ 

KLM 

F.P. 

KOLAST 

/PARTAB/ 

KOLBIT 

/PARTAB/ 

KOLCOL  /PARTAB/ 

KOLFST  /PARTAB/ 


SYMUPD 


(GTD,  INPUT,  MOM,  OUTPUT) 


KOLLNK 

/PARTAB/ 

;-v-- 

KOLLOC 

/PARTAB/ 

-  •  .  < 

KOLNAM 

/PARTAB/ 

KOLROW 

/PARTAB/ 

’*  "y 

LUPRNT 

/ADEBUG/ 

-  ■H 
.4 

J’*  _  , 

NATSYM 

F.P. 

V-p 

NOATBL 

/PARTAB/ 

NDFILE 

/IOFLES/ 

\'4 

NEMDAT 

F.P. 

NPDATA 

/PARTAB/ 

.  ■  •/ 

B.  OUTPUT 

LOCATION 

IERRF 

/ADEBUG/ 

•VA 

NDATBL 

/PARTAB/ 

:;:§ 

6.  CALLING  ROUTINES*: 

feegaj 

BANDIT  (3) 

PUTSEG  (1,2,3) 

EXCDRV  (2,3) 

SOLDRV  (3) 

_  < 

FLDDRV  (2,3,4) 

SUBPAT  (1) 

GEODRV  (1) 

TSKXQT  (1,2, 3,4) 

;Nv>: 

\s‘V 

LODDRV  (3) 

ZIJORV  (2,3) 

•  V  1 

LUDDRV  (3) 

*  •% 

7.  CALLED  ROUTINES: 

ASSIGN 

STATIN 

CONVRT 

STATOT 

ERROR 

WLKBCK 

■V  ■% 

♦1-INPUT 

2-GTD 

3-MOM 

4-OUTPUT 

D 


(GTD,  INPUT,  MOM,  OUTPUT) 


YES 


OUTPUT 

CALLING 

ARGUMENTS 


LOAD  NEWDAT 
IN  COLUMN 
SPECIFIED 


LOAD  1 
AS  FIRST 
WORD 
ADDRESS 


LOAD  CURRENT 
FILE  LENGTH  AS 
LAST  WORD 
ADDRESS 


NAME:  SYSCHK 


CM 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  Determine  If  time  for  checkpoint  has  passed. 

METHOD:  The  current  time  is  retrieved  and  if  less  than  the  next 
checkpoint  time,  control  is  returned  to  the  calling  subroutine.  If 
It  is  greater  than  the  next  checkpoint  time,  the  checkpoint  time  is 
incremented  and  subroutine  WRTCHK  is  called  to  write  a  checkpoint. 


INTERNAL  VARIABLES: 


VARIABLE 


ETIME 


FLTINC 

TIMCHK 


TLAST 


I/O  VARIABLES: 
A.  INPUT 
CHKPNT 
COMPLT 


INCCHK 


LUPRNT 

TIMTGO 


3.  OUTPUT 
COMPLT 
IERRF 


DEFINITION 
Delta  time 

Elapsed  time  from  last  time  check 
Total  elapsed  time  from  beginning  of  run 
Time  Increment 

Internal  variable  for  checking  time  elapsed 
Time  of  last  time  checkpoint 
Current  clock  time 


LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/SYSFIL/ 

/ADEBUG/ 

/ADEBUG/ 

/SYSFIL/ 

/ADEBUG/ 


LOCATION 

/SYSFIL/ 

/ADEBUG/ 


V>A'V- 
■  v  -  • 


i.-  ..  •.*  •/  *.'•."  .*■ 


5 


SYSCHK 


(GTD,  INPUT,  MOM,  OUTPUT) 


CALLING  ROUTINES*: 
OECOMP  (3) 

FLDORV  (2) 

TSKXQT  (1,2, 3, 4) 
ZIJDRV  (2,3) 

CALLED  ROUTINES: 

ASSIGN 

ERROR 

STATIN 

STATOT 

TICHEK 

WLKBCK 

WRTCHK 


♦1-INPUT 

2- GTD 

3- MOM 

4- OUTPUT 
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SYSCHK 


NAME:  SYSRTN 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  A  system-dependent  subroutine  to  return  various  auxiliary 
information  depending  on  local  subroutine  library  capability. 

METHOD:  The  quantity  desired  is  indicated  by  the  input  argument  I, 
and  the  quantity  to  be  returned  is  stored  in  argument  variable  J. 


INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 


DJ  =  IJ 
I 

IDATE 

ITIME 

J 

JHOURS 

JMINIT 

TIME 

NOTE: 


Time  of  day 

Input  argument  indicating  information 
desired 

Date  of  execution 

Intermediate  value  of  time  of  day  In 
minutes 

Output  argument  containing  information 
desired. 

Intermediate  value  of  time  of  day  (hours) 
Intermediate  value  of  time  of  day  (minutes) 
Time  of  day 

Explicit  form  of  this  subroutine  depends  on  local  library 
subroutines  available. 


I/O  VARIABLES: 

A.  INPUT  LOCATION 

I  F.P. 

B.  OUTPUT  LOCATION 

J  F.P. 

CALLING  ROUTINE: 

GEMACS 


PURPOSE:  To  compute  vectors  from  a  source  that  are  tangent  to  the 
cylinder  in  the  x-y  plane. 

METHOD:  The  unit  tangent  vectors  are  determined  by  solving  a  set  of 
equations  found  by  setting  the  incident  vector  from  the  source  equal 
to  the  general  unit  tangent  vector  to  the  elliptic  surface.  Details 
are  given  in  pages  90-93  in  reference  A.  General  tangents  and 
tangent  points  are  shown  in  figure  1. 


x 


XT2  =  J  A  cos  (VT(2>)  +  y  B  sin  (VT<2)) 
=  x  BT  (1)  +y  BT  (2) 
t2  =  xBT(3)  +  y  BT  (4) 


Figure  1.  Geometry  of  Source  Vectors  Tangent  to  the  Cylinder 
in  the  X-Y  Plane 


INTERNAL  VARIABLES: 


VARIABLE 

DEFINITION 

A 

Radius  of  cylinder  along  x  axis 

AA 

Distance  from  source  to  tangent  point 

AL 

Computational  variable 

B 

Radius  of  cylinder  along  y  axis 

BB 

Distance  from  origin  to  tangent  point 

BEI 

Computational  variable 

BT 

X  and  V  components  of  tangent  unit  vectors 
in  reference  coordinate  system 

CV 

Cosine  of  tangent  point  elliptical  angle 

CVE 

Cosine  of  V E 

DPR 

Degrees  per  radian  (  =  180. / -rr ) 

DT 

Dot  product  of  unit  vectors  of  the  two 
source  rays  tangent  to  the  cylinder  (2-D) 

DVl 

Angle  VI  in  degrees 

DV2 

Angle  V2  in  degrees 

El 

Error  detection  variable 

E2 

Error  detection  variable 

RHOE 

Distance  from  z  axis  to  point  where  ray 
from  origin  to  source  intersects  the 
cylinder 

RHOS 

Distance  from  source  to  z  axis 

SV 

Sine  of  tangent  point  elliptical  angle 

SVE 

Sine  of  VE 

SX 

X  component  of  ray  from  tangent  point  to 
source 

3 


Y  component  of  ray  from  tangent  point  to 
source 

X  component  of  tangent  ray  unit  vector 

(tangent  point  1) 

Y  component  of  tangent  ray  unit  vector 

(tangent  point  1) 

X  component  of  tangent  ray  unit  vector 

(tangent  point  2) 

Y  component  of  tangent  ray  unit  vector 

(tangent  point  2) 

Elliptical  angle  defining  tangent  point  1 

Elliptical  angle  defining  tangent  point  2 

Elliptical  angle  of  ray  from  origin  to 
source 

Elliptical  angle  defining  tangent  point 
location  in  RCS  x-y  plane 

Source  location 

X  component  of  tangent  point  location 
Computational  variable 

Y  component  of  tangent  point  location 


I/O  VARIABLES: 


INPUT 


LOCATION 


LUPRNT 


/GEOMEL/ 

/GEOMEL/ 

/PIS/ 

/ADEBUG/ 


'VVWV 


craw 


B.  OUTPUT 


LOCATION 


BT 

F.P. 

OT 

F.P. 

VT 

F.P. 

CALLING  ROUTINES: 

CYLINT 

GEOMC 

GEOMPC 

RPLSCL 

SCLRPL 

SCTCYL 

CALLED  ROUTINE: 

BTAN2 

REFERENCE: 

A.  R.  J.  Marhefka,  "Analysis  of  Aircraft  Wing-Mounted  Antenna  Pat¬ 
terns,"  Report  2902-25,  June  1976,  The  Ohio  State  University 
ElectroScience  Laboratory,  Department  of  Electrical  Engineer¬ 
ing;  prepared  under  Grant  No.  NGL  36-008-138  for  National  Aero¬ 
nautics  and  Space  Administration. 


RETURN 


NAME:  TICHEK 


(GTD,  INPUT,  MOM,  OUTPUT) 


PURPOSE:  To  obtain  the  current  clock  time  and  the  time  elapsed 
since  the  last  call  to  the  subroutine. 

METHOD:  The  time  is  initialized  to  zero  and  a  library  subroutine  is 
called  to  retrieve  the  current  processor  time.  The  elapsed  time  is 
computed  and  the  current  time  Is  saved  for  the  next  call. 

INTERNAL  VARIABLES: 


VARIABLE 

DEFINITION 

DT 

Output  argument  for  elapsed  time  since  last 
call 

T 

Output  argument  for  current  processor  time 

TLAST 

Saved  value  of  current  processor  time  to 
compute  DT  on  next  call 

TS 

Current  time  in  hours 

I/O  VARIABLES: 

A.  INPUT 

LOCATION 

None 

B.  OUTPUT 

LOCATION 

DT 

F.P. 

T 

F.P. 

CALLING  ROUTINES:* 

DECOMP  (3) 

WRTCHK  (1,2, 3, 4) 

SYSCHK  (1,2, 3,4) 

ZGTDRV  (2) 

TSKXQT  (1, 2,3,4) 

ZIJSET  (3) 

CALLED  ROUTINES: 

None. 

*1  -  INPUT  I  pREV»OUS 

2  -  GTD  \  is  blank — 

3  -  MOM 

4  -  OUTPUT 


(GTD,  INPUT,  MOM,  OUTPUT) 


Si 


TNtFLD 


(MOM) 


/  \  "jkr: 

'  V  * 


-  kA/2 ) 


}  sin 


sin  kA/2 
cos  kA/2 


+  (j  +  ^l)7^?(kz'  +  k4/2>  { - 

■  ‘o  (-1  2)^  TTT^  [ 


(P*  2’) 


sin  kA/2 
cos  kA/2 


cos  kA/2  \ 
-sin  kA/2  | 


-  (kz'  +  kA/2) 


-  (kz'  -  kA/2) 


cos  kA/2 


-jkr2  (  . 

,  e  Jsin 

+  — .  V 


kA/2 


kr^  I  sin  kA/2  j  kr2  |  cos  kA/2 


-jkr. 

(3  +  k'l)  TkTTr 


_  )sin  kA/2 
2  cos  kA/2 


Tc7~  | 'cos  kA/2  |  +  +  k4/2>2  (i  +  k^) 


-jkr. 


( k  r  x ) ' 


-sin  kA/2 
cos  kA/2 


The  expression  for  the  field  of  a  constant  current  distribution 
Involves  an  Integral  of  exp  (- jkr)/r  which  must  be  evaluated 
numerically.  The  field  components  for  a  current  I0  are: 


Ep  (P1*  *'>  -  X0 


TNEFLD 


(MOM) 


These  expressions  are  separated  into  real  and  imaginary  parts  for 
evaluation  In  the  program.  The  coordinate  p'  for  a  wire  segment  is 
taken  as  the  distance  from  the  observation  point  to  a  point  on  the 
side  of  the  segment  as  shown  in  figure  2. 


Figure  2.  Geometry  for  the  Determination  of  p' 


Also,  the  component  Ep  Is  multiplied  by  p/p'  to  account  for  the 
evaluation  1  n'lNEFLO .r8C*  CUrrent’  ,s  set  t0  °"e  for 


INTERNAL  VARIABLES: 


VARIABLE  DEFINITION 

B  Radius  of  wire  segment 

BK  kb 


CINT 

COINC 

CR1 

CR1R 

CR1RR 

CR2 

CR2R 

CR2RR 

CST 

ERIC 

ERIK 

ERIS 

ERRC 

ERRK 

ERRS 

ETA 

EZIC 

EZIK 

EZIS 

EZRC 


r  A/2 

/  cos  (kr)/r  dz 

J  -a/2 
p/p' 

(n/2)  cos  (kr^)/(krj) ,  (n  =  Jv0/to) 
CRl/(kr!) 

CRl/(krx)2 

(n/2)  cos  (kr2)/(kr2) 

CR2/(kr2) 

CR2/(kr2)2 
cos  (kA/2) 

Imaginary  part  of  Ep  for  cosine  current 
Imaginary  part  of  Ep  for  constant  current 
Imaginary  part  of  Ep  for  sine  current 
Real  part  of  Ep  for  cosine  current 
Real  part  of  Ep  for  constant  current 
Real  part  of  Ep  for  sine  current 

d  =  S^TTo 

Imaginary  part  of  Ez  for  cosine  current 
Imaginary  part  of  Ez  for  constant  current 
Imaginary  part  of  Ez  for  sine  current 
Real  part  of  Ez  for  cosine  current 


TNEFLD 


(MOM) 

Real  part  of  Ez  for  constant  current 
Real  part  of  Ez  for  sine  current 
Flag  for  numerical  integration 
A  flag  indicating  a  patch  observation  point 
P 

kp 
kp ' 

(kp*)2 

kri 

(kri)2 

kr2 

(kr2)2 

Length  of  segment 
/•A/2 

I  sin  (kr)/r  dz 

J- A/2 

kA/2 

(n/2)  sin  (kri)/(kri) 

SRl/(kn) 

SRl/(kri)2 

(n/2)  sin  (kr2)/(kr2) 

SR2/(kr2) 

SR2/(kr2)2 
sin  (kA/2) 


TNEFLD  (MOM) 


T1 

T1S 

T2 

T2S 

T3 

T3S 

T4 

T4S 


Temporary  storage  of  terms  in  electric 
field  expressions 


WAVLGH 

X 

WAVNUM 

k  =  2ir/X 

ZD1 

C\l 

< 

+ 

N 

ZD2 

kz'  -  kA/2 

ZP 

z' 

ZPK 

kz1 

11 

I/O  VARIABLES: 

n/2,  (n  =  J  y0/e 0) 

A.  INPUT 

LOCATION 

B 

/AMPZIJ/ 

ETA 

/AMPZIJ/ 

IJX 

F.P. 

ISOFF 

/ADEBUG/ 

RH 

F.P. 

S 

/AMPZIJ/ 

WAVLGH 

/AMPZIJ/ 

WAVNUM 

/AMPZIJ/ 

TNEFLD  (MOM) 


OUTPUT 

LOCATION 

ERIC, EZ IC 

F.P. 

ERIK.EZIK 

F.P. 

ERIS.EZIS 

F.P. 

ERRC.EZRC 

F.P. 

ERRK.EZRK 

F.P. 

ERRS.EZRS 

F.P. 

IJ 

/TMI/ 

I PATCH 

/TMI/ 

RHK 

/TMI/ 

RKB2 

/TMI/ 

ZPK 

/TMI/ 

CALLING  ROUTINES: 

NERFLD 

NTRPLT 

CALLED  ROUTINES: 

ASSIGN 

ROMBNT 

STATIN 

STATOT 

WLKBCK 

REFERENCE: 

A.  Stratton,  J.  A.  Electromagnetic  Theory.  McGraw  Hill  Book  Co. 
New  York,  1941,  p.  454. 


TNEFLD 


CALCULATE 

INDIVIDUAL 

TERMS 


CALCULATE 

E7  AND  E 
z  P 

FOR  SINE 


CALCULATE 

Ez  AND  Ep 

FOR  COSINE 

~T“ 


CALCULATE 

E  FOR 
P 

CONSTANT 

CALL  ROMBNT 
FOR  CONSTANT 
CURRENT 

~r 


CALCULATE 
Ez  FOR 

CONSTANT 


I 


RETURN 


1.  NAME:  TNHFLD 


(MOM) 


2.  PURPOSE:  To  compute  the  near  magnetic  field  at  a  point  in  space  due 
to  the  current  on  a  wire  segment.  The  field  is  computed  for  three 
current  distributions:  sine,  cosine,  and  constant  functions,  each 
of  unit  amplitude. 

3.  METHOD:  The  wire  segment  is  considered  to  be  located  at  the  origin 
of  a  local  cylindrical  coordinate  system  with  the  point  at  which  the 
field  is  computed  being  (p1,  4>',  z‘).  The  geometry  for  a  filament 
of  current  of  length  A  is  shown  in  figure  1. 


segment 


Figure  1.  Geometry  for  Fields  Due  to  a  Filament  of  Current 


For  a  sine  or  cosine  current  distribution  the  field  can  be  written 
in  closed  form.  The  p  field  component  for  a  current 


TNHFLD  (MOM) 


"♦  (P'-  2'>  “4^  {e«P(-jkr2)  -  expl-jkr^ 

G 


j ( z '  -  A/2 ) 


exp( - jkr2 ) 


s in ( k A/2 ) 
cos( kA/2 ) 


exp{ - jkr .  ) 

+  j  (  z  1  +  A/2)  - - - — 

rl 


-sin(kA/2)  ) 
cos(kA/2)  / 


where  I0  =  1  is  assumed  in  this  routine. 

For  small  values  of  p  with  |z'|  >  A/2,  this  equation  may 
produce  large  numerical  errors  due  to  cancellation  of  large  terms. 
Hence,  for  z1  >  0  and  p'/(z'±A/2)  <  10~3t  a  more  stable  approxima¬ 
tion  for  small  p ' /(z ' ±A/2)  is  used: 


Ip' 

H<j,  (p>'  z,)  =  “alT  exPMkz’) 


( z ’ +A/2 )  ( z 1 -A/2 ) 


[-.] 


exp( jkA/2 ) 

/ sin  ( kA/2 ) \ 

_  expl- jkA/2) 

-J 

/-sin  (kA/2)\ 

\ 

_(z'-A/2)2 

\  cos  (kA/2 )  ) 

( z ' +A/2 ) 2 

\  cos  (kA/2  )J 

i 

For  z'  <  0,  the  above  equation  is  evaluated  for  H^(p',  -  z').  The 
field  of  a  sin  kz  current  is  multiplied  by  -1  in  this  case,  since  it 
is  an  odd  function  of  z. 

The  field  due  to  a  constant  current  is  obtained  by  numerical 
integration,  which  is  performed  by  subroutine  ROMBNT.  If 

p1  is  zero,  all  field  quantities  are  set  to  zero,  since 

is  undefined. 

4.  INTERNAL  VARIABLES: 


VARIABLE 

CDHK 

CONST 

C0NST1 

GNST2 


DEFINITION 
cos(kA/2) 
l/(4irp  1 ) 

P '  /8tt 
p '  k^/4-n 


1200 


TNHFLD 


(MOM) 


HPIC.HPRC 


HPIK,HPRK 


HPIS.HPRS 


I  PATCH 


cos(kri) 
cos(kr2) 
cos(kz' ) 

A/2 

kA/2 

The  imaginary  and  real  parts  of  the  <J> 
component  of  the  magnetic  field  due  to  a 
cosine  current 

The  imaginary  and  real  parts  of  the  <\> 
component  of  the  magnetic  field  due  to  a 
constant  current 

The  imaginary  and  real  parts  of  the  4> 
component  of  the  magnetic  field  due  to  a 
sine  current 

The  sign  of  z 

A  flag  indicating  that  the  magnetic  field 
at  a  patch  is  to  be  computed 

P1 

p'k 

P7(z'  -  A/2) 

rl 

kri 

r2 

kr2 

A,  the  segment  length 
sin  (kA/2) 
sin  (krjJ 
sin  (kr2) 


1201 


SZPK 

T1  \ 

T1K  I 

T2  / 

T2  ) 

T2S  \ 

T3C  I 

T3S  / 

WAVNUM 

ZD1 

ZD2 

ZP 

ZPK 

ZPSV 

I/O  VARIABLES 


sin  (kz1) 


Temporary  storage  variables  used  in 
Ccomputing  the  magnetic  field 


2tt/X 
z'  +  A/2 
z'  -  A/2 
z1 
kz1 

Save  value  of  z1 


INPUT 

LOCATION 

ISON 

/AOEBUG/ 

RH 

F.P. 

S 

/AMPZIJ/ 

TWOPI 

/AMPZIJ/ 

WAVNUM 

/AMPZIJ/ 

ZP 

F.P. 

OUTPUT 

LOCATION 

HPIC.HPRC 

F.P. 

HPIK,HPRK 

F.P. 

-  A''* 
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TNHFLO 


(MOM) 


HPIS.HPRS 
I PATCH 


/TMI  / 
/TMI/ 
/TMI/ 


CALLING  ROUTINE: 
NTRPLT 

CALLED  ROUTINES: 


ASSIGN 


ROMBNT 


STATIN 

STATOT 


WLKBCK 


1 


NAME:  TPNFLD 


(GTD) 


2.  PURPOSE:  To  calculate  the  theta  and  phi  unit  vectors  for  the  near¬ 
field  observation  direction. 

3.  METHOD:  Vector  algebra  is  used  to  compute  the  two  unit  vectors. 
Figure  1  shows  the  geometry  required. 


4.  INTERNAL  VARIABLES: 
VARIABLE 
DP 
DT 
PH 
TH 


DEFINITION 
Phi  unit  vector 
Theta  unit  vector 
Phi  angle 
Theta  angle 


5.  I/O  VARIABLES: 

A.  INPUT 
PH 

TH 

B.  OUTPUT 
DP 

DT 


LOCATION 

F.P. 

F.P. 

LOCATION 

F.P. 

F.P. 


6.  CALLING  ROUTINES: 


DIFPLT 

END  IF 

REFCYL 

DPLRCL 

RCLDPL 

RPLDPL 

DPLRPL 

RCLRPL 

RPLRCL 

SCLRPL 

7.  CALLED  ROUTINE: 
NONE 
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COMPUTE  3  COMPONENTS 
FOR  THETA 


COMPUTE  2  COMPONENTS 
FOR  PHI 


NAME:  TRCEBK  (GTD,  INPUT,  MOM,  OUTPUT) 

PURPOSE:  To  print  out  the  table  of  subroutines  generated  by  WLKBCK 
for  locating  a  fatal  error. 

METHOD:  Prints  out  the  table  of  subroutines  called  before  the  fatal 

err0r  was. grated  by  WLKBCK.  This  table  is  contained  in 

array  NAMRTN  and  is  indexed  by  INDXWB. 


INTERNAL  VARIABLES: 
VARIABLE 

11 

12 

NAMSUB 

I/O  VARIABLES: 


DEFINITION 

Internal  variable  set  to  Hollerith  name  of 
subroutine  for  error  output  message 

Internal  variable  set  to  Hollerith  name  of 
subroutine  for  error  output  message 

Internal  variable  set  to  Hollerith  name  of 
subroutine 


i 


sS 

M 

A.  INPUT 

LOCATION 

1 

INDXWB 

/ADEBUG/ 

& 

LUPRNT 

/ADEBUG/ 

NAMRTN 

/AOEBUG/ 

Ip 

B.  OUTPUT 

LOCATION 

• .  *, 

INDXWB 

/ADEBUG/ 

<*;«•. 

Caui 

6.  CALLING  ROUTINES: 

m 

ERROR 

!  •  *  « 

WLKBCK 

m 

7.  CALLED  ROUTINES: 

NONE 

1 


PREvipuTp^~ 

IS  BUMK 
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(GTD,  INPUT,  MOM,  OUTPUT) 


1.  NAME:  TRNLAT 


(INPUT) 


PURPOSE:  Translates  a  point  to  or  from  the  origin  of  a  coordinate 
system. 

METHOD:  The  point  is  translated  along  its  cartesian  coordinate 
axes,  and  the  operation  code  specifies  whether  it  is  a  translation 
to  or  from  the  origin. 


INTERNAL  VARIABLES: 


VARIABLE 


I/O  VARIABLES: 
A.  INPUT 


DEFINITION 


The  amount  of  the  translation  along  x  axis 

The  amount  of  the  translation  along  y  axis 

The  amount  of  the  translation  along  z  axis 

Translation  operation  code.  If  greater 

than  zero,  a  translation  to  origin.  If 
less  than  zero,  a  translation  from  the 
origin 

Input/output  of  x  coordinate 
Input/output  of  y  coordinate 
Input/output  of  z  coordinate 


LOCATION 


->  -,v  vN*N-  ■  ;  •  v 


TRNLAT 


(INPUT) 


B.  OUTPUT 


LOCATION 


CALLING  ROUTINES: 
COORDS 


WYRDRV 


CALLED  ROUTINES: 
ASSIGN 
STATIN 
STATOT 


WLKBCK 


*•  .  -  A  . 


NAME:  TSKXQT  (GTD) 


PURPOSE:  To  read  the  task  list  and  call  the  appropriate  processors 
to  execute  the  tasks. 

METHOD:  The  task  list  is  scanned  twice:  during  the  first  scan  the 
subroutines  necessary  to  execute  the  tasks  are  called  in  order  to 
initialize  the  required  parameters.  During  the  second  pass  the 
subroutines  are  called  to  perform  the  tasks  as  specified  by  the 
user. 

Task  execution  normally  begins  with  the  first  task  in  the  task  list 
and  proceeds  sequentially  through  the  list  unless  a  LABEL  task  is 
encountered.  The  LABEL  task  will  redirect  execution  to  its  associa¬ 
ted  LOOP  task  until  the  required  number  of  LOOP/LABEL  loops  has  been 
fulfilled.  Task  execution  terminates  when  an  END  command  is  encoun¬ 
tered,  the  end  of  the  task  list  is  reached,  or  an  error  occurs  in 
executing  a  task. 

If  the  task  list  has  been  generated  by  a  RSTART  command,  execution 
may  not  necessarily  begin  at  the  top  of  the  task  list.  Normally, 
restart  is  begun  from  the  task  which  wrote  the  checkpoint  read  in  to 
generate  the  task  list.  In  modules  subsequent  to  the  one  which  gen¬ 
erated  the  checkpoint,  execution  can  begin  at  the  top  of  the  task 
list  (if  the  preceeding  run  did  not  complete  its  execution)  or  at 
the  restart  task  (if  execution  was  successfully  completed). 

The  following  tasks  are  active  in  the  GTD  module: 


FORTRAN 

LABEL 

TASK 

NAME 

GTD  MODULE  FUNCTION 

120 

BACSUB 

Link  data  set  of  solution  vector  to  inter¬ 
action  matrix  data  set  and  identify  it  as  a 
solution  data  set 

130 

BAND 

Link  banded  matrix  data  set  to  full  matrix 
data  set 

150 

CHKPNT 

Retrieve  timed  checkpoint  parameters  or 
write  a  command  checkpoint 

180 

DEBUG 

Turn  off  or  on  the  debug  flags 

190 

DECOMP 

Link  data  sets  of  decomposed  matrix  to  its 
parent  data  set 

END 

Terminate  module  execution 

250  BHI 


260 

LABEL 

270 

LOOP 

390 

RESTRT 

410 

SOLVE 

440 

WIPOUT 

480 

GMDATA 

490 

ZGEN 

530 

EFIELD 

540 

DMP 

550 

ESRC.VSRC 

570 

SETINT 

INTERNAL  VARIABLES: 

VARIABLE 

CPFRWD 

DBGPRT 

DT 

IBIT 

I  BITS 

INCCHK 


I5>KXQI  (ulD) 


Link  data  set  of  solution  vector  to  Inter¬ 
action  matrix  data  set  and  Identify  It  as  a 
solution  data  set 

Decrement  loop  counter  and  branch  to  LOOP 
If  positive 

Initialize  loop  counter 

Process  RSTART  command  error 

Link  solution  to  excitation  and  excitation 
to  Interaction  matrices 

Process  W I POUT  command  error 

Advance  edition  of  geometry  data  set  and 
reinitialize  GTD  geometry  data 

Call  ZIJDRV  to  generate  GTD  interactions 

Call  FLDDRV  to  generate  incident  field 
matrix  and  scattered  field  Green's  function 
matrix 

Call  DMPDRV  to  process  direct  manipulations 

Call  EXCORV  to  generate  GTD  excitation 

Call  SET  to  select  GTD,  MOM,  and  incident 
field  interactions 


DEFINITION 

Checkpoint  file  rewind  flag 
Debug  print  flag 

Time  Interval  between  calls  to  TICHEK 
Attribute  word 

Attribute  word  for  geometry  data  set 
Checkpoint  time  increment 


1216 


6 


Ml 

Ml 


mi 


SI 

Ml 

M 

>2 1 

m 


M 

I 

If 


INDXA.INDXB,  INOXC 
INDXG,  INOXX 

IOCKPT 

ITASK 

JTASK 

KOUNT 


LINDX 


LINKA 


LOCARG 

LOCNXT 


LOCTP1 

LOCTSK 

LSTARG 


LSTTPF 


LTRACE 


NAMDAT 


NAME B.NAMEC, NAME X 
NAMGEO 


Pointer  to  symbol  table  entry  for  a  data 
set 

Logical  unit  number  of  checkpoint  file 

Pointer  to  task  in  task  list  being  executed 

Internal  variable  equal  to  ITASK 

The  number  of  times  the  loop  terminating  on 
the  reference  label  has  been  executed 

Index  to  the  loop  table  entry  currently 
being  executed 

Pointer  to  data  set  linked  to  data  set 
pointed  to  by  INDXA 

Pointer  to  task  argument  in  NARGTB 

Pointer  to  NARGTB  for  the  next  task  to  be 
executed 

Pointer  to  first  argument  for  a  given  task 

Locution  of  task  parameter  in  NARGTB 

Location  In  argument  list 

Pointer  to  last  task  executed  for  a  restart 
job 

Trace  flag  for  debug 
Loop  Index 

User-assigned  name  of  geometry  data  set 

User-assigned  names  of  INTARG  data  sets 

Pointer  to  default  geometry  data  set  name 
In  NCODES 

Index  to  NCODES  array  for  the  task  name 
mnemonic 

No  operation  flag 


*  ■/.  ,  ,,kuW,t.‘Awa 


. .  .A’VA/. -V  V 


*•  v. 

I  ’  O  i  *  «  *  »  S  ’  *  »  *  M  ■  «  »  «  t  .  »  . 


TSKXQT 


(GTD) 


NOSTAT 

Logical  flag  set  if  statistics  have 
been  requested 

not 

NPRREC 

Number  of  words  per  geometry  data  set  used 

NT 

Hollerith  name  of  task 

NUMARG 

Number  of  INTARG  arguments  for  a  task 

NUMTSK 

Task  Identification  number 

NXTTSK 

Pointer  to  the  next  task  to  be  executed 

TNOW 

Current  processor  time 

TRACST 

Logical  flag  set  If  trace  statistics 
desired 

are 

YSSTAT 

Logical  flag  set  If  statistics  have 
requested 

been 

I/O  VARIABLES: 

A.  INPUT 

LOCATION 

CHKPNT 

/SYSFIL/ 

CHKWRT 

/SYSFIL/ 

COMPLT 

/SYSFIL/ 

D8GPRT 

/ADEBUG/ 

ISOFF 

/ADEBUG/ 

ISON 

/ADEBUG/ 

KBGEOM 

/PARTAB/ 

KBREAL 

/PARTAB/ 

KBSOLN 

/PARTAB/ 

KOLBIT 

/PARTAB/ 

KOLCNT 

/PARTAB/ 

KOLCOL 

/PARTAB/ 
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TSKXQT 


KOLLNK 

KOLNAH 

KOLTIM 

KOLTSK 

KWOFF 

KWON 

KWSTAT 

KWTRAC 

LOOPMX 

LSTTPF 

LUPRNT 

NAXBLK 

MAXSEG 

MXARGS 

NAMTSK 

NARGTB 

NCOOES 

NDATBL 

NLOOPS 

NOGOFG 

NOPCOO 

NOSTAT 

NPRSEG 

NPTASK 


(STD) 


/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/ADEBUG/ 

/SEGMNT/ 

/SEGMNT/ 

/ARGCOM/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 

/SEGMNT/ 

/ADEBUG/ 


NTINT 


M  ^  V  V 


NTSKTB 

NXTTSK 

RSTART 

RSTRTA 


f 


OUTPUT 

CHKWRT 

CPFRWD 

DBGPRT 

IERRF 

INCCHK 


INTARG 


IOCKPT 


IPASS 


L* 

•X. 


LSTTPF 

MAXBLK 


NOGOFG 


HOSTAT 

HUMARG 

NXTTSK 

RSTART 


6.  CALLING  ROUTINE: 


TSKXQT  (GTD) 

/PARTAB/ 

/ADEBUG/ 

/SYSFIL/ 

/SYSFIL/ 

LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/AOEBUG/ 

/ADEBUG/ 

/SYSFIL/ 

/ARGCOM/ 

/SYSFIL/ 

/AOEBUG/ 

/SYSFIL/ 

/SEGMNT/ 

/AOEBUG/ 

/AOEBUG/ 

/ARGCOM/ 

/AOEBUG/ 

/SYSFIL/ 


GEMACS 


TSKXQT  (GTD) 


CALLED  ROUTINES: 

ASSIGN 

GETGEO 

SYSCHK 

CONVRT 

GTDDRV 

TICHEK 

DMPORV 

SET 

WLKBCK 

ERROR 

STATIN 

WRTCHK 

EXCORV 

STATOT 

ZIJDRV 

FLOORV 

SYMDEF 

GETARG 

SYMUPD 

NAME:  TSKXQT  (INPUT) 


PURPOSE:  To  read  the  task  list  and  call  the  appropriate  processors 
to  execute  the  tasks. 

METHOD:  The  task  list  is  scanned  twice:  during  the  first  scan  the 
subroutines  necessary  to  execute  the  tasks  are  called  in  order  to 
initialize  the  required  parameters.  During  the  second  pass  the 
subroutines  are  called  to  perform  the  tasks  as  specified  by  the 
user. 

Task  execution  normally  begins  with  the  first  task  In  the  task  list 
and  proceeds  sequentially  through  the  list  unless  a  LABEL  task  Is 
encountered.  The  LABEL  task  will  redirect  execution  to  its  associa¬ 
ted  LOOP  task  until  the  required  number  of  LOOP/LABEL  loops  has  been 
fulfilled.  Task  execution  terminates  when  an  END  command  is  encoun¬ 
tered,  the  end  of  the  task  list  Is  reached,  or  an  error  occurs  In 
executing  a  task. 

If  the  task  list  has  been  generated  by  a  RSTART  command,  execution 
may  not  necessarily  begin  at  the  top  of  the  task  list.  Normally, 
restart  Is  begun  from  the  task  which  wrote  the  checkpoint  read  in  to 
generate  the  task  list.  In  modules  subsequent  to  the  one  which 
generated  the  checkpoint,  execution  can  begin  at  the  top  of  the  task 
list  (If  the  preceding  run  did  not  complete  its  execution)  or  at  the 
restart  task  (If  execution  was  successfully  completed). 

The  following  tasks  are  active  in  the  INPUT  module: 


FORTRAN 

TASK 

LABEL 

NAME 

INPUT  MODULE  FUNCTION 

120 

BACSUB 

Link  data  set  of  solution  vector  to  data 
set  of  Interaction  matrix 

130 

BAND 

Link  banded  matrix  data  set  to  full  matrix 
data  set 

150 

CHKPNT 

Recover  the  checkpoint  file  rewind  flag 

180 

DEBUG 

Turn  off  or  on  the  debug  flags 

190 

DECOMP 

Link  data  sets  of  decomposed  matrix  to  Its 
parent  data  set 

200 

END 

Terminate  module  execution 

230 

INPUT 

Call  GEODRV  to  process  Input  data 

TSKXQT 


(INPUT) 


250 

BMI 

Link  data  set  of  solution  vector  to  data 
set  of  interaction  matrix 

260 

LABEL 

Decrement  loop  counter  and  branch  to  LOOP 
if  positive 

270 

LOOP 

Initialize  loop  counter 

390 

RESTRT 

Process  RSTART  command  error 

410 

SOLVE 

Link  data  set  of  solution  vector  to  data 
set  of  Interaction  matrix 

440 

WIPOUT 

Process  WIPOUT  command  error 

480 

GiMDATA 

Call  GEODRV  to  process  geometry  Input 

490 

ZGEN 

Link  data  set  of  Interaction  matrix  to 
geometry  data  set 

530 

EFIELD 

Call  EFDGEO  to  assure  that  EFIELD  argument 
is  linked  to  geometry  data  set 

540 

DMP 

Call  DMPDRV  to  process  direct  manipulations 

550 

ESRC.VSRC 

Link  data  set  of  excitation  vector  to  geom¬ 
etry  data  set 

INTERNAL  VARIABLES: 

VARIABLE  DEFINITION 

CPFRWD  Checkpoint  file  rewind  flag 

DBGPRT  Debug  print  flag 

DT  Time  Interval  between  calls  to  TICHEK 

INDEX  Pointer  to  geometry  data  set 

INDXAtINDXB(INDXX  Pointer  to  symbol  table  entry  for  a  data 

set 

IOCKPT  Logical  unit  number  of  checkpoint  file 

Temporary  storage  for  argument  2  of  INTARG 
array 


ISAV2 


TSKXQT 


(INPUT) 


ISAV3 

ITASK 

KOUNT 

LINDX 

LINKA 

LOCARG 

LOCNXT 

LOCTPl 

LOCTSK 

LSTARG 

LSTTPF 

LTRACE 

N 

NAMEA,NAMEB 

NAMGEO 

NOX 

NOP 

NOSTAT 

NPRREC 

NT 


Temporary  storage  for  argument  3  of  INTARG 
array 

Pointer  to  task  in  task  list  being  executed 

The  number  of  times  the  loop  terminating  on 
the  reference  label  has  been  executed 

Index  to  the  loop  table  entry  currently 
being  executed 

Pointer  to  data  set  linked  to  data  set 
pointed  to  by  INDXA 

Pointer  to  task  argument  in  NARGTB 

Pointer  to  NARGTB  for  the  next  task  to  be 
executed 

Pointer  to  first  argument  for  a  given  task 

Location  of  task  parameter  in  NARGTB 

Location  in  argument  list 

Pointer  to  last  task  executed  for  a  restart 
job 

Trace  flag  for  debug 
Loop  Index 

NAMEX  User-assigned  names  of  INTARG  data  sets 

Pointer  to  default  geometry  data  set  name 
In  NCODES 

Index  to  NCODES  array  for  the  task  name 
mnemonic 

No  operation  flag 

Logical  flag  set  If  statistics  have  not 
been  requested 

Number  of  words  per  geometry  data  set 
record 

Hollerith  name  of  task 


TSKXQT 


(INPUT) 


NUMARG 

Number  of  INTARG  arguments  for  a  task 

NUMTSK 

Task  Identification  number 

NXTTSK 

Pointer  to  the  next  task  to  be  executed 

TNOW 

Current  processor  time 

TRACST 

Logical  flag  set  if  trace  statistics 
desired 

YSSTAT 

Logical  flag  set  if  statistics  have 
requested 

I/O  VARIABLES: 

A.  INPUT 

LOCATION 

CHKWRT 

/SYSFIL/ 

COMPLT 

/SYSFIL/ 

DBGPRT 

/ADEBUG/ 

ISOFF 

/ADEBUG/ 

ISON 

/ADEBUG/ 

KOLCNT 

/PARTAB/ 

KOLLNK 

/PARTAB/ 

KOLNAM 

/PARTAB/ 

KOLTIM 

/PARTAB/ 

KOLTSK 

/PARTAB/ 

KWGEOM 

/PARTAB/ 

KWNAME 

/PARTAB/ 

KWOFF 

/PARTAB/ 

KWON 

/PARTAB/ 

KWSTAT 

/PARTAB/ 
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KWTRAC 


LOOPMX 
LSTTPF 
LUPRNT 
MXARGS 
NAMTSK 
NARGTB 
NCODES 
NDATBL 
N LOOPS 
NOGOFG 
NOPCOD 
NOSTAT 
NPTASK 
NTSKTB 
NXTTSK 
RSTART 
RSTRTA 

OUTPUT 

CHKPNT 

CHKWRT 

CPFRWD 

DBGPRT 

IERRF 

INTARG 


/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/ADEBUG/ 

/ARGCOM/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/SYSFIL/ 

/SYSFIL/ 

LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/SYSFIL/ 

/ADEBUG/ 

/ADEBUG/ 

/ARGCOM/ 

/SYSFIL/ 


IOCKPT 


TSKXQT 


(INPUT) 


I  PASS 

/ADEBUG/ 

LSTTPF 

/SYSFIL/ 

LTRACE 

/ADEBUG/ 

NOSTAT 

/ADEBUG/ 

NUHARG 

/ARGCOM/ 

RSTART 

/SYSFIL/ 

TRACST 

/ADEBUG/ 

CALLING  ROUTINE: 

GEMACS 

CALLED  ROUTINES: 

ASSIGN 

STATIN 

CONVRT 

STATOT 

DMPDRV 

SYHUPO 

EFDGEO 

SYSCHK 

ERROR 

TICHEK 

GEOORV 

WLKBCK 

ZZXDUM 


*-ASL 


;vai 


NAME:  TSKXQT 


(MOM) 


PURPOSE:  To  read  the  task  list  and  call  the  appropriate  processor 
to  execute  the  tasks. 

METHOD:  The  task  list  Is  scanned  twice:  during  the  first  scan  the 
subroutines  necessary  to  execute  the  tasks  are  called  in  order  to 
initialize  the  required  parameters.  During  the  second  pass  the 
subroutines  are  called  to  perform  the  tasks  as  specified  by  the 
user. 

Task  execution  normally  begins  with  the  first  task  in  the  task  list 
and  proceeds  sequentially  through  the  list  unless  a  LABEL  task  is 
encountered.  The  LABEL  task  will  redirect  execution  to  its  asso¬ 
ciated  LOOP  task  until  the  required  number  of  LOOP/LABEL  loops  has 
been  fulfilled.  Task  execution  terminates  when  an  END  command  is 
encountered,  the  end  of  the  task  list  is  reached,  or  an  error  occurs 
In  executing  a  task. 

If  the  task  list  has  been  generated  by  a  RSTART  command,  execution 
may  not  necessarily  begin  at  the  top  of  the  task  list.  Normally, 
restart  is  begun  from  the  task  which  wrote  the  checkpoint  read  in  to 
generate  the  task  list.  In  modules  subsequent  to  the  one  which 
generated  the  checkpoint,  execution  can  begin  at  the  top  of  the  task 
list  (If  the  preceeding  run  did  not  complete  its  execution)  or  at 
the  restart  task  (If  execution  was  successfully  completed). 

The  following  tasks  are  active  in  the  MOM  module. 


FORTRAN 


MOM  MODULE  FUNCTION 


BACSUB  Call  SOLDRV  to  back  substitute  to 

find  solution  vector 

BAND  Call  BANDIT  to  band  a  matrix 

CHKPNT  Retrieve  timed  checkpoint  parameters 

or  write  a  command  checkpoint 

DEBUG  Turn  off  or  on  the  debug  flags 

OECOMP  Call  LUDDRV  to  decompose  matrix  Into 

upper  and  lower  triangular  matrices 

END  Terminate  module  execution 

BMI  Call  SOLDRV  to  perform  banded  matrix 

iteration 


PREVIOUS  PAGE 
IS  BLANK 


TSKXQT 


(MOM) 


260 

LABEL 

Decrement  loop  counter  and  branch  to 
LOOP  if  positive 

270 

LOOP 

Initialize  loop  counter 

310 

WRITE 

Call  PRTSYM  to  write  symbol  data  to 
output  file 

340 

PRINT 

Call  PRTSYM  to  write  symbol  data  to 
output  file 

350 

PURGE 

Purge  symbol  from  NDATBL  and  close 
data  file 

390 

RESTART 

Process  RSTART  command  error 

400 

SET 

Call  SETDRV  to  set  data  set  entries 

410 

SOLVE 

Call  LUDORV  and  SOLORV  to  obtain 
solution  vector 

440 

WIP0UT 

Process  WIPOUT  command  error 

460 

ZSET 

Call  SETDRV 

480 

GMDATA 

Advance  edition  of  geometry  data  set 

490 

ZGEN 

Call  ZIJDRV  to  generate  MOM  inter¬ 
action  matrix 

520 

ZLOAOS 

Call  LODORV  to  generate  load  vector 

530 

EFIELD 

Call  FLDDRV  to  generate  total  field 
from  incident  and  scattered  fields 

540 

DMP 

Call  DMPDRV  to  process  direct 
manipulations 

550 

ESRC,  VSRC 

Call  EXCDRV  to  generate  MOM  excita¬ 
tion 

570 

SETINT 

Call  SET  to  select  GTD,  MOM,  and 
Incident  field  Interactions 

INTERNAL  VARIABLES: 

VARIABLE 

DEFINITION 

CPFRMO 

Checkpoint  file  rewind  flag 

D6GPRT 

Debug  print  flag 
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OT 

I 

IBITS 

INCCHK 

INOXG 

IOCKPT 

ISAV2 

ISAV3 

ITASK 

KOUNT 

LINDX 

LOCARG 

LOCFIL 

LOCNXT 

LOCSYM 

LOCTP1 

L0CTP2 

LOCTSK 

LSTARG 

lsttpf 


TSKXQT  (MOM) 


Time  Interval  between  calls  to  TICHEK 
Loop  Index 

Attribute  word  for  geometry  data  set 

Checkpoint  time  Increment 

Pointer  to  symbol  table  entry  for  a  data 
set 

Logical  unit  number  of  checkpoint  file 

Temporary  storage  for  argument  2  of  INTARG 
array 

Temporary  storage  for  argument  3  of  INTARG 
array 

Pointer  to  task  In  task  list  being  executed 

The  number  of  times  the  loop  terminating  on 
the  reference  label  has  been  executed 

Index  to  the  loop  table  entry  currently 
being  executed 

Pointer  to  task  argument  In  NARGTB 

Logical  file  associated  with  symbol  to  be 
purged 

Pointer  to  NARGTB  for  the  next  task  to  be 
executed 

Location  pointer  for  a  symbol  name 

Pointer  to  first  argument  for  a  given  task 

Pointer  to  second  argument  for  a  given  task 

Location  of  task  parameters  In  NARGTB 

Location  In  argument  list 

Points  to  last  task  executed  for  a  restart 
job 


LTRACE 


Trace  flag  for  debug 


NAHDAT 

NAH6E0 


NOSTAT 


NPRREC 


NUHARG 

NUMTSK 

NXTTSK 


TRACST 


YSSTAT 


I/O  VARIABLES: 
A.  INPUT 
CHKURT 
COHPLT 
OBGPRT 
ISOFF 


KBGEOH 


TSKXQT 


(MOM) 


Loop  Index 

User-assigned  name  of  geometry  data  set 

Pointer  to  default  geometry  data  set  name 
in  NCOOES 

Index  to  NCODES  array  for  the  task  name 
mnemonic 

.No  operation  flag 

Logical  flag  set  if  statistics  have  not 

been  requested 

Number  of  words  for  geometry  data  set 

record 

Hollerith  name  of  a  task 

Number  of  INTARG  arguments  for  a  task 

Task  Identification  number 

Pointer  to  the  next  task  to  be  executed 

Current  processor  time 

Loclcal  flag  set  if  trace  statistics  are 

desired 

Logical  flag  set  If  statistics  have  been 

requested 


LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 
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SB 

bbK 


s.s.V.VJ 
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(MOM) 


KBREAL 

KOLCNT 

KOLCOL 

KOLLOC 

KOLNAM 

KOLTIM 

KOLTSK 

KWOFF 

KWON 

KWSTAT 

KWTRAC 

LSTTPF 

LOOPMX 

LUPRNT 

MAXBLK 

MAXSEG 

MXARGS 

NAMTSK 

NARGTB 

NCOOES 

NDATBL 

NLOOPS 

NOGOFG 

NOPCOO 


/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/PARTAB/ 

/ADEBUG/ 

/SEGMNT/ 

/SEGMNT/ 

/ARGCOM/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/ADEBUG/ 

/ADEBUG/ 

/ADEBUG/ 


NOSTAT 


TSKXQT 


(HON) 


NPRSEG 
NPTASK 
NTINT 
NTSKTB 
NXTTSK 
RSTART 
RSTRTA 
B.  OUTPUT 
CHKPNT 
CHKWRT 
CPFRWO 
OBGPRT 
IERRF 
IHCCHK 
INTARG 
IOCKPT 
IPASS 
LSTTPF 
LTRACF 
MAXBLK 
NOGOFG 
HOSTAT 
NUMARG 
NXTTSK 
RSTART 


/SEGHNT/ 

/PARTAB/ 

/AOEBUG/ 

/PARTAB/ 

/AOEBUG/ 

/SYSFIL/ 

/SYSFIL/ 

LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/SYSFIL/ 

/AOEBUG/ 

/AOEBUG/ 

/SYSFIL/ 

/ARGCOM/ 

/SYSFIL/ 

/AOEBUG/ 

/SYSFIL/ 

/SYSFIL/ 

/SEGHNT/ 

/AOEBUG/ 

/AOEBUG/ 

/ARGCOM/ 

/AOEBUG/ 

/SYSFIL/ 


1236 


1.  NAME:  TSKXQT 


(OUTPUT) 


2*  PURPOSE:  To  read  the  task  list  and  call  the  appropriate  processor 
to  execute  the  tasks. 

3.  METHOD:  The  task  list  is  scanned  twice:  during  the  first  scan  the 
subroutines  necessary  to  execute  the  tasks  are  called  in  order  to 
initialize  the  required  parameters.  During  the  second  pass  the 
subroutines  are  called  to  perform  the  tasks  as  specified  by  the 
user. 

Task  execution  normally  begins  with  the  first  task  In  the  task  list 
and  proceeds  sequentially  through  the  list  unless  a  LABEL  task  Is 
encountered.  The  LABEL  task  will  redirect  execution  to  its  asso¬ 
ciated  LOOP  task  until  the  required  number  of  LOOP/LABEL  loops  has 
been  fulfilled.  Task  execution  terminates  when  an  END  command  Is 
encountered,  the  end  of  the  task  list  is  reached,  or  an  error  occurs 
In  executing  a  task. 

If  the  task  list  has  been  generated  by  a  RSTART  command,  execution 
nmy  not  necessarily  begin  at  the  top  of  the  task  list.  Normally, 
restart  Is  begun  from  the  task  which  wrote  the  checkpoint  read  in  to 
generate  the  task  list.  In  modules  subsequent  to  the  one  which 
generated  the  checkpoint,  execution  can  begin  at  the  top  of  the  task 
list  (If  the  preceedlng  run  did  not  complete  Its  execution)  or  at 
the  restart  task  (If  execution  was  successfully  completed). 

The  following  tasks  are  active  In  the  OUTPUT  module: 


FORTRAN 

TASK 

LABEL 

OUTPUT  MODULE  FUNCTION 

150 

CHKPNT 

Retrieve  timed  checkpoint  parameters  or 
write  a  command  checkpoint 

180 

DEBUG 

Turn  off  or  on  the  debug  flags 

200 

END 

Terminate  module  execution 

260 

LABEL 

Decrement  loop  counter  and  branch  to 
LOOP  If  positive 

270 

LOOP 

Initialize  loop  counter 

390 

RESTART 

Process  RSTART  command  error 

440 

WIPOUT 

Process  WIPOUT  command  error 

480 

GMDATA 

Advance  edition  of  geometry  data  set 

1239 


-■  *  '.v 1  v.'V •'-/  w, 

> 1  •  *  v 


TSKXQT 


(OUTPUT) 


530 

EFIELD 

Call  FLDDRV  to  print  and  plot 
and  Incident  fields 

scattered 

540 

DMP 

Call  DHPDRV  to  process  direct 
tions 

manlpula- 

570 

SETINT 

Call  SET  to  select  GTD, 
Incident  field  Interactions 

MOM,  and 

When  the  OUTPUT  module  execution  Is  complete,  the  alternate  restart 
flag  is  turned  off  and  the  pointer  to  the  last  task  executed  Is  set 
to  the  last  task  of  the  run. 


INTERNAL  VARIABLES: 

VARIABLE 

CPFRWD 

DBGPRT 

DT 

IBITS 

INCCHK 

INOXG 

IOCKPT 

ITASK 

KOUNT 

LINDX 

LOCARG 

LOCNXT 

L0CTP1 

LOCTSK 


DESCRIPTION 

Checkpoint  file  rewind  flag 
Debug  print  flag 

Time  interval  between  calls  to  TICHEK 

Attribute  word  for  geometry  data  set 

Checkpoint  time  Increment 

Pointer  to  symbol  table  entry  for  a  data 
set 

Logical  unit  number  of  checkpoint  file 

Pointer  to  task  In  task  list  being  executed 

The  number  of  times  the  loop  terminating  on 
the  reference  label  has  been  executed 

Index  to  the  loop  table  entry  currently 
being  executed 

Pointer  to  task  argument  In  NARGTB 

Pointer  to  NARGTB  for  the  next  task  to  be 
executed 

Pointer  to  first  argument  for  a  given  task 
Location  of  task  parameters  In  NARGTB 


.*•! 


A" 
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LSTARG 

LSTTPF 


LTRACE 


NAMDAT 

NAMGEO 


Location  in  argument  list 

Pointer  to  last  task  executed  for  a  restart 
job 

Trace  flag  for  debug 
Loop  Index 

User-assigned  name  of  geometry  data  set 

Pointer  to  default  geometry  data  set  name 
In  NCOOES 

Index  to  NCODES  array  for  the  task  name 
mnemonic 


NPRREC 


No  operation  flag 

Number  of  words  per  geometry  data  set 
record 


NOSTAT 


NUMARG 

NUMTSK 


Logical  flag  set  If  statistics  have  not 
been  requested 

Hollerith  name  of  a  task 

Number  of  INTARG  arguments  for  a  task 

Task  Identification  number 


NXTTSK 


Pointer  to  the  next  task  to  be  executed 


TRACST 


Current  processor  time 

Logical  flag  set  If  trace  statistics  are 
desired 


YSSTAT 


Logical  flag  set  If  statistics  have  been 
requested 


I/O  VARIABLES: 
I.  INPUT 
CHKURT 
COMPLT 
OBGPRT 


LOCATION 

/SYSFIL/ 

/SYSFIL/ 

/ADEBUG/ 
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A 


S 


ISOFF 

ISON 

KBGEOM 

KBREAL 

KOLCNT 

KOLCOL 

KOLNAM 

KOLTIH 

KOLTSK 

KWOFF 

KWON 

KWSTAT 

KWTRAC 

LOOPMX 

LSTTPF 

LUPRNT 

HAXBLK 

MAXSEG 

MXARGS 

NAMTSK 

NARGTB 

NCOOES 

NDATBL 

N LOOPS 

NOGOFG 


/ADEBUG/ 

/ADEBUG/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/SYSFIL/ 

/ADEBUG/ 

/SEGMNT/ 

/SEGMNT/ 

/ARGCOM/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 

/PARTAB/ 
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Tih 


/ADEBUG/ 


/ADEBUG/ 


/ADEBUG/ 


/SEGHNT/ 


/PARTAB/ 


/ADEBUG/ 


/PARTAB/ 


/ADEBUG/ 


/SYSFIL/ 


/SYSFIL/ 


LOCATION 


/SYSFIL/ 


/SYSFIL/ 


/SYSFIL/ 


/SYSFIL/ 


/ADEBUG/ 


/ADEBUG/ 


/SYSFIL/ 


/ARGCOH/ 


/SYSFIL/ 


/ADEBUG/ 


/SYSFIL/ 


/SEGHNT/ 


/ADEBUG/ 


/ADEBUG/ 
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NUMAR6 

HXTTSK 

RSTART 

RSTRTA 

6.  CALLING  ROUTINE: 
GEMACS 

7.  CALLED  ROUTINES: 
ASSIGN 

CONVRT 

DMPDRV 

ERROR 

FLDORV 

GETARG 

GETGEO 

SET 

STATIN 

STATOT 

SYHDEF 

SYHUPO 

SYSCHK 

TICHEK 

WLKBCK 

WRTCHK 


TSKXQT  (OUTPUT) 

/ARGCOM/ 

/ADEBUG/ 

/SYSFIL/ 

/SYSFIL/ 
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MISSION 
of 

Home  Air  Development  Center 

.  KAVC  plan*  and  executes  research,  development,  tut  and 
[  erected  acquisition  programs  in  support  oi  Command ,  Control 
i  Communication*  and  Intelligence  (C3IJ  activities.  Technical 
(  and  engineering  support  within  autos  oi  technical  competence 
u  provided  to  ESP  Program  OUices  f  PO-4  J  and  other  ESP 
elements .  The  principal  technical  mission  alien*  ate 
anmam'mtions,  electromagnetic  guidance  and  control ,  sur- 
vtAUance  oi  ground  and  aerospace  objects ,  intelligence  data 
collection  and  handling ,  iniormation  system  technology , 
ionospheric  propagation ,  solid  state  sciences,  microwave 
physics  and  electronic  reliability ,  maintainability  and 
compatibility. 
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